AD-A008  035 


A HIGH  ALTITUDE  INFRARED  RADIANCE  MODEL 
Thomas  C.  Degges 
Visidyne,  Incorporated 


Prepared  for: 

Air  Force  Cambridge  Research  Laboratories 
Defense  Advanced  Research  Projects 


27  December  1974 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (Whit  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


• REPORT  NUMBER 

AFCRL-TR  -74-0606 


4*  TITLE  (and  Subtitle) 


A HIGH  ALTITUDE  INFRARED  RADIANCE 
MODEL 


t.  author*-*;  


Thomas  C.  Degges 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


*•  RECIPIENT'S  CATALOG  NUMBER 

/3  5 


5.  TYPE  OF  REPORT  A PERIOD  COVERED 

Final,  1 June  1972 
through  30  April  1974 


6 PERFORMING  ORG.  REPORT  NUMBER' 

VI-236 


S.  CONTRACT  OR  GRANT  NUMBERCa; 

F19628-72-C-0330 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Visidyne,  Inc. 

19  Third  Avenue,  N.  W.  Industrial  Park 

r-> li i _ % « i n 


II  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Cambridge  Research  Laboratories 
Hanscom  AFB,  Massachusetts  01731 
Contract  Monitor:  Bertram  D.  Schurin/OPI 

~4  MONITORING  AGENCY  NAME  S AOORESSCIf  dlllerent  Irom  Controlling  Ollice) 


16.  DISTRIBUTION  STATEMENT  (of  I hi.  R.porl) 


io.  program  element,  project,  task 

AREA  6 WORK  UNIT  NUMBERS 


8692  N/A  N/A 


12.  REPORT  DATE 

December  27,  1974 


<3.  NUMBER  OF  PAGES 

346 


IS.  SECURITY  CLASS,  (ot  this  report) 


Unclassified 


15«.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited 


IT.  DISTRIBUTION  STATEMENT  (o I the  abstract  entered  In  Block  20.  If  different  from  Report) 


18  SUPPLEMENTARY  NOTES 


Reproduced  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

U S Department  ol  Commerce 
Springfield,  VA.  22151 


PRICES  SUBJECT  10  CHANGE 


This  research  was  sponsored  by  Defense  Advanced  Research 
Projects  Agency.  ARPA  Order  No.  1366 


19.  KEY  WOROS  (Continue  on  r verse  aide  if  necessary  and  Identify  by  block  number ) 

Atmosphere,  Infrared,  Radiance,  Model 


20.  ABSTRACT  (Continue  on  revered  tide  II  nece.tery  end  Identity  by  block  number ) 

A knowledge  of  the  natural  infrared  radiance  originating  in  the 
earth's  upper  atmosphere  is  of  interest  for  systems  design,  military 
surveillance  and  the  advancement  of  knowledge  about  physical  pro- 
cesses in  the  upper  atmosphere.  A physical  model  that  includes 
experimental  data  on  and  theoretical  estimates  of  excitation  processes 
that  lead  to  emission  of  infrared  radiation  has  been  implemented  in  a 


EDITION  OF  t NOV  65  IS  OBSOLETE 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  iDw  D.r.  Entered) 


Unclassified 


FC.UhlTV  CLASSIFICATION  OF  THIS  PAt  Dalm  Kntsnwfl 


20.  computer  program  that  computes  infrared  radiances  for  an  earth 
limb  viewing  geometry.  The  nominal  spectral  region  of  this  study 
lies  between  2.7  and  25  micrometers  and  emphasis  is  placed  on 
radiation  originating  at  altitudes  between  70  and  500  km.  The  phys.-al 
model  is  described,  with  emphasis  on  the  changes  required  in  ex- 
tending its  usefulness.  Application  of  the  computer  program  is  de- 
scribed and  estimates  are  given  of  uncertainties  in  results  due  to 
assumptions  made  in  the  model  and  lack  of  data  on  actual  atmospheric 
composition. 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfRTiMi  Data  Hnlmrnl) 


TABLE  OF  CONTENTS 


SECTION 

I 

H 

III 

IV 

V 

APPENDIX 


AIR  FORCE  Cl)  FEBRUARY  25, 


Page  No. 


Introduction  1 

Radiative  Processes  and  Modelling  5 

Atmospheric  Composition  39 

Results  from  Computational  Model  55 

High  Altitude  Radiance  Programs  77 

References  261 

Chemistry  Program  SN\PS2  273 


1975  — 2 11 


iii 


[ 

t 


I 


SECTION  I 


INTRODUCTION 


The  research  reported  here  has  had  as  its  major  objective 
the  development  of  a computer  program  to  simulate  the  natural 
infrared  radiance  background  of  the  earth's  upper  atmosphere. 

The  nominal  spectral  region  under  study  lies  between  2.  7 and  25 
micrometers  and  emphasis  is  placed  on  radiation  originating  at 
altitudes  between  70  and  500  kilometers.  The  general  problem 
area  is  of  interest  for  systems  design,  military  surveillance  and 
the  advancement  of  knowledge  about  physical  processes  in  the  upper 
atmosphere.  The  immediate  application  of  this  work  will  be  to  aid 
in  developing  optimum  infrared  background  measurements  programs 
and  in  interpreting  the  results  of  such  measurements. 

This  work  is  an  extension  of  the  study  of  Corbin,  et  al.  (1969) 
and  Degges  (1972).  The  former  investigated  the  natural  infrared 
background  of  the  earth  in  the  5 to  25  micrometer  spectral  region, 
with  the  goal  of  estimating  earth  limb  viewing  radiances  for  tangent 
heights  from  the  surface  to  500  km  altitude.  For  convenience,  their 
study  divided  the  atmosphere  into  two  regions  with  a division  at  70  km. 
Below  70  km  the  atmosphere  wj  s assumed  to  be  in  thermal  equilibrium. 
Above  70  km  explicit  calculations  were  made  of  processes  which 
excite  and  de- excite  molecular  vibrational  and  rotational  levels  which 
are  the  source  of  infrared  radiation.  Their  study  concentrated  on 
radiation  from  water  vapor,  carbon  dioxide,  ozone,  nitric  oxide  and 
nitrous  oxide,  which  are  the  principal  radiating  species  in  the 
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spectral  region  considered.  In  addition,  nitric  acid  was  included 
in  the  lower  atmosphere  work  and  estimates  were  made  of  radiation 
to  be  expected  from  particulate  matter  suspended  in  the  atmosphere. 

Corbin,  et  al.  (1969)  presented  models  for  the  lower  atmos- 
phere for  a wide  range  of  seasonal  and  latitudinal  conditions.  This 
was  not  possible  for  the  abundances  of  most  minor  neutral  species. 

More  data  has  since  become  available,  particularly  for  nitric  oxide 
and  the  hydroxyl  radicil,  but  at  present  it  appears  that  the  best 
means  of  estimating  abundances  of  important  infrared  emitting 
species  is  chemical  rate  equation  integrations  including  molecular 
diffusion  and  eddy  mixing.  Even  calculations  involving  transport 
properties  cannot  always  be  accepted  because  the  values  of  eddy 
mixing  coefficients  are  to  a large  extent  only  informed  guesses  and 
published  calculations  often  use  outdated  rate  coefficients.  Degges 
(1972)  reported  a computer  program  with  which  to  determine  diurnal 
variations  in  abundances  of  minor  species,  to  investigate  the  effects 
of  changes  in  assumed  eddy  mixing  coefficients,  and  to  estimate  the 
effects  of  new  determinations  of  chemical  rate  coefficients. 

A second  area  of  study  involved  in  improving  the  radiance  model 
includes  the  physical  processes  that  control  the  population  of  infrared 
emitting  states  of  atmospheric  molecules.  Except  for  the  pure  ro- 
tational radiation  from  molecules  such  as  water,  the  degree  of  excitation 
of  vibrational  levels  determines  the  radiation  from  infrared  emitting 
molecules.  The  most  important  mechanisms  are  collisional  exci- 
tation and  de-excitation  and  absorption  and  re-emission  of  electro- 
magnetic radiation. 

Jn  the  troposphere  and  lower  stratosphere,  collisional  processes 
are  rapid  enough  to  control  the  population  of  vibrational  levels. 
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f bove  30  to  50  km,  however,  collisional  excitation  becomes  less 
efficient  and  radiative  processes  become  important.  The  combined 
effects  of  collisional  and  radia.tive  processes  must  therefore  be  con- 
sidered, Below  an  altitude  of  about  90  km  molecular  nitrogen  and 
oxygen  are  the  most  important  collision  partners.  Above  that 
altitude,  atomic  oxygen  becomes  important,  both  in  exciting  nitric 
oxide  and  in  determining  molecular  oxygen  and  nitrogen  vibrational 
temperatu  res. 

In  determining  the  effects  of  radiation  on  the  populations  of 
vibrational  levels  of  infrared  active  molecules,  it  is  necessary  to 
separate  the  radiation  of  a single  change  in  vibrational  quantum 
numbers  from  the  rest  of  the  radiation  field.  The  previously  re- 
ported studies  did  this  by  assuming  a Doppler  line  shape  for  the 
individual  rotational  lines  of  a band  and  were  able  to  obtain  adequate 
numerical  approximations  for  radiative  transfer  functions  appropriate 
to  single  bands  of  linear  molecules,  and  less  accurately,  for  water 
vapor  and  ozone  bands  lying  in  the  spectral  region  of  interest. 

The  principal  areas  of  research  and  program  revision  reported 
here  include  a further  study  of  band  radiance  modeling  and  the 
effects  of  line  shape  on  radiance  computations.  In  addition  to  some 
exploratory  calculations  using  a Voigt  line  profile  instead  of  a 
Doppler  line  shape,  some  asymptotic  analytical  results  have  been 
obtained  which  provide  an  estimate  of  the  accuracy  of  the  numerically 
obtained  band  functions  at  large  optical  thicknesses.  It  is  found  that 
when  a Voigt  profile  is  used,  the  connection  between  the  band  models 
and  functions  for  a single  line  is  more  readily  seen  than  when  a 
Doppler  line  shape  is  used.  This  material,  and  a discussion  of 
changes  made  in  the  numerical  methods  is  described  in  Section  II. 
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The  status  of  current  knowledge  of  atmospheric  composition 
is  reviewed  in  Section  III.  Recent  measurements  of  chemical  rats 
coefficients  have  resulted  in  changes  in  estimates  of  abundances 
of  some  molecular  species.  An  attempt  is  made  to  set  upper  and 
lower  bounds  for  mixing  ratios  at  high  altitudes. 

Results  of  calculations  are  given  in  Section  IV,  and  estimates 
are  made  of  the  validity  of  the  results. 

The  computer  programs  are  described  in  Section  V. 


SECTION  II 


RADIATIVE  PROCESSES  AND  MODELLING 
A.  INTRODUCTION 

The  physical  processes  and  assumptions  used  in  developing  a 
model  for  transport  of  infrared  radiation  in  the  upper  atmosphere 
are  discussed  below,  beginning  with  a general  overview  and 
discussion  of  a simpler  problem.  Next,  the  interaction  of  electro- 
magnetic radiation  with  an  isolated  atom  or  molecule  is  described. 

An  understanding  of  the  basic  processes  involved  is  the  first  step 
in  any  computation  of  t'  e radiation  from  a gas  or  mixture  of  gases. 
Subsequently,  absorption  from  a continuum  source,  a spectrally 
varying  source  and  fluorescent  processes  are  treated.  The  effects 
of  optical  thickness  when  the  path  of  the  radiation  includes  a large 
number  of  molecules  of  the  same  species  is  considered  for  spectral 
lines  with  a Voigt  profile,  which  includes  Doppler  (Gaussian)  and 
Lorentz  shapes  as  limiting  cases.  Finally,  a model  for  the  transport 
of  radiation  for  a single  vibration- rotation  band  is  described,  and  the 
effects  of  temperature  estimated. 

Figure  II- 1 illustrates  the  geometry  employed  in  the  current 
modeling  study.  A molecule  at  point  P absorbs  radiation  from  three 
sources;  (1)  a collimated  beam  from  the  sun,  (2)  the  lower  boundary 
which  radiates  as  a black  body,  and  (3)  other  molecules  between  the 
lower  and  upper  boundaries.  Radiative  transfer  functions  are  needed 
to  compute  the  effects  of  absorption  of  radiation  by  other  molecules 
between  each  of  the  sources  and  the  molecule  at  point  P. 

In  the  absence  of  collisional  excitation  and  de- excitation,  the 
simplest  integral  equation  that  may  be  written  to  describe  radiative 
transport  in  a finite  plane  parallel  atmosphere  is 

rT 

. I max 

J(T)=2  / K1(|T  - 1 1)  J(t)dt  + S(T)  (2-1) 

J o 
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BLACK  BODY  SURFACE,  TOP  OF  LOWER 
ATMOSPHERE 

Geometry  assumed  in  deriving  the  radiative  transfer  functions 


FIGURE  II- 1 
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This  equation  assumes  that  there  is  only  one  radiative  transition 
involved  and  that  the  radiation  emitted  by  a molecule  after  absorption 
is  independent  of  the  direction  of  incidence.  Here  V is  a measure  of 
the  optical  path  length  in  the  vertical  direction.  J(t)  is  a dimension- 
less measure  of  the  population  of  the  upper,  emitting  molecular  state. 
K1  is  radiative  transfer  function  which  accounts  for  absorp- 

tion of  radiation  emitted  by  molecules  at  a vertical  path  length  T 
distant  from  the  absorbing  molecule.  S(r  ) is  an  auxiliary  function 
which  includes  the  effects  of  radiation  originating  outside  the  atmos- 
phere. If  there  is  only  a black  body  source  located  at  the  lower 
boundary  from  which  optical  path  length  is  measured  ( T = 0 at  the 
lower  boundary)  then 

S(r)  = 2 No  K2  (r>  (2-2 

Here  is  the  radiance  of  the  lower  boundary  in  arbitrary  units 
such  as  watts  cm  - ster  or  the  spectral  radiance  in  units  such 

,,  -2-2  -1  i 

as  watts  cm  cm  -ster  (wavelength  interval) " . 

In  the  gray  or  spectrally  independent  case,  when  equation  (2-2) 
holds,  the  kernel  function  (X  ) is  the  first  exponential  integral 
Ej  (T)  and  the  function  describing  absorption  from  the  black  body 
source  is  the  second  exponential  integral  Ez  ('f  ).  These  functions 
are  described  in  detail  in  such  works  as  Kourgar.off  (1963).  A 
complete  analytic  solution  is  available  and  the  auxiliary  functions 
required  for  a numerical  solution  have  been  tabulated  by  King  (1956). 
Figure  (2-2)  plots  the  solution  of  equation  (2-1)  when  equation  (2-2) 
holds  for  the  gray  case  for  values  of  total  optical  thickness  T 

max 

0,  0.  1,  1.0,  and  10.  When  the  total  optical  thickness  is  small,  the 


effective  excitation  of  the  molecule  at  any  point  in  the  atmosphere 

is  half  that  of  a molecule  at  the  (vibrational)  temperature  of  the 

black  body.  As  the  total  optical  thickness  increases,  the  excitation 

of  the  molecules  near  the  black  body  increases,  while  that  of 

molecules  near  the  top  of  the  atmosphere  decreases.  This  occurs 

because  as  the  optical  thickness  increases,  molecules  near  the 

lower  boundary  absorb  an  increasing  amount  of  radiation  re-radiated 

from  molecules  above  them.  At  the  same  time,  less  radiat'on 

penetrates  to  higher  altitude.  The  physics  >f  the  problem  require 

that  if  a second  black  body  of  the  same  temperature  be  placed  at 

the  top  of  the  atmopshere,  the  population  of  radiators  at  all  points 

would  be  uniformly  that  corresponding  to  the  black  body  temperature. 

This  symmetry  and  the  linear  character  of  the  integral  equation  (2-1) 

results  in  the  relation  J(T)  + J(T  -T)  = N 

max  o' 

When  a collimated  beam  of  light  is  incident  on  the  top  of  the  atmosphere, 
another  term  must  be  added  to  the  function  S(T)  in  equations  (2-1)  and 
(2-2)  to  account  for  this.  If  the  flux  is  incident  at  a zenith  angle  with 
cosine/*and  has  an  intensity  in  arbitrary  units  such  asTTF  watts  cm"2 
normal  to  its  direction,  the  term  added  to  S('t')  has  the  form 

FMl[(Tmax  " W/*]‘  In  the  gray  case,  M^T)  is  simply  exp(-T). 

The  corresponding  functions  for  individual  spectral  lines  and  bands 
will  be  described  after  a closer  look  at  the  interaction  of  electro- 
magnetic radiation  with  a single  molecule. 

B.  INTERACTION  OF  RADIATION  WITH  A SINGLE  MOLECULE 

When  treating  a problem  in  which  the  number  of  photons  emitted 
or  absorbed  by  an  atom  or  molecule  must  be  considered  rather  than 
total  energy  emitted,  it  is  convenient  to  examine  the  physics  of  the 
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interaction  of  radiation  with  matter  through  use  of  the  Einstein  A and 
B coefficients.  In  treating  a radiative  problem  such  as  the  current 
study  of  high  altitude  infrared  radiances  which  also  involves  transfer 


of  energy  between  molecules  through  collisions,  this  approach  is  a 
necessity.  In  the  follow,  rj,  discussion,  cgs  units  are  used  except 
where  explicitly  stated. 

The  Einstein  A coefficient  is  the  rate  at  which  an  atom  or 
molecule  in  an  excited  state  spontaneously  emits  radiation.  It  has 
the  units  photons / sec -molecule.  In  an  assemblage  of  molecules  all 
in  a state  which  can  radiate,  the  average  rate  of  radiation  is  given 
by  the  Einstein  A coefficient  for  the  particular  radiative  transition 
from  a state  of  higher  energy  (subscript  u)  to  one  of  lower  energy 
(subscript  1). 

In  addition  to  spontaneous  emission,  two  other  processes  occur. 

These  are  absorption,  characterized  by  the  Einstein  coefficient  for 

absorption  B^,  and  induced  emission,  characterized  by  the  Einstein 

coefficient  for  induced  emission  B^.  In  a black  body  radiation  field 

with  spectral  density,  yOi  , the  probability  that  an  atom  or  molecule 

lu 

absorbs  a quantum  of  radiation  in  unit  time  is  B,  a . . When  radiation 

lu  rj 

lu 

density  is  expressed  in  c.  g.  s.  units,  ergs- sec/cm3,  the  coefficient 
3 

B[u  has  units  cm  /ergs-sec.  Similarly,  an  atom  or  molecule  in  a 

radiation  field  with  the  above  spectral  density  and  in  an  excited  state 

emits  radiation  at  the  rate  B ^ in  addition  to  that  characterized 

ulrV. 

lu 

by  the  coefficient  A . The  Einstein  coefficients  are  properties  of 

ul 

an  atom  or  molecule.  If  one  coefficient  is  known,  the  others  may  be 
determined  by  the  relations 

3 


A = 87ThB  /A' 

ul  ul 


(2-3 


and 
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Thus,  it  is  possible  to  treat  the  problem  knowing  only  the 

Einstein  A , coefficient, 
ul 

Experimental  determinations  of  the  interaction  of  infrared 
radiation  with  gases  are  usually  made  by  measurements  of  light 
absorption.  Experimental  results  may  be  expressed  as 
integrated  absorption  from  a continuum  source. 


Siu  = N1 


A A2  g 
ul  lu  5u 
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or  as  line  or  band  strengths 


— 

/ 

1 - exp 

he 

A.  kT 
!u  j 

(2-5) 


lu 


sA„/pc 


(2-6) 


where  p is  the  pressure. 

In  the  visible  and  ultraviolet  regions  of  the  spectrum,  it  is 
customary  to  use  oscillator  strengths  to  characterize  the  interaction 
of  radiation  and  matter.  In  this  case,  the  defining  equation  is 


TTe 


lu 


me 


Ki£i„ 


/ 


1 - exp 


he 

A,  kT 

lu 


A 


647T 


g, 


ul 


3h 


A3  g 

' u 


(2-7) 


The  results  of  theoretical  computations  of  the  interaction  of 
matter  with  radiation  are  frequently  expressed  as  dipole  moments, 
//ul-  These  are  related  to  the  Einstein  A coefficient  by 


ul 


(2-8) 


The  rate  at  which  a single  atom  or  molecule  absorbs  and 


reradiates  light  from  a parallel  beam  from  a continuum  source  at 
a wavelength  ^ corresponding  to  a resonance  transition  is 
(Mitchell  and  Zemansky,  1961): 


TT  e 


me 


I f 
o o 


(photons/sec) 


= 8.  852  x 10"13I  f ^ 1 
o o o 


(photons/sec) 


(2-9) 


where  I is  the  continuum  source  flux  and  f is  the  oscillator 
o ° 2 

strength.  If  the  incident  flux  is  expressed  in  photons/cm  -^im, 

and  oscillator  strength  is  replaced  by  the  equivalent  line  or  band 

-2  - 1 

strength  S expressed  as  cm  atm  , then 

G - 2.  72  x 10"24I  sA2 
o 


(2-10) 


Implicit  in  this  expression  is  the  convenient  relation  between  the 
Einstein  A coefficient  and  line  or  band  strength  with  wavelenth 
dimensions  yum 


ul 


= 2.  804S(cm  2atm  ^ 


)//\2(yu.  m) 


(2-11) 


When  the  source  is  an  extended  black  body  surface  with  a 
spectral  radiance  , the  G factor  is  obtained  by  integration 
over  the  solid  angle  subtended  by  the  source.  For  radiation  from 
a plane  parallel  surface  of  infinite  extent, 

G = 1.  169  x 10'23N^  sA2  (2-12) 


Another  purely  radiative  mechnaism  of  importance  in  upper 
atmosphere  radiance  computations  is  fluorescence.  A striking  example 


occurs  for  the  weak  10.  4 fXm  CC>2  band.  Absorption  of  solar 
radiation  by  CO^  4.  3 yl(m  band  maintains  the  vibrational 
mode  at  a vibrational  temperature  of  about  290  K.  During  the 
day,  this  leads  to  an  enhancement  of  emission  in  the  10.  4 yUm 
band  by  as  much  as  two  orders  of  magnitude  when  viewed  along 
the  earth's  limb  from  outside  the  atmosphere. 

Fluorescence  occurs  when  the  upper  state  of tne  resonance 
transition  reradiates  into  other  than  the  initial  states.  The  relative 
probability  of  a radiative,  transition  to  a lower  state,  n from  the 
upper  state,  u when  there  are  k possible  lower  state  is: 


P 

un 


A 

un 

k 


1 


A . 
uj 


j = 1 


where  the  A^  are  the  Einstein  A coefficients  for  the  spontaneous 
emission  probabilities  for  a transition  from  the  upper  state  to  the 
jth  lower  state. 

In  the  case  of  molecular  resonance-fluorescent  scattering  at 
low  temperature,  only  the  lowest  vibrational  level  (v"  = 0)  in  the 
ground  electronic  state  is  populated  significantly.  Absorption  of  a 
photon  excites  the  molecule  to  a vibrational  level,  v'  in  the  same 
or  a higher  electronic  state.  Radiation  is  possible  from  this  state 
to  the  initial  level  (v"  £ 0)  or  to  other  vibrational  levels  (v"  = 0)  in 
the  ground  electronic  state.  The  late  of  resonance  fluorescent 
scattering  in  a given  electronic  vibrational  transition  may  be 
obtained  by  combining  Equations  (2-9)  and  (2-13)  and  introducing 
appropriate  subscripts: 
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-13) 
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(2-14) 
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1 , „ = 8.  853  x 10  I f 
v'v"  v'o  v o v'o 


Here  the  subscript  v'v"  refers  to  a transition  from  an  upper  state 
with  vibrational  level  v'  to  a lower  state  with  vibrational  level  v". 

The  subscript  v'o  refers  to  absorption  of  a photon  by  a molecule  in 
the  lowest  vibrational  level  of  the  ground  state,  exciting  the  molecule 
to  a vibrational  level  v1  in  the  same  or  a higher  electronic  state. 


C.  RADIATIVE  TRANSFER  FUNCTIONS  FOR  LINES  AND  BANDS 

In  developing  a realistic  model  of  the  interaction  of  radiation 
with  matter,  it  is  necessary  to  take  into  account  the  details  of 
absorption.  Line  shape  must  be  considered  because  it  determines 
the  amount  of  radiation  excaping  from  an  optically  thick  layer  of  gas. 
Thus,  for  large  optical  thickness  (greater  than  a magnitude  of  about 
10),  the  total  emission  from  an  isolated  spectral  line  with  Lorentz 
shape  is  proportional  to  the  square  root  of  the  optical  thickness, 
while  that  from  a line  with  Doppler  shape  is  proportional  to  the 
square  root  of  the  logarithm  of  the  optical  thickness.  In  addition, 
absorption  of  radiation  emitted  at  one  place  in  the  gas  by  a molecule 
at  another  location  cannot  be  computed  without  knowledge  of  the  line 
shape. 

At  high  altitudes,  the  profile  of  the  spectral  absorption  curve 
of  individual  line  in  a vibration-rotation  band  is  predominantly 
due  to  velocity  or  Doppler  broadening.  Kuhn  and  London  (1969),  who 
investigated  the  heat  budget  of  the  atmosphere  between  30  and  110  km, 
report  that  radiative  transfer  calculations  require  the  use  of  the 
Voigt  profile  (including  both  collisional  and  Doppler  broadening)  at 
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altitudes  between  20  and  70  km.  Alternatively,  the  Doppler  profile 
is  adequate  at  higher  altitudes  where  the  Doppler  shape  is  the 
largest  contributor  to  line  broadening.  The  original  computational 
model  therefore  assumed  a Doppler  profile. 

The  discussion  below  assumes  a Voigt  profile  for  a single  line, 
and  it  will  be  shown  that  for  small  optical  thicknesses  and  low 
pressures,  the  difference  between  calculations  with  Doppler  and 
Voigt  profiles  is  small.  Asymptotic  values  for  l.irge  optical 
thicknesses  are  quite  different,  however. 

The  radiative  transfer  functions  required  to  treat  transfer 
of  resonance  radiation  in  a single  line  are  well  known  for  the 
Voigt  profile  and  its  two  limiting  cases,  the  Doppler  and  Lorentz 
profiles.  Avrett  and  Hummer  ( i ? A 5 ) review  the  properties  of 
four  of  these  functions,  giving  series  expansions  for  small  values 
of  optical  thickness  and  asymptotic  expressions  for  large  values 
of  the  optical  thickness.  The  notation  used  for  these  functions  varies 
widely.  That  used  here  is  adapted  from  that  used  by  Ivanov  and 
Shcherbakov  (1965  a,  b)  in  their  tabulation  of  functions  for  the 
Doppler  profile. 

The  best  known  of  these  functions  is  probably  chat  relating 
to  the  equivalent  width  of  a line  in  absorption  or  to  the  total  radiation 
emitted  in  a line  from  a source  in  thermal  equilibrium  at  a constant 
temperature,  the  integrated  absorbtance,  defined  as 


L(a.  T) 


Here  f(a,  x)  is  used  to  denote  the  Voigt  profile, 


(2-15) 
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f(a,  x)  = 


..  .2  2 

(t-x)  + a 


The  parameter  a is  proportional  to  the  ratio  of  the  Lorentz  half 
width  at  half  maximum  b and  the  Doppler  half  width  at  half 

i-J 

maximum 


(ln2) 


The  parameter  x is  proportional  to  the  difference  between  the 
spectral  frequency  of  the  line  center  and  that  at  any  other 
point  in  the  line  profile,  l) 

_ 2(\)-^)0)  1/2 

x (ln2) 

bD 

The  quantityTis  the  optical  thickness  at  the  center  of  a line  with 
Doppler  breading  alone. 


r = - fS.  {lnZ)m 

me  bD 


where  N is  the  number  density  of  ground  state  molecules  and  1 is 
the  path  length.  The  actual  optical  thickness  at  the  center  of  a line 

2 

with  Voigt  profile  is  the  product  of  theTof  Equation  (2-19)  and  exp  (a  ) 
times  the  co-error  function  of  the  parameter  a,  erfc(a).  In  general, 
Equation  (2-15'.  must  be  evaluated  numerically.  For  small  values 
of  T,  a series  expansion  may  be  obtained  in  the  form 


(2-20) 
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where 


$n(a)  = v M 


oo 


[f(a,  x)]ndx 


(2-21) 


00 


For  all  values  of  a,  (a)  is  one.  When  (^^(a)  is  divided  by  erfc  a, 
a smooth  variation  between  a = o (Doppler  profile)  and  a = qo 
(Lorentz  profile)  is  obtained  as  f hown  in  Figure  (II- 3)  for  small 
n.  For  large  values  of  T*,  asymptotic  series  are  available,  as 
given  for  example  by  van  Trigc  (1968). 


L(o,r)~2[ 


1/2 


, . .2886 

1 + — r-=r  +.  . . 


In  T 


(2 -22a) 


L(a*o,T) 


2 (a 


(7 r) 


r)1/4  L 8 aT  J 


(2- 22b) 


Equation  (2-10)  gives  the  rate  at  which  a single  molecule 
absorbs  radiation  from  a parallel  beam.  If  the  molecule  is  located 
within  an  optically  thick  medium,  the  rate  of  absorption  of 
radiation  is  reduced  by  a factor  which  depends  on  the  optical  thick- 
ness, 

roo 


Mi(a,T)  = ;?r/ 


f(a, x)e 


- Tf  (a,  x) 


dx 


(2-23) 


J-c 


-Oo 


For  small  values  of?*,  a series  expansion  is 


CO 


.n 


M 


i(a’T)'Z  ^fJ-<l(a> 

n+1 


(2-24) 


n=o 


Asymptotic  expressions  for  large  values  of  are 


Mi<0’T)~  Wf^r 


[ 


1 - 


0. 2886 


InT 


• • • • 


] 


(2-25a) 
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FIGURE  II- 


<2-25b> 

The  function  is  the  deriviative  with  respect  to  "?*of  the  function 
L.  In  addition  to  its  use  in  modifying  the  G factor  of  Equat 
(2-10),  it  gives  the  probability  of  escape  of  a photon  in  a direction 
with  optical  thickness  TT. 

A third  function  is  defined  as 


r°° 

M2  (a’T)  = W/ 

'-OO 


..  2 - Tf(a,  x) 

f(a,x)  e dx 


(2-26) 


and  has  a series  expansion 


M2(a 


OO 

■ r>  ■ I *3*- § »«<*> 

n=o 


(2-27) 


This  function  is  proportional  to  the  rate  at  which  a molecule 
absorbs  radiation  from  another  molecule  at  a optical  distanced. 
This  function  would  be  used  as  the  kernel  function  in  Equation 
(2-1)  if  one  dimensional  radiative  transport  were  being  treated. 

It  is  included  here  because  of  its  use  in  generating  the  kernel 
function  for  a plane  parallel  geometry. 

When  the  source  of  radiation  is  an  infinitely  extended  plane 
black  body,  the  function  which  is  used  to  determine  the  rate  at 
which  a molecule  absorbs  radiation  is  given  by 


N12(a,T) 


(2-28a) 


(2-28b) 
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(2-28c) 


4 |cos  & |S” 


6 d$ 


The  aeries  expansion  for  this  function  is  given  by 


I12(a,T)  =<|1(a)  + (YT-T+  rinT)<f>2(a) 


-I2 

n=2  (n-l)n!  I 


Here  a new  integral  is  defined. 


w/ 


[f (a,  x)]  In  f(a,  x)dx 


and  ^ is  the  Euler-Mascheroni  constant  0.  5772.  . . The  asymptote 


expansions 


N12(o’T)SVff  <0-5-2^+----) 


N12W 


may  be  obtained  from  Equations  (2-29b)  and  (2-25). 

The  kernel  functn  or  a plane  parallel  geometry  is  given  by 

r 

N2!<a'  = y=p=/  ff(a'x)]2  EjTffa.x)]  dx 

-oo 


(2-3la) 


dt 

M2(t)  -p 


r 

. 7772 


[ X sin  <9  d$ 
2 cos$  cos$ 


The  series  expansion  useful  for  small  T is 


oo 

n=  1 n*  n!  X ] 


ani  the  asymptotic  expansions  are 


(2  - 3 lb ) 


(2-31c) 


(2-32) 


n21(o.t)s;; 


rn  r 0193 
( • - lnT  + ) 


(2-33a) 


a1/2  /l 

N2!(a^°'  77‘l/4’X'3/2  (T  + ] (2-33b) 


The  recurrence  relation 


Nkn(a.TI  = 


*^r  [Mk<a’ 


T)-TNk+i,  „-!<*•  T> 


(2-34 


carries  over  into  the  corresponding  band  functions  and  provides 
a partial  check  on  numerical  calculations. 

In  extending  this  treatment  of  radiative  transfer  tc  a vibration- 
rotation  band,  in  addition  to  assuming  that  each  line  has  the  same 
Voigt  profile  shape,  the  following  assumptions  are  made: 

(1)  Rotational  levels  within  a vibrational  level  remain  in 
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equilibrium  with  the  translational  temperature. 

This  is  valid  up  to  at  least  150  km  according  to 
Goody  (1964). 

(2)  There  is  complete  redistribution  of  frequency  within 
a single  line  following  absorption  of  radiation. 

(3)  There  is  complete  redistribution  of  absorbed  energy 
among  all  rotational  levels. 

(4)  Line  strengths  are  used  as  if  each  line  were  at  the 
band  center,  neglecting  wavelength  variation  across 
a band. 

(5)  There  is  no  overlap  of  adjacent  lines. 

(6)  Temperature  variations  at  different  levels  do  not 
affect  the  radiative  transfer. 

Let  g be  the  fraction  of  total  band  strength  associated  with  a 

J 

single  line,  so  that 


j ei=i 


(2-35) 


and  the  sum  over  J may  involve  complete  summation  over  rotational 
quantum  number  for  two  or  more  branches.  The  optical  thickness 


at  the  center  of  a single  line  will  be  given  byTg  , where  is  the 

J 


optical  thickness  that  would  result  if  the  entire  band  strength  were 
in  a single  line.  Then  the  radiative  transfer  functions  for  a band 
may  be  defined  as  functions  of  radiative  transfer  functions  for  a 
single  line.  First, 


to  \ 

L (a,  T ) = / L(a,  Tg  T) 
J J 


where  the  superscript  b refers  to  an  entire  band. 
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(2-36) 


In  order  that  the  correct  derivative  relations 


Mib(a>r)  =7T  Lb(a*  r) 


M2  (a ,T)  = 


Tr  Mib(a’T) 


and  the  integral  relations  of  Equations  (2 -28b)  and  (2 -3 lb)  hold. 


Mr(a’*n  = Z gJ 

J 


M2b(a’T)=  Z gj  M2(a’^gj) 

j 


N12  (a,  r 


•I 


gJ  N12(a’  TgJJ 


2 


N21(a’r'=  X*3  Nzi^'Tgj) 

J 


While  these  functions  must  be  evaluated  numerically,  it  is 
poHdible  to  perform  this  evaluation  in  two  different  ways  for 
small  values  ofT.  Making  the  definitions 


!n  = X *J 


G 


L 
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(2-37a) 


(2-37b) 


(2-38a) 


(2- 38b) 


(2-38c) 


(2-38d) 


(2-39a) 


(2-3  9b ) 
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The  following  series  expansions  are  obtained: 


(2-40a) 


Mj(a,r) 


(2-40b) 


(2-40c) 


Ljta.TI^G^^alT  + ~ 


G2inr  + gl  + y - y 


$2(a» 


+ 02  f L(a) 


2 

r + 


oo 

5.  <Hi£  Gn  $„<*> 

n=3 


(2-40d) 


N^a.T)  = G^'  j(a)  + G2  ([  2(a)Tln T + f GL(J>2( 


(a) 


+ (V-1)g2  C]  2(a)  + GLCf>L(a) 


r - 


-I 

n — > 


irjri 


n 


(n-l)n! 


n+1 


\+i(a» 


(2-40e) 
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(2-41a) 


Parallel 

Lb(o ,T)  - V TF€q  lnT 

x\ 

Mb(o,T)  V7TTqr/T 
Mb(o,T)  - VWQr/T2 
Lb(o,T)~  ~V7T€Qr  lnT 

N b2  (°,  T)  -*  4 V 7T  Cqr  / r 

Nb  j (o,  r)  - jvw Qr  /r2 


Perpendicular 

-f  V 7TCQ  lnT 
2 R 


fVTTFcyT 


7 V TTfQj,  lnT 


-VWq  /t 

4 R7  C 


A VWq  /r2 

4 R 


(2-41b) 


(2-41c) 


(2-41d) 


(2-41e) 


(2-  41f) 


Here  Q is  the  rotational  partition  function  and  £"  ir  hcB  /kT,  the 
R v 

dimensionless  ratio  of  rotational  quantum  energy  to  kinetic 
temperature. 

The  first  term  in  the  asymptotic  expansions  for  Voigt  profile 

b 

band  functions  is  easily  obtained,  as  is  illustrated  for  the  L (a, 7") 


function.  At  large  optical  paths,  L(a,T)  for  a single  line  is  from 
Equation  (2-21b)  approximately  27T*^(aT)*  Since  the  optical 
thickness  at  the  center  of  a line  in  a band  is  simply  gj'f’.  then 


b,  oTT-1/4  1/2  V 1/2 

L (a,T)  2/7  a (Tg_) 

J J 

V 1/2 

or,  denoting  /jgJ  by  G^2 
J 


and 


L (a,TW 


7T 


1/4 


1/2 

(aT)  G 


1/2 


(2-42a) 


(a,TK 


(a/T) 


1/2 


rr 


1/4  1/2 


(2-42b) 


, i/2 

M2  (a,D^a  1/2 


27f1/4r3/2 


(2-42c) 


L0(a,T)^ 


4(aT) 


1/2 


srr 


1/4  1/2 


(2  -42d) 
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(a,T) 


G 

1/2 


(2-42e) 


Nn  <a-T)' 


a 1 /2 

3rr1/4T3/2 


1/2 


(2-42f) 


In  order  to  simplify  programming,  a single  set  of  functions 
is  used  for  all  parallel  bands  of  linear  molecules  and  a single  set 
is  used  for  perpendicular  bands.  These  functions  vary  with  the 
rotational  temperature  and  in  the  program,  in  order  to  partially 
account  for  temperature  variation,  the  parameter  T is  not  used, 
but  rather  a derived  parameter  T , which  is  the  optical  thickness 
at  the  center  of  the  strongest  line  in  the  band.  This  is  either  an  R 
branch  line  for  a parallel  band  or  a Q branch  line  for  a perpendicular 
band. 

The  effects  of  temperature  variation  will  be  illustrated  by  the 
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results  of  calculations  for  the  L (o,T*)  function  for  a perpendicular 

band  of  a linear  molecule.  The  discussion  assumes  the  Doppler 

profile  because  it  is  possible  to  obtain  a series  expansion  for  the 

high  temperature  limit  in  the  variable  T . Use  of  double  ore- 

o r 

cision  computation  permits  accurate  calculation  into  the  region 

b 

where  the  logarithmic  approximatio n to  L (o,T)  becomes  good. 

The  band  used  as  an  example  is  a "perpendicular"  transition 

with  upper  state  a 7 T level  and  lower  state  a level.  It  is 

assumed  that  the  molecule  is  not  symmetric  and  there  are  no 

nuclear  spin  effects.  The  relative  strength  of  each  line  is  the 

product  of  the  exponential  part  of  the  Boltzmann  factor  of  the 

rotational  energy  level  of  the  lower  state  and  the  appropriated 

Hoenl- London  transit. on  probability  or  factor.  For  the  Jth  rotational 

level,  the  Boltzmann  factor  is  exp[-  (f  J(J+1)1  where  € = hcB  /kT  and 

v 

h is  Planck's  constant,  B^  is  the  rotational  energy  constant  for  the 
lower  vibrational  state,  and  k is  Boltzmann's  constant.  The 
Hoenl- London  factors  for  the  P,  Q,  and  R branches  are 


Q _ 2J  + 1 
J 2 


CR  _ J + 2 
"j  ~ 2 


(2-43a) 


(2-43b) 


(2-43c) 


and  have  been  normalized  so  that  the  sume  for  a given  value  of  J 


The  rotational  partition  function  is  the  sum  over  all  rotational 
levels  of  the  product  of  the  rotational  degeneracy  factor  2J  + 1 
and  the  Boltzmann  factor  for  each  level.  At  high  temperatures, 
an  asymptotic  approximation  is  available  for  the  partition 
function. 

OO 


The  J value  of  the  strongest  line  in  the  Q branch  may  be 
found  by  setting 

J (2J  + 1)  expf-  Cj(j  + 1)]=  0 
and  solving  for  J.  The  result  is 


Adopting  ~ 1 /(•'  , the  fraction  of  band  strength  in  the 
strongest  line  is  V?/2  exp(-  1/2  + C/4)  and  the  optical  thickness 
at  the  center  of  this  (fictitious)  line  is 

T0  = W C/2  exp  (-  1/2  + C/4) 

using  k = T \fe  / 2 exp  (1/2  - C/4),  the  optical  thickness  at  the 
o 

center  of  any  line  in  any  branch  is 

= k (J  - 1)  exp  [-  C J(J  +1)] 

J o 


(2-44a) 


(2 -44b) 


(2-45) 


(2-46) 


(2-47) 


(2-47a) 
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(2-47b) 


7~?  = k (2J  + 1)  exp  [-  Cj(J  + 1)] 

J 


T^=  k(J  + 2)  exp  [-  CJ(J  + 1)] 

J 


(2-47c) 


These  expressions  are  then  substituted  in  Equations  (2-36) 
and  2 - 1 5 ) to  obtain 


L (0,7';  €~ 


r°°r-- 

1 kT<P,Q,R)e-X  \ 

r L 

1-e"  J 

1 

J J;P,  0,  R 
•'-00 

i 

dx 


(2-48) 


“ X 

where  f(o,  x)  = e is  used  for  the  Doppler  profile.  The  sum  over 
J may  be  replaced  by  an  integral  over  z = J(J  +1)  and  the  summands 
1 - exp[-  kf  (z,  x)]  expanded  in  series  and  integrated  over  x and  z. 
After  simplification, 


co  , — .n-1  (n-  1 ) / 2 

V ' ( ” C _ ) e 


L (o,r  r L n;Y^ 

n=  1 » 


(2n“1+  l)(2n)"(n+1)/2 


r 


n+1 


2 

(2-49) 


Use  of  double  precision  arithmetic  on  the  CDC  6600  permits 


evaluation  of  Equation  (2-49)  for  values  of  7?  up  to  about  50. 


Define  the  functior. 


S(T)  = Lb  (o,T) 


(2-50) 


This  function  is  plotted  in  Figure  (H-4)  for  values  ofTbetween 
zero  and  ten,  using  the  series  ct  Equation  (2-49)  for  £ = 0 and 
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numerical  integration  for  € = 0.  002,  0.  01  and  0.  05.  On  the 
scale  of  this  drawing,  the  curves  of  £ = 0 and  £ = 0.  002  are 
indistinguishable,  the  values  at  t=  10  being  0.  3501  and  0.  3483 
respectively.  For  larger  values  of  “X,  it  is  convenient  to  plot 
TS(T)/lnT  versus  (ln*tT)  This  is  shown  in  Figure  (II - 5 ) for 
values  of'l'between  10  and  50  for  £ = 0,  and  between  10  and  10^^ 
for  £ = 0.  002,  0.  01  and  0.  05.  At  T?  = 50,  the  value  for  £ = 0 is 
about  one  percent  higher  than  that  for  £ = 0.  002.  The  value  of 
€ = 0,05  is  about  eight  percent  lower  than  that  for  £ = 0.  002. 

The  finite  values  of  6 used  here  correspond  to  rotational  tempera- 
tures of  about  280K,  56K,  and  1.  1 K for  carbon  dioxide  and  1400  K, 
280K,  and  56K,  for  carbon  monoxide. 

The  asymptotic  behavior  of  the  other  radiative  transfer 
functions  is  similar  in  the  high  temperature  limit,  taking  into 
account  the  differences  in  the  leading  term  of  the  expansion  for 
large  optical  depth.  In  the  upper  atmosphere,  the  region  in  which 
radiative  transfer  is  the  most  significant  factor  affecting  vibrational 
populations  is  at  such  an  altitude  that  the  Doppler  shape  is  the  main 
contributor  to  the  effects  of  radiative  transfer. 

The  next  step  is  to  fit  these  radiative  transfer  functions  into 
the  context  of  obtaining  a solution  to  Equation  (2-1)  and  the  extensions 
of  Equation  (2-1)  when  the  effects  of  collisional  excitation  and  de- 
excitation are  added. 

A definite  lower  limit  to  the  source  function  at  the  far  boundary 
may  be  obtained  by  assuming  the  source  function  to  be  linear  in***. 

In  the  gray  case,  this  leads  to 


J(T)  = 


1 - 2 C3(T o) 

-T  ) 

2(1  + T - e °; 
o 


(2-51) 
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where  T ^ again  refers  to  the  total  optical  path  length  through 
a finite  plane  layer.  This  expression  is  never  in  error  by  more 
than  13.4%,  when  compared  with  King's  (1956)  exact  solution. 
Making  the  same  approximation  for  line  or  band  radiative  transport 
functions  gives  the  result 


J(T  ) 

O 


2L  (T  ) 

3 o 

r + L(r  ) 

o o 


(2-52) 


This  becomes  increasingly  inaccurate  as  the  optical  path  through 

the  layer  increases,  because  of  the  long  range  effects  of  line  shape. 

Ivanov  (1972)  has  obtained  an  symptotic  solution  to  two  r^.uiative 

transfer  problems  involving  line  radiation  with  arbitrary  line  shape. 

The  first  is  that  of  a plane  layer  of  finite  optical  thickness  Tf  in 

o 

which  one  boundary  layer  is  in  radiative  quilibrium  and  is  the  sole 
source  of  radiative  flux  to  the  medium.  The  corresponding  integral 
equation  is 


) = 
o 


1 1 -T  | ) dt 


+ 


A 

2 


n21(T) 


(2-53) 


1 


The  parameter  ^ allows  for  loss  of  radiation  by  "true"  absorption 
or  collisional  de  - excitation.  The  asymptotic  solution  to  this 
equation  when  *£*  is  very  large  is 
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(2-54) 


<£a(T,To) 


2 N2i(Y) / 1 -A  + A Ni 2('^o)  ' 

[1  - A f A N12(T)]3/2  1 1-A  + * N12(ro"T)j 


The  second  problem  is  that  of  a plane  layer  with  a constant  excitation 
mechanism  and  no  external  source  of  radiation.  The  integral 
equation  is 


S(T) 


(T-  t)  S(t)dt  + 1 


(2-55) 


and  the  asymptotic  solution  is 

sa(T)  = [i  -A  + Ani2(T)]'1/2  [i  -A  + A n12(To-T)]"1/2 

(2-56) 

This  last  result  may  then  be  used  to  provide  an  asymptotic  solution 

to  the  earthshine  problem  with  no  collisions,  in  which  case  A is  equal 

to  one.  A lower  limit  to  the  excitation  from  below  is  -7-  N (*J*  ). 

2 12  o 

From  symmetry, 


$,(To/2'  T„!  = 7 <2-57a> 

$(7.ro)  = i -$rro  -r,ro)  (2-s7b> 


A lower  (in  the  asymptotic  case  of  large  X ) limit  to  the  source 

o 

function  for  ^greater  than  X / 2 is  then 

o 


tVTo» 


tN12(T>  N12«To  -T'l 


,r 


r / 2 

o 
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(2-58) 


This  gives  low  values  in  the  interior,  and  a bridging  formula 

which  gives  a good  initial  solution  which  can  be  refined  by  iteration 
is 


^a(7'ro)  = itN12(T)/N12(T0  -T,]1/2’  T'ro n < 

This  can  be  used  with  Equation  (2-57b)  to  cover  the  entire  range 

of  T' from  zero  to*?  . 

o 

In  practice,  rather  than  formally  integrate  a modified  Equation 

(2-1)  including  aTdependent  A to  account  for  collisional  de-excitation, 

and  an  additional  source  term  to  account  for  collisional  excitation, 

the  computer  program  sets  an  initial  vibrational  population  at  each 

altitude  and  computes  the  total  radiative  excitation  at  each  level 

which  results  from  the  sum  of  radiation  from  the  earth  and  lower 

atmosphere,  the  sun  and  other  levels  in  the  atmosphere. 

The  band  functions  are  used  to  determine  the  effects  of 

absorption  between  a molecule  at  a given  altitude  and  the  various 

sources  of  radiation  affecting  a given  band  in  the  following  way. 

The  vertical  optical  paths  to  the  lower  boundary  and  to  each  level 

in  the  atmosphere  are  computed  and,  for  daytime  conditions,  the 

optical  path  to  the  assumed  top  of  the  atmosphere  in  the  direction 

of  the  sun  using  the  secant  of  the  solar  zenith  angle  as  a multiple. 

The  effective  radiation  from  the  lower  boundary  is  determined  by 

multiplying  the  surface  radiancy  by  the  appropriate  Nb,  function. 

The  effective  radiation  from  the  sun  is  determined  by  use  of  the 

proper  function.  The  effective  radiation  from  the  other  levels 

of  the  atmosphere  may  be  obtained  by  integrating  the  product  of 

emission  at  each  level  and  the  band  function  Nb 

21’ 


This  is  no 


i 

I 


i 

I 

i 

[ 

I 


longer  done  directly,  however,  as  it  was  discovered  that  the 
earlier  integration  subroutine  gave  incorrect  results.  Instead, 
the  assumption  is  now  made  that  the  relative  vibrational  uopu- 
lation  between  two  points  varies  linearly  with  optical  thickness. 
Then  the  integral  relations  between  the  band  functions  are  used 
in  the  form 


(A  + Bt)  t^dt 


' J<Tl!  NU  (a'Tl>  - J(T2>  Nl\ 


(a,T2) 


J(T2) 


J(Tl) 


T2-Tl 


L3(A,rL) 


b _ 

L3(a,T1) 


(2-60) 


Vibrational  populations  are  determined  by  balancing  colli  - 
sional  and  radiative  excitation  and  de-excitation  rates.  Total 
band  intensities  at  an  observation  point  inside  or  outside  the 
atmosphere  are  computed  by  determining  the  optical  path  from 
each  level  to  the  observation  point  and  integrating  the  product  of 
local  emission  rate  and  the  band  function  (a,TT).  The  spectral 
distribution  of  radiation  is  determined  by  finding  a mean  effective 
temperature  and  line  of  sight  column  density  for  the  radiating 
molecules,  computing  intensity  for  each  line  in  the  band  with 
the  line  function  L(a,'T),  summing  the  line  emission  rates  and 
scaling  each  line  emission  rate  so  that  the  total  band  intensity  is 
obtained. 


3 

\ 
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The  above  analysis  has  been  presented  for  the  general  Voigt 
profile.  Actual  use  of  the  radiative  transfer  functions  in  the  pro- 
grams has  been  restricted  to  the  Doppler  profile.  This  restriction 
is  adequate  for  the  water  vapor  and  ozone  bands  which  in  the  con- 
centrations used  in  the  program  have  a maximum  vertical  optical 
thickness  at  line  center  less  than  3.  For  limb  viewing  of  carbon 
dioxide  bands  the  Voigt  profile  should  be  used.  Efforts  to  find  a 
general  way  of  setting  up  the  radiative  transfer  functions  to  use  a 
Voigt  profile  that  varies  with  altitude  have  not  been  successful. 
Exploratory  calculations  show  that  the  computed  radiance  values 
for  strongest  bands  at  a tangent  height  of  60  kilometers  may  be  in 
error  by  20  — 50  %.  In  view  of  the  uncertainties  in  chemical 
abundances,  this  error  is  not  felt  to  be  serious. 

The  same  radiative  transfer  functions  have  been  used  for 
the  non-linear  molecules  as  were  used  for  the  linear  molecules. 

The  match  between  radiative  transfer  functions  has  been  made  by 
determining  a single  "effective"  rotational  constant  for  the  non-linear 
tri-atomic  molecules.  This  was  done  by  computing  actual  radiative 
transfer  functions  for  a limited  range  of  optical  thicknesses  adequate 
to  cover  that  needed  in  the  model.  The  computed  functions  were 
compared  with  those  for  the  linear  molecules  and  the  effective 
rotational  constant  needed  to  produce  an  approximate  match  at 
optical  thickness  at  line  center  of  about  2 was  used. 


SECTION  III 


ATMOSPHERIC  COMPnsmnM 


There  have  been  relatively  few  measurements  of  the  variation 
with  altitude  of  infrared  emitting  species  at  altitudes  above  40  km,  the 
l.mit  of  direct  sampling  by  balloons.  While  ozone  measurements  by 
different  workers  agree  reasonably  well  with  each  other  and  with  theo- 
retical computations  at  altitudes  up  to  100  km,  and  nitric  oxide  profiles 
appear  reasonably  well  known  between  about  80  and  150  km  at  middle 
lattitudes,  corresponding  data  are  lacking  for  other  species  which  are 
important.  There  is  in  some  cases  indirect  data  available  which  when 
combined  with  computational  models  can  yield  limits  to  concentrations 
Thus.  Anderson  ,1,71  a,  has  been  able  to  put  an  upper  limit  on  total 
hydroxyl  radical  abundance  above  60  km  from  measurement  of  resonance 
fluorescence  and  with  this  result  computes  possible  water  vapor  pro- 
files between  50  and  100  km  using  a photochemical  model  that  includes 
vertical  transport  by  eddy  diffusion.  In  a later  paper,  Anderson  (1971  b, 

gives  a measured  altitude  profile  for  hydroxyl  with  error  limits  t 120 
per  cent. 

At  present,  then,  estimates  of  chemical  abundances  of  minor 
species  required  to  predict  infrared  radiation  require  the  use  of 
photochemical  computations.  Even  with  the  inclusion  of  best  estimates 
of  the  effects  of  such  transport  processes  as  eddy  mixing,  these 
computations  must  be  used  with  caution,  even  when  they  appear  to 
agree  with  the  limited  concentration  determinations  that  are  available, 
Wofsy  and  McElroy  (1,73)  give  a discussion  of  methane  concentrations 
in  the  stratosphere  which  illustrates  the  interplay  of  experiment  and 


modelling  studies,  showing  how  estimates  of  eddy  mixing  coefficients 
may  be  obtained. 

A knowledge  of  the  rates  of  the  appropriate  chemical  reactions 
is  necessary  in  order  to  estimate  the  abundance  of  trace  constituents 
at  high  altitudes  in  the  earth's  atmosphere.  A set  of  reaction  rates, 
with  estimates  of  probable  error  in  the  rates  was  listed  in  the  final 
report  on  the  previous  research(Degges,  1972).  This  list  was  neces- 
sarily incomplete  and  several  of  the  rates  quoted  were  only  estimates. 
Several  reactions  have  been  measured  only  at  room  temperature  and 
others  only  at  flame  temperatures.  In  the  latter  case,  only  upper  or 
lower  bounds  may  be  known. 

Since  that  report  was  completed,  a series  of  evaluations  of 
chemical  reaction  rate  coefficients  has  been  published  by  workers 
at  the  National  Bureau  of  Standards  as  a part  of  the  Department  of 
Transportation  Climatic  Impact  Assessment  Program.  The  most 
recent  of  these  compilations,  edited  by  Garvin  and  Hampson  (1974) 
provides  a thorough  summary  and  evaluation  of  the  accuracy  of  rate 
coefficients  for  neutral  chemical  relations  and  photochemistry, 
vibrational  energy  transfer  processes,  high  temperature  air  reactions, 
and  ion-molecule  reactions.  These  reports  are  the  main  source  of 
data  used  in  updating  the  chemistry  program. 

Among  the  most  important  new  measurements  which  have  been 
made  are  those  of  Davis,  Herron  and  Huie  at  temperatures  between 
200  and  350  K for  the  reactions 

O + NO  ► NO  + O 

Ld  Li 

and 

O + Oz  + Ar ► 03  + Ar 
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(3-i) 

(3-2) 


Reaction  (3-1)  is  found  to  be  temperature  independent  with  a rate 


-12  3 

constant  9.  1 x 10  cm/molec-sec  in  contrast  to  the  exp  (-300/T) 


temperature  dependence  of  the  analysis  of  Baulch,  et  al  (1970)  which 
we  have  used  previously.  Reaction  (3-2)  is  found  to  have  a tempera- 
ture dependence  of  exp  (500/T)  rather  than  the  exp  (1060/T)  given  by 
Johnston  (1968)  and  is  slower  at  low  temperatures  than  previously 
believed.  This  implies  slower  ozone  formation  and  slightly  greater 
atomic  oxygen  concentrations  at  high  altitudes. 

Many  reaction  rates  remain  poorly  determined,  and  the  compilation 
of  Garvin  and  Hampson  points  these  out  clearly.  One  such  reaction  is 


OH  + HO. 


H2°  + °2 


(3-3) 


- 1 1 - 12 

where  measurements  of  the  rate  range  between  10  and  2x10 


3 -1-1 

cm  molecules  s . They  suggest  that  model  calculations  should 


include  the  extreme  ranges. 

Water  vapor  abundances  are  not  well  known,  even  in  the  stratos- 
phere. Balloon  measurements  are  limited  to  altitudes  below  40  km 
and  are  subject  to  considerable  uncertainty  as  to  the  amount  of  water 
vapor  carried  to  high  altitude  by  the  balloon.  Direct  measurements 
put  lower  and  upper  limits  of  the  volume  mixing  ratio  at  two  and  ten 
parts  per  million  in  the  upper  stratosphere.  A recent  optical  measure- 
ment from  a balloon  has  been  made  by  Goldman,  et  al  (1973)  who 
observed  thermal  emission  from  the  rotational  spectrum  near  25 
micrometers.  Near  and  below  30  km,  the  volume  mixing  ratio  is 
nearly  three  parts  per  million.  In  the  region  between  30  and  36  km, 
above  the  peak  altitude  reached  by  the  balloon,  the  mixing  ratio 
appears  to  be  two  or  three  times  as  large.  Scholz,  et  al  (1970) 
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obtained  comparable  results  in  a direct  sampling  from  a rocket  near 
50  km.  Hunten  and  Stroebel  (1974)  have  reviewed  measurements  and 
computations.  They  present  a model  which  assumes  a mixing  ratio 
of  four  parts  per  million  at  30  km.  This  includes  the  water  producer 
by  oxidation  of  methane  in  the  stratosphere.  Figure  III- 1 shows  the 
results  of  their  computation  (solid  line)  and  the  water  vapor  profile 
used  by  Degges  (1972),  (dashed  line).  Hunten  and  Strobel  (1974)  show 
that  their  results  are  relatively  insensitive  to  variations  in  the  eddy 
mixing  processes.  Except  for  very  large  changes  in  the  mixing  ratio 
at  the  30  km  level,  water  vapor  concentrations  at  all  altitudes  will 
vary  directly  with  the  low  altitude  mixing  ratio.  They  find  that  better 
agreement  can  be  obtained  with  Anderson's  (1971  a,  b)  measurements 
if  half  as  much  water  is  assumed  to  be  present  in  the  atmosphere. 

Figure  III- 2 shows  a typical  nitric  oxide  concentration  curve, 
based  on  mid-latitude  measurements. 

While  ozone  abundances  below  about  50  km  must  take  into  account 
lateral  transport  as  well  as  vertical  transport  by  eddy  mixing  and 
molecular  diffusion,  photochemical  models  are  adequate  above  that 
altitude,  and  agree  with  observations  obtained  by  measuring  absorption 
of  solar  and  stellar  ultraviolet  radiation  and  those  deduced  from 
measurements  of  the  molecular  oxygen  1.  27 yUm  infrared  atmospheric 
band.  The  exact  profile  depends  on  assumed  eddy  diffusion  coefficients, 
however,  principally  through  atomic  hydrogen  and  oxygen  concentration, 
abundances  used  here  are  shown  in  Figure  III-3. 

Carbon  dioxide  has  a lifetime  of  about  10  days  at  the  top  of  the 
atmosphere,  where  it  is  dissociated  by  solar  radiation.  It  is  well 
shielded  by  molecular  oxygen  below  about  120  km,  however,  and 
should  remain  well  mixed  below  100  km.  The  only  good  measurement 
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above  100  km  remains  that  of  Offerman  and  von  Zahn  (1971),  who 
find  that  between  120  and  140  km  the  mixing  ratio  is  not  more  than 
50  per  cent  less  than  the  tropospheric  value.  Above  140  km,  molec- 
ular diffusion  is  rapid  enough  compared  to  photodissociation  to  main- 
tain diffusive  equilibrium.  Figure  III- 4 is  a typical  estimate  of  CO, 
concentrations. 

Recent  work  by  Liu  and  Donahue  (1974)  has  resulted  in  a model  for 
H20,  CH^,  H2  and  odd  hydrogen  that  properly  relates  the  measured 
mixing  ratios  in  the  stratosphere  to  escape  of  atomic  hydrogen  from 
the  atmosphere  by  the  mechanisms  of  the  Jeans  thermal  flux,  charge 
exchange  and  the  polar  wind.  Their  water  vapor  concentrations 
derived  for  a 1000  K exospheric  temperature  differ  from  that  of 
Figure  III- 1 by  less  than  a factor  of  two  for  altitudes  below  100  km. 

At  higher  altitudes,  their  values  become  somewhat  larger  than  those 
of  the  drawing. 

The  implication  is  that,  since  the  chemistry  program  used  here 
does  not  include  the  proper  boundary  conditions  for  hydrogen  escape  at 
high  altitudes,  long  time  calculations  for  hydrogen  containing  species 
may  be  in  error.  Future  modifications  to  the  chemistry  program  should 
include  a provision  for  adding  the  hydrogen  escape  flux  as  part  of  the 
high  altitude  boundary  conditions. 

The  chemistry  program  is  described  and  listed  in  Appendix  A,  in 
almost  the  same  form  as  used  by  Degges  (1972).  In  addition  to  modi- 
ficat  on  of  the  boundary  conditions,  it  may  be  possible  to  improve  the 
efficiency  of  the  integration  method.  The  remainder  of  this  section 
describes  the  results  of  an  investigation  of  a modification  of  the  basic 
algorithm. 
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The  numerical  method  used  to  integrate  the  chemical  rate 
equations  in  the  chemistry  program  assumes  that  the  rate  equation 
for  a single  molecular  species  may  be  written 


x - e^<  x 


where  q,  J] , and  cX  are  time  dependent  functions  of  the  rate  co- 
efficients and  the  concentrations  of  the  other  molecular  species  in  the 
system.  The  analytical  solution  of  equation  (3-4)  for  q,  Jj , and  <=>< 


constant, 


x(t)  = 


[(/}  +Cf)exp(-  Cft)  + (Cf  + 2q[l  - exp(-  Cft)] 

”(2o <x+J]  )[1  - exp(-Cft)]  + Cf [1  + exp(-Cft)] 


where 


c f = (y3 2 + 4o<q) ] 


is  applied  within  the  framework  of  the  second  order  improved  tangent 
Runge  - Kutta  method  as  described  by  Degges  (1972).  This  algorithm 
is  stable  for  any  value  of  the  time  t and  avoids  the  necessity  of 
assumptions  such  as  that  some  reactants  with  extremely  short  time 
constants  are  in  chemical  equilibrium. 

When  the  product  (ft  is  large  and  any  of  the  quantities  q,  JJ , 
and  cX  are  changing  with  time,  however,  true  time  dependence  is 
lost.  For  large  Cft,  equation  (3-5)  becomes 

x(t)  = = J Ijj  ’ Cft  large  (3_6 
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and  the  value  computed  for  the  time  t actuallv  corresponds  to  time 
t/2  when  Equation  (3-5)  is  embedded  in  any  second  order  Runge  - 
Kutta  method.  In  effect,  the  basic  step  is  only  first  order  for  molecular 
species  for  which  O^t  is  large.  The  extrapolation  algorithm  does  con- 
verge to  a value  for  the  concentrations  of  such  components,  but  the 
process  may  require  more  steps  than  if  the  basic  step  were  truly 
second  order,  as  it  is  for  small  values  of  Cft,  or  whenq  , JJ , and 
are  constant. 

Mod'fications  of  the  algorithm  has  been  attempted  as  follows. 
Consider  the  differential  equation 


which  is  Equation  (3-4)  with  the  x term  absent.  The  formal  solution 


for  q and  JJ  constant  is 

x(t)  = x exp(-yj  t)  + (q /J}  )[1  - exp(- fl  t)]  (3-8) 

If  q and  p>  are  functions  of  time, 

q = qo  + qlfc  + q2fc2  + (3-9a) 

ft*  ftQ  + At  + /^2t2+ (3’9b) 

The  improved  tangent  method  applied  over  an  interval  h results  in: 

x( h/2)  = x + (h/2)(q  - B x ) (3-10) 

O 0^00  ' 
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(3-11) 


dx 

dt 


= qQ  + (h/2)qi  + (h/2)  q2  + 


- [/?o  + (h/2)/^  + (h/2 )Zfl2  + ...][xq  + (h/2)(qo  +/?oXq;] 


and 


x(h)  = x + 
o 


h — 
hdt 


h/2 


* *o  + h%  + (h2/2)1i  - hflo\  - ,h2/2)  lXo 
Z 2 3 4 

- (h  / 2)  o xq  + terms  in  h , h , ...  (3-12) 

which  is  correct  through  in  terms  in  h1*. 

When  Equation  (3-8)  is  embedded  in  the  improved  tangent  method, 
the  result  is: 


x(h/2)  = xoex-'.  pQt)  + (qQ//?o)[l  - expf-^t)]  (3-13) 


dx 

dt 


= qo  + (h/2)qi  + (h/2)2q2  + 

h/2 


(3-14) 


- [ J3Q  + (h /Z)fi1  + (h/2 )ZJ3Z  +....]  x (h/2) 


and 


x(h)  = xq  exp[=  - (h/2) ] 

qQ  + (h/2)q^  + . . . . 


(3-15) 


+ <h/2 )/?!  + ..  1 ■ exp  hA,(h/2)A  - •••] 
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A series  expansion  of  Equation  (3-15)  is  identical  to  Equation  (3-12) 
through  terms  in  second  order.  The  higher  order  terms  are 
different,  but  the  method  is  still  second  order.  For  large  /J h, 
however,  Equation  (3- 1 5)  becomes 


x(h) 


qQ  + (h/2)q1  + (h/2)2q2  + . . . 

fiQ  + (h/2  )fi1  + (h/2  )ZJ3Z  + ... 


(3-16) 


and  not  the  correct  result 


x(h)~ 


q + hq  th  q + .... 
o 1 22 

ftQ  + + h 2 + * * * * 


(3-17) 


When  h is  large,  the  algorithm  may  be  altered  as  follows  to 
retain  second  order  accuracy  and  obtain  a better  appro  ximation  to 
the  endpoint  value.  Let: 


x = q /ft  + 0 


Then 


d0  _ _d 
dt  dt 


= 4>-M 


(3-18) 


and 


0 - 0oexP  (-JJ t)  + (0/0)[  1 - exp(-/?t)] 


(3-19) 
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further,  let 


0 = 0//?  +¥ 


-R-fiV 


(3-20) 


4/  *yoexp(-j3t)  + (R/yQ)[l  - exp(-y^t)] 


(3-21) 


Only  terms  in  q^,  q^,  Q>  and  J3^  enter  into  the  second  order  terms 
in  the  solution  of  Equation  (3-7)  given  in  Equation  (3-12).  If  q and  J3 
are  assumed  linear  in  the  interval  of  integration  h,  q = q(t  = 0), 

Po  = P{t  = 0)’  = 2[q(h/2)  - qo]/h,/?j  = 2[  (h/2)  -/]Q]/h,  and  if 

4>o  * xo  - %'Po’ 


x(h/2)  = q J fio  +0oexp(-^oh/2) 


(3-22) 


Then  if  i'o=0c 


/},q(h/2)  - qj  /?Q(h/ 2) 


and  the  term  R is 


/T(h/?> 

neglected,  '^(h)  = ^oexp[-h^(h/2)]  and  the  resulting  expression 
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...  qo  + hqi  A%-Aqi  . „ 

x(h)  ~~ 7p^J} TIT3  [1 ' exp('  h/?°  ‘ h /?l/2,] 


+ (x° ' zf  ’ exp('  h^°  ’ h2A/2) 


(3-; 


2 

when  expanded  in  series,  agrees  through  terms  in  h with  Equations 
(3-12). 

The  same  analysis  may  be  applied  to  Equation  13-4).  In  this  case 
the  appropriate  substitution  is; 


JLa_. 


x = CT+  b + 


0 


and 


^ = " dt  0+  b ‘ CT0 


= 0-CT0 


(3-; 


setting 


A =CT  + 4c <(p  , make  the  further  substitution 


0 = 1'  + 


A+cf 


to  obtain 


v/  - _ -d  ) 2 0 

dt  iA+/?j 


A y -cx^1 


= R - A 4/  -<*!' 


(3-2 
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Again  neglecting  the  R term  in  Equation  (3-25),  the  method  becomes 


x(h/2)  = -.lqo  . CTq  0Qexp(-  0>/2) 

\Do  - exP(“  C^h/2 )]0O  + 

The^o'0^'  A'  A'  V ql'  and  ^ are  evaluated  as  before  frc 
the  values  available  at  t = 0 and  t = h/2  and; 


(3-26) 


:c(h)  __  2q(h)  , 2<$(h/2)  , A(h/2)^oexp[-  hA(h/2)] 

CT(h)+^(h)  A(h/2)+CT(h/ 2)  + ^(h/2{l  -exp[ - h^  (h/2)  \jTT+  ^(h/2) 


When  A2  is  negative,  the  soluti 


(3-27a) 


on  must  be  modified  to 


15(h)  . A(h/2)tan{arctan[2o<(h/2)^+^(h/2)/A(h/2)]-l/2  A(h/2)  A -rfth /2^ 

'~cT(h)+3ih)  + ‘ — — 

^ 2cx.(h/2) 


(3-27b) 


The  best  value  of/?h  or  Oh  to  divide  the  region  of  applicability, 
of  the  two  methods  was  studied  by  sample  computations,  but  results 
v/ere  not  appreciably  better  than  with  the  basic  algorithm. 
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SECTION  IV 

RESULTS  FROM  COMPUTATIONAL  MODEL 

A.  Introduction 

The  revised  high  altitude  spectral  radiance  programs  can 
produce  spectral  radiance  profiles  for  wavelengths  between  2.  5 and 
25  micrometer >j  or  between  400  and  4000  wave  numbers.  Results  are 
presented  here  based  on  estimated  abundance  profiles  of  the  four 
species  CO^,  H^O,  ^3  anc^  The  temperature  and  pressure  pro- 

files used  were  taken  from  the  U.  S.  Standard  Atmosphere,  1962. 
Upwelling  terrestrial  and  lower  atmosphere  radiances  are  assumed  to 
be  those  corresponding  to  a black  body  source  at  the  60  km  temperature. 

Spectral  radiance  plots  are  shown  for  typical  tangent  heights 
for  day  and  night  conditions.  Features  of  the  spectral  radiance  curves 
are  explained  with  the  aid  of  altitude  profiles  for  total  band  radiances 
of  individual  vibrational- rotational  bands. 

B.  Limb-Viewing  Spectral  Radiances 

Representative  limb  viewing  radiance  profiles  are  shown  in 
Figures  IV-1  through  IV-6.  These  are  based  on  the  concentration 
profiles  of  Section  III,  except  that  the  experimental  curve  of  Roble 
and  Hays  (1974)  has  been  used  for  night  time  ozone  and  a constant 
mixing  ratio  has  been  used  for  CO^.  The  temperatures  and  number 
densities  used  are  those  of  the  U.  S.  Standard  Atmosphere,  1962. 

The  results  shown  are  quite  different  from  those  of 
Degges  (1972)  for  many  of  the  carbon  dioxide  bands,  where  the 
earlier  version  of  the  program  predicted  appreciably  lower 
radiance  values  at  hign  altitudes.  The  night  time  ozone  radiance 
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FIGURE  IV- 1 
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FIGURE  IV- 3 
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FIGURE  IV -4 
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FIGURE  IV-6 


values  are  also  higher  in  spite  of  a reduction  in  total  column  count 
2 

{molecules/cm  ) above  60  km.  The  differences  in  the  CO  values 

2 

are  primarily  because  the  early  program  did  not  correctly  inte- 
grate the  radiation  field  produced  by  molecules  radiating  above 
60  km. 

Comparison  of  day  and  night  profiles  for  a given  tangent 
height  shows  that  the  low  energy  (small  wave  number  and  long  wave 
length)  band  radiances  do  not  vary  appreciably  during  the  day.  The 
high  altitude  high  energy  bands,  however,  do  vary  in  intensity  by 
as  much  as  two  orders  of  magnitude  as  the  result  of  resonant  and 
fluorescent  emission  following  absorption  of  solar  radiation.  This 
is  shown  clearly  in  the  CO^  4.  3 micrometer  bands  and  in  the  water 
vapor  complex  at  2.  7 micrometers.  A faint  indication  of  this  effect 
is  visible  in  the  5.  4 micrometer  nitric  oxide  fundamental  where  the 
day  time  curve  is  broadened  because  of  the  fluorescent  contribution 
from  absorption  of  sunlight  at  2.  7 micrometers  and  in  the  ultra- 
violet. 

C.  Individual  Band  Radiances 

For  some  purposes,  it  may  be  more  convenient  to  have  available 
the  total  radiance  of  an  individual  band  rather  than  the  spectrum  of  the 
region  where  the  band  is  located.  This  may  aid  in  making  allowances 
for  species  abundances  different  from  those  assumed  in  a computation 
or  in  assigning  an  experimentally  observed  radiance  value  to  one 
species  rather  than  another. 

Figures  IV- 7 and  IV- 8 show  the  day  and  night  total  band 
radiances  for  the  ten  water  vapor  bands  included  in  the  spectra  of 
figures  IV-1  through  IV-6,  The  6,3  micrometer  V band  is  the 
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strongest,  day  or  night.  During  the  day,  the  V ^ band  at  2.  66  micro- 
meters strongly  absorbs  solar  radiation  and  becomes  almost  as 
bright  at  a tangent  height  near  90  km.  At  night  this  band  is  nearly 
four  orders  of  magnitude  fainter. 

Figures  IV- 9 and  IV- 10  show  night  and  day  limb  radiance 
values  for  seven  bands  of  the  CO^  16  micrometer  band  complex,  the 
4.  3 micrometer  V band  and  one  of  the  two  fluorescent  bands 
originating  from  the  state,  the  10.42  micrometer  00011  - 10001 
Band.  The  15  micrometer  band  radiances  vary  little  between  day  and 

night.  During  the  day,  however,  the  radiance  of  the  4.  3 micrometer 

-9  -8  2 

band  rises  from  its  night  value  of  between  10  and  10  watts/cm  -ster 

- 8 - 7 2 

to  a value  of  between  10  and  2 x 10  watts/'cm  -ster.  This  band  is 

thick  in  the  limb  viewing  geometry,  and  most  of  the  increase  is  due  to 

increased  vibrational  temperatures  at  altitudes  above  90  km.  The 

10.41  micrometer  band,  hich  has  much  lower  optical  thickness, 

gives  a better  indication  of  the  altitudes  at  which  the  solar  excitation 

- 8 2 

is  important.  The  60  km  value  is  about  4 x 10  watts/cm  -ster 
day  or  night.  Figure  IV- 11  shows  night  limb  radiance  values  for 
the  same  CO^,  bands  when  the  mesopause  temperature  is  raised  to 
230  K.  The  main  effect  is  to  double  or  triple  the  band  radiance 
values,  indicating  the  importance  of  the  temperature  Drofile  in 
computing  radiance  values. 

D.  Comparison  with  Observations 

A limited  amount  of  data  from  the  ICECAP  program  is  avail- 
able for  comparison  with  the  results  of  the  radiance  computer  pro- 
gram. Figure  IV-12  compares  peak  spectral  radiance  values  for 
the  15  micrometer  band  measured  by  Stair,  et  al  (1974)  with  those 
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Zenith  peak  spectral  radiance  at  15/J.m  measured  with 
LWIR  spectrometer  aboard  rocket  flown  at  Poker  Flat, 
Alaska  on  22  March  73  compared  to  theoretical  model. 
The  solid  curve  is  the  measured  profiles  (ascent  and  de- 
scent) while  the  dashed  curves  show  model  calculations 
of  Kumer  and  James  (KJ)  and  Corbin  and  Degges  (D). 
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FIGURE  IV- 12 


ZENITH  SPECTRAL  RADIANCE 


computed  by  the  radiance  program  (D  Model)  using  the  U.  S. 

Standard  Atmosphere,  1962  temperature  profile  witli  180  K meso- 
pause  temperature  and  that  computed  by  J.  B.  Kumer  (K.  J.  Model) 
using  the  Kumer  and  James  (1973)  model  radiance  program.  The 
differences  between  the  results  of  the  two  computer  programs  are 
due  principally  to  differing  temperature  profiles  and  assumed  CO 

2 

mixing  ratios.  The  degree  of  agreement  between  the  two  models 
and  the  experimental  data  give  an  indication  of  what  can  be  expected 
from  use  of  the  computer  programs  and  the  need  to  improve  both  the 
computational  algorithms  and  knowledge  of  atmospheric  composition. 

Figure  IV-13  compares  the  radiance  model  prediction  with 


the  data  of  Stair,  et  al  (1974),  this  time  for  the  ozone  V ^ band.  The 


amount  of  ozone  used  for  the  calculation  is  optically  thin  in  the  up- 
ward direction,  and  doubling  the  amount  would  double  the  computed 
radiance.  The  effective  width  of  the  observed  ozone  band  is  about 
one-third  micrometer,  and  the  plotted  experimental  data  agrees  with 
the  measurement  to  within  a factor  of  two.  The  shape  of  the  observed 
radiance  curve  supports  the  accuracy  of  the  ozone  measurement  of 
Hays  and  Roble  and  their  interpretation  that  their  observed  ultra- 
violet absorption  is  due  to  ozone  and  not  a mesospheric  aerosol 
layer.  Figure  IV-13,  as  well  as  Figure  IV-12,  shows  that  at  least 
some  results  of  the  radiance  program  are  in  agreement  with  measured 
radiances,  and  suggests  that  with  more  accurate  data  on  atmospheric 
temperature  and  composition,  observations  will  either  confirm  model 
predictions  or  show  how  modeling  algorithms  must  be  modified. 


E. 


Model  Evaluation 
The  most  important  general  areas  of  investigation  in  this 
research  effort  have  been  (a)  the  chemistry  and  spatial  variation 
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Zanith  peak  spectral  radiance  at  9.6^Am  measured 
with  LWIR  spectrometer  aboard  rocket  flown  at 
Poker  Flat,  Alaska  on  22  March  73  compared  to 
theoretical  model.  The  solid  curve  is  the  measured 
profile  while  the  dashed  curve  shows  the  D Model 
calculations  in  watts  cm-^  ster-^. 


FIGURE  IV- 1 3 


ZENITH  SPECTRAL  RADIANCE 


of  atmospheric  species  important  to  infrared  radiation  and  (b)  the 
improvement  of  models  for  the  transport  of  infrared  radiation. 

Section  III  dealt  with  atmospheric  composition  and  modifications 
to  the  chemistry  program  written  during  a previous  research  effort 
(Degges,  1972).  This  program  remains  limited  in  its  usefulness 
to  short  time  extrapolation  from  a given  set  of  atmospheric 
composition  profiles.  It  can  be  used  to  give  limited  information 
on  such  problems  as  auroral  chemistry,  but  requires  knowledge 
of  the  limitations  imposed  by  the  limited  chemistry  set  and  un- 
certainties in  global  transport  of  such  important  species  as  atomic 
oxygen. 

More  important  advances  have  been  made  in  the  radiative 
transport  program.  A serious  programming  error  has  been 
corrected,  and  an  improved  formulation  of  the  transport  problem 
has  been  coded.  The  program  produces  results  which  agree  with  a 
limited  amount  of  observational  data.  The  agreement  with  the  carbon 
dioxide  radiance  model  of  James  and  Kumer  (1973)  is  believed  to  be 
satisfactory.  There  remain  areas  in  which  results  were  not  as  good 
as  those  anticipated  when  the  work  was  initiated. 

The  program  for  computing  band  radiances  requires  an  excessive 
amount  of  computer  time  for  bands  that  are  optically  thick  in  the  upward 
direction  above  the  altitude  at  which  collisions  are  important  in 
determining  vibrational  populations.  The  spectral  radiance  program 
is  much  slower  now  that  a relatively  complete  set  of  spectral  lines 
is  used  for  the  water  vapor  and  ozone  bands.  The  first  of  these 
problems  can  be  resolved  by  using  a matrix  inversion  to  obtain  a 
first  approximation  to  the  vibrational  populations  of  the  optically 
thick  bands,  as  is  done  in  the  program  developed  by  James  and 
Kumer.  The  second  problem  is  a result  of  trying  for  completeness. 


and  a smaller  set  of  lines  may  be  selected  and  used  for  some 
purposes,  using  more  complete  results  as  a guide  to  what  may  be 
safely  omitted. 

The  program  could  be  improved  in  some  other  ways  by  adding 
features  from  programs  such  as  that  of  Whiting,  et  al  (1969),  which 
was  written  to  compute  line-by-line  spectra  of  diatomic  molecules. 

Use  of  this  program  with  either  of  the  Voight  profile  subroutines 
used  in  the  current  research  effort  would  permit  improved  calculation 
of  limb  radiances,  provided  composition  data  are  sufficiently  accurate 
to  justify  the  additional  programming  effort  and  computer  time. 

Finally,  we  use  Figures  IV-14  and  IV-15  to  illustrate  one 
remaining  problem.  These  figures  give  night  and  day  spectral 
radiance  values  for  a limb  view  with  a tangent  height  of  f JO  km. 
Comparison  with  Figures  IV-1  through  IV-6  shows  that  the  band 
widths  are  appreciably  wider  than  for  the  lower  atmosphere  where 
the  temperature  is  lower.  The  program  assumes  that  the  effective 
radiation  received  from  the  lower  atmosphere  and  from  the  lower 
parts  of  the  upper  atmosphere  is  such  as  to  produce  a constant 
vibrational  temperature  above  110  or  120  km.  Above  these  altitudes, 
collisions  play  no  part  and  excitation  is  entirely  by  radiative  transfer. 
The  widening  of  the  bands  at  very  high  altitudes  will  permit  absorption 
of  radiation  in  the  wings  of  the  band,  and  this  effect  is  not  properly 
accounted  for  in  the  program.  Additionally,  at  an  altitude  of  500  km, 
it  is  probable  that  the  rotational  temperature  is  not  in  equilibrium  with 
the  kinetic  temperature.  Further  work  is  required  to  investigate  the 
actual  rotational  temperature  when  collisions  are  not  sufficiently  rapid 
to  maintain  rotational  temperature  in  equilibrium  with  kinetic  temperatu 
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FIGURE  IV- 15 
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SECTION  V 


HIGH  ALTITUDE  RADIANCE  PROGRAMS 

The  high  altitude  (>60  km)  radiance  program  set  includes  four 
programs,  LOWRAD,  BCKGND,  SPCTRA  and  TRYCRT.  Each 
program  may  be  used  separately  or  in  appropriate  combination  with 
one  or  more  of  the  others  within  the  CDC  6600  operating  system. 

LOWRAD  is  used  only  to  provide  an  approximation  to  the  spectral 
radiance  in  the  upward  direction  from  the  lower  atmosphere.  Input 
to  LOWRAD  is  from  two  separate  files,  TAPE  1 which  holds  the 
premanent  data  and  data  for  five  model  atmospheres  and  TAPE  5 = 
INPUT  which  directs  the  choice  of  one  of  the  five  atmospheres  or 
another  atmosphere,  in  which  case  data  for  this  atmosphere  must 
also  be  on  TAPE  5.  Printed  output  is  written  to  TAPE  6 = OUTPUT 
and  output  in  a form  suitable  for  use  by  BCKGND  is  written  to 
TAPE  2,  which  may  be  a temporary  or  permanent  disc  file  or  a 
magnetic  tape. 

BCKGND  computes  vibrational  populations  and  total  band 
radiances  for  an  earth's  limb  viewing  geometry  at  5 km  invervals 
between  60  and  120  and  at  25  km  intervals  between  120  and  500  km. 
Control  information  is  read  from  TAPE  5 = INPUT  and  permanent 
data  is  read  from  TAPE  1.  Low  altitude  spectral  radiances  ere 
read  from  TAPE  2.  Atmospheric  composition  at  2 km  intervals 
between  60  and  150  km  and  at  5 km  intervals  between  150  and  700  km 
may  be  read  either  from  a sample  transferred  from  the  UPDATE  set 
to  TAPE  1 or  from  TAPE  5.  Output  for  printing  is  written  to  TAPE  6 
OUTPUT  and  to  TAPE  3 for  storage  or  for  use  as  input  to  program 
SPCTRA. 

Preceding  page  blank 
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SPCTRA  computes  spectral  radiances  corresponding  to  the 
total  band  radiances  output  by  BCKGND.  Control  information  is 
read  from  TAPE!  5 = INPUT.  Permanent  data  is  read  from  TAPE  1 
and  the  results  of  the  BCKGND  computation  are  read  from  TAPE  3. 
Output  data  is  sent  to  TAPE  6 = OUTPUT  for  printing  and  to  TAPE  4 
for  use  by  the  plotting  program.  Optional  output  of  unfiltered  spectral 
radiances  may  be  placed  on  TAPE  7 and  saved  for  later  use. 

TRYPEN  and  TRYCRT  plot  spectral  radiance  curves  using 
values  computed  by  program  SPCTRA.  TRYCRT  produces  microfilm 
output;  TRYPEN  uses  a pen  plotter. 

The  programs  are  listed  and  described  below. 
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Program  LOWRAD 


LOWRAD  is  used  to  provide  an  approximation  to  the  spectral 
radiance  in  the  upward  direction  from  the  lower  atmosphere.  It  is 
based  on  the  computer  program,  LOWTRAN  1,  written  by  J.  E.  A. 
Selby  of  AFCRL.  LOWTRAN  1 provides  numerical  values  based 
on  the  atmospheric  transmission  model  of  McClatchey,  et  al  (1970). 

LOWTRAN  1 is  capable  of  computing  transmittances  with  20 
wavenumber  resolution  at  wavenumbers  between  350  and  40,000  cm-* 
(0.25  to  28.  5 micrometers).  Computations  are  restricted  to  points 
with  numerical  value  a multiple  of  five  wavenumbers.  LOWTRAN  1 
uses  atmospheric  parameters  (temperature,  pressure,  and 
composition)  at  one  kilometer  intervEtlS  between  0 and  25  km,  at 
five  kilometer  intervals  to  50  km,  and  at  70  and  100  km.  Variables 
required  at  intermediate  altitudes  are  obtained  by  interpolation. 

The  modifications  to  LOWTRAN  1 to  make  LOWRAD  c.re  simple 
and  straightforward.  Beginning  at  350  cm  , transmission  between 
60  km  and  each  lower  level  of  the  atmosphere  is  computed  for  5 
nadir  angles.  The  Planck  function  at  that  wavenumber  is  computed 
for  ground  level  temperature  and  the  mean  temperature  between 
altitude  levels  used  by  the  program.  The  ground  level  value  of  the 
Planck  function  is  multiplied  by  the  transmittance  between  ground 
level  and  60  krr  altitude.  This  and  the  product  of  the  Planck  function 
and  transmittance  differences  for  the  intermediate  levels  are  summed 
to  produce  spectral  radiance  values  at  each  of  the  given  angles.  The 
version  of  LOWRAD  listed  here  uses  the  five  angles  arcsin 
V (21  - 1 ) / 1 0 with  I taking  values  one  to  five.  Equal  weights  are 
assigned  to  the  radiances  computed  with  each  value  of  I,  and  the 
arithmetic  mean  is  taken  to  be  the  value  of  the  upwelling  spectral 
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radiance.  An  effective  radiation  temperature  is  computed  from  the 


final  spectral  radiance  value  for  use  (if  desired)  with  program 
BCKGND.  This  procedure  is  repeated  at  5 cm  1 intervals  until 


•1 


4000  cm  (2.  5 micrometers)  for  a total  of  731  output  radiance 
values.  These  output  values  are  written  both  to  the  line  printer 
and  to  a disc  file  for  later  use. 

Figures  V-l  through  V-5  show  the  resulting  values  of 
radiation  temperature  for  the  five  model  atmospheres  provided 
by  McClatchey,  e t al  (1970).  The  effects  of  the  strong  absorption 


bands  of  CO^,  H^O  and  are  easily  seen  in  the  lower  temperatures 


at  the  most  intense  parts  of  the  bands,  where  a detector  looking  down 
would  in  effect  see  the  cold  lower  stratosphere. 

Table  V-l  lists  radiation  temperature  values  for  the  five  model 
atmospheres,  (1)  tropical,  (2)  midlatitude  summer,  (3)  midlatitude 
winter,  (4)  sub-artic  summer  and  (5)  sub-artic  winter. 

The  radiance  values  and  effective  radiation  temperatures 
computed  by  LOWRAD  are  limited  in  usefulness  by  the  low  re- 
solution available.  A molecule  in  the  upper  atmosphere  acts  as 
its  own  spectrometer  with  a resolution  comparable  to  its  line  width 
which  is  of  the  order  of  between  one  hundredth  and  one  tenth  of  a wave 
number.  For  most  bands,  a better  value  of  effective  radiation 
temperature  is  that  of  the  kinetic  temperature  at  the  altitude  taken 
to  be  the  base  of  the  high  altitude  region.  This  is  60  km  for  the 
high  altitude  radiance  program  BCKGND. 

Input  to  LOWRAD  is  from  two  separate  files,  TAPE  5,  which 
hr  Ids  permanent  transmission  data  and  data  for  five  model  atmosphere 
and  TAPE  4 which  directs  the  choice  of  one  of  the  five  atmospheres  or 
another  atmosphere,  in  which  case  data  for  this  atmosphere  must  also 
be  on  TAPE  5.  Printed  output  is  written  to  TAPE  6 = OUTPUT  and 
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235. 75 
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7’7.hl 

254. 44 

122= 

351.°2 

.M3.?: 

76  5.4'. 
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244. 27 

243.41 

234.13 

2"  1 . 17 

226.90 

?34" 

243 . 51 

7*7 • 53 

233. 2 f 

2cn . 30 

226.21 

23-45 

243 . 33 

2*7. 35 

233.11 

260 . 08 

226.04 

2350 

244.21 

243. 35 

234.06 

251 . 11 

226.86 

2355 

2 44  . 85 

249. C 2 

234. 70 

2C1 . 78 

227.38 

2 3 60 

244.97 

249. 14 

234.91 

251 . 39 

227.49 

2 3 b5 

242 . 00 

246. C 5 

231 . 31 

?'*  * . ,3  2 

225.09 

2370 

231 . 24 

235,63 

223 .Qt 

239.11 

219.34 

2 3 75 

223.05 

232.53 

22'*.  19 

275 .89 

21,3.34 

2 3 8 0 

232 . 30 

234.35 

224,10 

775 . 18 

220.05 

23  3^ 

243 . 38 

247.05 

274.06 

247.70 

227.7? 

23  j 0 

262.73 

2b9.73 

244.15 

254 . 54 

235.74 

2395 

374.62 

270.69 

253.56 

76'4.53 

242.38 

2 4 j 0 

7*4. 50 

279. 85 

260.6? 

273. 18 

248.95 

2425 

239 . 28 

23  4.2  3 

264. ?7 

277.41 

251.81 

2410 

2°1 . 64 

256.4J 

265.06 

279.52 

253.21 

2415 

291 .04 

286.75 

265.29 

279  .73 

" 253.36 

2420 

392.21 

286,99 

266.48 

7 3 0 . 02 

253.49 

2425 

292 . 32 

287.09 

255.56 

28].  13 

253 . 54 

2h30 

29?  . 4-* 

237.21 

265 • 6C 

283 . 24 

253.59 

2 4 35 

?Q2 . 62 

237 , 35 

265.77 

7n0 . 39 

253.66 

2440 

292  . 74 

237. 46 

255 .85 

280 .50 

253. 70 

2445 

292 . 98 

287.  6 7 

267.57 

280 . 71 

~ 25'3.137 

2 453 

2°3 . 3G 

237.95 

257.25 

?30 . 99 

253.97 

2*55 

293 . 74 

2,8  8 . 37 

267,60 

281 . 39 

254.19 

2460 

294 . 16 

28  8.7-5 

267. 9T 

7 81  . 77 

254.40 

2 465 

294. 43 

238.99 

263.1 n 

287 .Ci 

254.53 

2470 

7 94 . 77 

289. 30 

263.35 

28? . 71 

254.69 

2475 

795 . io 

2o9.6T~" 

2ET75‘fT~' 

'7874  "51 

~2?4Tff4 

2460 

790.44 

289.91 

263.84 

2 8 ? , 91 

235.03 

2435 

295 . 88 

290.30 

257.17 

283.31 

255.27 

2490 

29  6. 22 

290.61 

269.42 

283 . 61 

25  5. 40 

2 495 

296.54 

290.50 

269.66 

283 .90 

235.57“ 
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M0P5L 

9 

7C  .0 

296.61 

231.2  + 

26  js 

? <7 . ’4 

2)1.53 

25ic 

?-“*7 . 33 

291.67 

26  16 

’97.46 

’31.75 

25  20 

?°7 . 50 

’31.73 

25  29 

•’97.39 

2 a 1 . r.  l 

?r  3? 

’97.15 

221. 4'’ 

2539 

’96.  «7 

231. 21 

’64? 

’96.57 

2°  : . p ♦ 

2 •“ -9 

’r'6. 75 

2 Jv . 7* 

?:pD 

? >6.24 

? -j  C .63 

2 5-35 

”6.27 

2 3 G . 7 j 

25  62 

296.31 

29  G • 73 

25  o'5 

296.16 

2 c«  i • 6 -> 

?r -r. 

’95 . 52 

’4  0.u? 

25  7C 

’ '>5 .66 

? * C . 2 3 

253' 

’95.67 

23G . 27 

?J>  =5 

’96.56 

29C.t? 

icjri 

296.55 

291.*) 

? = 'i5 

’97 . 22 

2J1.5* 

25  2 2 

297.32 

291, 3 ? 

25:5 

297.53 

292. C7 

261? 

56 

’=>2.  r ) 

25  15 

”7.57 

2a’. 12 

2522 

’97.3*; 

2?  1.6-5 

26  25 

297.25 

291. 83 

263? 

’97.15 

231.73 

25  35 

79^.92 

291.61 

26  4? 

’96 .73 

231. 4. 

26  45 

296.5? 

291.. 7) 

265? 

296.19 

291 . 15 

2655 

?9f .13 

291.21 

26  » 3 

’96.10 

2 3 ' . 9 ) 

26  63 

’96.15 

’°1.C? 

267'’ 

’96.39 

291.21 

28’3 

’’6.44 

?’l. 24 

26)0 

’96. 58 

291.75 

26  65 

’96 . 7? 

291.51 

269? 

296.61 

231.54 

2695 

296.91 

291.6? 

27  C ? 

296.96 

291.63 

27  j5 

296.61 

291.0) 

271? 

296.04 

29C.95 

AT:M> 

o4r^r 

3 

4 

5 

269.94 

284.24 

255.75 

270.1 7 

’"♦.64 

2-35.91 

270 . 75 

’“4.62 

255.95 

’73 . 36 

284.74 

256.02 

’77.45 

’14.76 

256.03 

271. 3? 

?«4 .69 

255.97 

773.16 

284.49 

255.34 

’63.96 

’84. ’4 

255.69 

’69.75 

287.99 

233.8? 

269.6? 

’87 ,60 

255.79 

’6  A .65 

’•3.72 

255.36 

’63.6  r 

283.77 

255.40 

263.67 

’•7 , o’ 

255.46 

269.6’ 

’•7.75 

25  5. 4> 

26  ) . 4 6 

2°7 .51 

255.30 

’6  9 . 3 C 

’•7.71 

’55.20 

269.36 

’87 . 35 

755.26 

263.71 

’37.73 

255.52 

’7j.it 

’“4.19 

255.87 

2 7 3 .54 

234.63 

256.13 

273.76 

234 .91 

256.36 

271.99 

’•5.12 

256.52 

271. *-7 

283.16 

256.59 

”1.  ]«■ 

’"5.19 

256.61 

271.04 

’ 0 5 . 35 

256.63 

273.9° 

284 . 97 

’56.53 

2’ 3.9 6 

?• 4 . al 

256.56 

27). 87 

234.74 

256.53 

271.6  0 

?»4 .5° 

256.50 

273.73 

’84.47 

256.43 

2”  .6? 

’34 .24 

256.43 

’7  3 • 6t 

’•4.20 

236.47 

270.59 

2 ? 4 . 1 6 

256.43 

’73.6’ 

2 • 4 .21 

256.44 

279.70 

2 3 4 . 73 

256.47 

273.79 

’“4.42 

256.43 

’71.77 

’“4.52 

256. 50 

273.85 

234.66 

256.54 

271 . 8 7 

234.69 

256.55 

27-3.91 

’84 . 7o 

’56.57 

273.07 

2“ 4. 75 

256.57 

’73.35 

’•4.55 

25  6 . 52 

270.61 

234,16 

756.44 
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2 

7 

4 

c 

7715 

295.35 

290.4+ 

270 . 37 

?p3 . 69 

256.34 

2721 

294 . 8C 

29  0 . Cl 

270 . IP 

2*3 . 32 

256.26 

2729 

295 . 06 

290. 22 

270 . 26 

“28  7.4  !> 

276.29 

2 7 T0 

7^5.95 

29  0.33 

270 . 5 6 

284  . P 8 

256 • 41 

2735 

296 .66 

231.41 

2-0 . 79 

2*4 .56 

256.49 

2740 

297.14 

291.76 

273.94 

284 . 88 

256.54 

2745 

297. 20 

291.89 

270.95 

2*4 .92 

256.54 

2750 

2.97 . 01 

291.63 

27  } . 8 8 

284 . 77 

256. 50 

2755 

296. 80 

291.47 

270.80 

“7*4 .63“ 

“236  .“47 

2760 

296.56 

291.32 

270.72 

2*4.48 

256.47 

2765 

296.32 

291.12 

'’7  3,61 

284.29 

256.37 

2773 

295. 05 

29 j.9j 

270.50 

294.09 

256.71 

2775 

7«5 . 79 

29C.73 

270 .38 

2*7.90 

236.25 

27*2 

295.55 

290.51 

270.27 

2C>7  . 72 

25  6.13 

2785 

295.37 

29  C . 35 

77  0 n^- 

783.59 

2 55  .“15 

2790 

995.24 

2 9 0 . 2 7 

270 . 16 

7*3.51 

256.14 

2 7 ‘ J 5 

295 . 29 

293.31 

775.19 

2*2.35 

256.15 

2fc  C 0 

295 . 29 

290.34 

570.23 

287.59 

256.21 

28  05 

295.27 

290. 34 

270 .25 

287.59 

256.24 

2810 

295 . 17 

290.27 

270  .?  3 

2*3.53 

256.25 

2915 

295.01 

■ 79  on* 

~ 370 .17  ' 

“7  67;  4 2“ 

7255.77 

2820 

294.63 

269.8’ 

270.03 

2*7 . 16 

256.16 

2«25 

294.67 

239.6 ) 

270.04 

2*7. 18 

256.16 

2830 

294. 75 

289.95 

270.07 

2°  3.24 

256.13 

28  35 

294.97 

290. 13 

270.16 

2*5 ,4C 

256.22 

2640 

295.37 

29 j.42 

27  0 . 3 5 

28  7 , 68 

256.29 

2 845 

295 .73- 

““295".  79' 

“ 270 .46- 

— 7ir3".7T5 — 256.76  - 

2P50 

796. C2 

29  3.91 

270.57 

7*4.11 

256. 38 

2855 

296.20 

291.01 

270 . 55 

2*4.21 

256.37 

2860 

295 . 17 

290.93 

27  T . 51 

7*4. 17 

256. 37 

2665 

296 . 12 

290.94 

27‘J  .45 

2*4 . 12 

256.29 

2872 

295.95 

290.3] 

27}  . 34 

2*7.98 

256.21 

”“28  77" 

' 395V63" 

297.  49  ' 

‘ 77T.T7 

75T.”69 

335T7£T“ 

2880 

295. 20 

29  G . 1 5 

269. 86 

2*3.36 

255. 90 

2635 

294 . 85 

*26  9. 35 

769.66 

2*  T. 08 

255.74 

2890 

29u.24 

289. 3? 

269.30 

2*2.59 

255.51 

2395 

293.48 

288.71 

269.9? 

2 8 ? , 0 1 

255.27 

2900 

293. 29 

236.51 

268 .73 

2*1  .81 

255,10. 

2937 

~2  9 ty  rg™  Tsirrrr 

*767T5TT~ 

? r i . g,  7 

’25  4795  "* 

2910 

2«2. 79 

236.03 

268.26 

2*1  . 32 

254.72 

2915 

292.45 

237.71 

263  . t?  6 

2*lV55 

254.59 

2923 

291 . 85 

287. 23 

267.75 

2*0.59 

254.40 

2925 

290 . 81 

786. 37 

267.25 

?’,9. 80 

254,09 

— . . 
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Wfl  \/5 

000*1 

NJ.IBF7 

1 

0 

u 

7 

4 

5 

29  39 

239. 24 

235.1’ 

266.57 

?73 .64 

753.  73 

2935 

288.37 

23^.40 

?6d . 2 1 

?~7 .96 

257.55 

29  93 

237.59 

2 3 3 . 6 ♦ 

26  5.  ? 6 

l7’ ,43 

253.33 

29U9 

237 . Cb 

23  3.3? 

765.55 

375 . 95 

733.16 

2990 

235.59 

2 3 2 . 9 i 

26 1 . 3 17 

3-3 . 58 

233.05 

29  53 

236.47 

34?  . 84 

265.3  2 

776. 50 

253.05 

29  oO 

'3P5 . 12 

232.5’ 

265.2? 

376 ,?4 

253.03 

2969 

7«3 . 74 

282. ?7 

265.1  r> 

271,96 

25  3.0  1 

2970 

234.95 

231.61 

264 .37 

?75 .40 

35?.9.3 

2979 

233.93 

288. 8+ 

264.5? 

’*"•*  . 64 

252.87 

2983 

233.40 

250.-2 

254  . .3  6 

2" 4. ?6 

252.34 

2 9 5 

2 5 3.  78 

280.75 

264,6? 

2^4.58 

253.02 

2990 

233.94 

280.9] 

264.74 

3 7'4.  71 

253.11 

2999 

234.  m 

2 31. C 2 

264. 7C 

3-4.81 

253.11 

It  01 

2P? ,?4 

279.31 

263 . 3 T 

377.15 

231.6? 

30  y3 

279 . 33 

27  5 . 4o 

260.39 

2’ 3 . 36 

25  0.05 

3 1 13 

276.33 

274.11 

259.0  7 

263 . 12 

348.65 

30  15 

279.11 

27  0 • 63 

255.39 

264.36 

246.66 

3020 

272 . 65 

27  0.37 

256.5  7 

364.53 

246.97 

30  25 

278  . 14 

272. cb 

253.99 

?66 . 99 

2h9.11 

30  33 

77 -.61 

275.20 

261.21 

?69 . 35 

25  0 . 94 

30  35 

2 79 . 50 

27  7 . C2 

262.48 

77] .01 

751.93 

3090 

281.25 

37  6.6  1 

263. 5 p 

27?. 53 

252.56 

3 0h5 

231 .25 

278.57 

263.57 

3*’?. 50 

232.4? 

30  53 

?«1 . 41 

278.61 

263.48 

?’?  .59 

252.40 

30  55 

231.30 

278.64 

263.4  4 

?7^>.56 

252.44 

30  60 

°51 . 16 

"278.53 

263.39 

7-?. 43 

232.43 

30  b5 

73. .7b 

273 . 11 

26?.:  T 

7-3.15 

252.23 

3 0 70 

2.31  .ug 

27  6. 8? 

26  3 , 4T 

?73 ,7? 

252.15 

30  75 

2*2.4*) 

279.6? 

267.8? 

?‘,3.30 

252.51 

3 C 0 

23? .40 

27  9. 5 3 

263 . 7^ 

273.40 

252.40 

3085 

281 . 92 

279.11 

263.44 

377.01 

252.24 

32  90 

737 , u4 

27  7 7 85 

262.74 

731 . 32 

251.40 

3 0 95 

2 7‘.  .91 

276. 51 

262.0? 

371.60 

251,54 

31C* 

276.62 

274.51 

260 . 9? 

263.70 

253.99 

3105 

272 . 91 

271.24 

251.01 

265.54 

250.00 

3110 

271 . 34 

369.79 

253 .15 

264.16 

249.55 

3115 

271. Q7 

26  9.64 

253.0? 

2*7.91 

249.43 

T123 

272 . 67 

271.  C-* 

?59  • 07  “ 

' 265',  75 

25  0;T3~ 

31  ’5 

374.63 

272.84 

263.15 

267.07 

250.79 

3130 

277. 70 

275.60 

261.6  7 

263.71 

251. 72 

3135 

231 . 52 

278.^1 

263.77 

272 . 85 

252.72 

3140 

2 =*4.47 

281.43 

265.20 

?’’5 , 23 

253.51 
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ATMOS0 

NUMBER 

1 

? 

3 

4 

5 

3145 

2 P.6 .52 

233.17 

265.21 

275 . 96 

254.07 

315C 

287. h4 

293.93 

266.77 

277.66 

254.43 

3155 

287 . 78 

234. 31 

267. Co 

2T7 . 98 

75 4756 

3160 

?87. 00 

283.72 

256.87 

27’  . 43 

254. 68 

3166 

2°  5 . 35 

262.38 

265.24 

276.14 

254.45 

3170 

284,46 

2 0 1 • 6 3 

263.44 

273.53 

254.37 

3175 

283. 70 

281. 04 

265.61 

??4 . 93 

254.24 

3183 

783.63 

230.99 

265.61 

274 . 36 

254. 2o 

3185- 

283 .69 

231.05" 

767.6  7 

?7V. 95 

254. Tr 

3190 

7P3 . 74 

261. 11 

765.72 

275  . QO 

25*+. 35 

3195 

2 63.91 

281.25 

263.81 

273 .15 

254.40 

32  00 

?«7 . 37 

230.79 

765.56 

274.71 

254.33 

3205 

281 .90 

279.54 

264.83 

273.53 

253.98 

3210 

280 . 51 

278. 34 

264. 2? 

272 . 36 

253. 35 

‘ 3215 

279 .93 

276.91 

263.47" 

■2°T.TI3~ 

253.2V 

3220 

278  . G7 

276.18 

262.97 

27] . 32 

253.01 

3225 

277.23 

775.47 

262.5? 

269.59 

252.78 

3230 

276. 90 

275. I1 

262. 3 r- 

269. 30 

252.68 

3235 

776.90 

275.12 

262 .35 

259  . ?1 

252.69 

32  hG 

777.43 

275  . c3 

262.63 

269 .77 

252.64 

3245 

277.77 

779  797 

218  2 ."8  7" 

"2TT.07 

252.93 

3250 

278.63 

276. 63 

263. 2 6 

?7] . 31 

253.17 

3255 

778.64 

276.69 

263.2  7 

77? . 81 

257.17 

3265 

278 . 48 

276.55 

263. IP 

779 .65 

253.13 

3265 

27? .15 

276.25 

263.01 

2^0 .40 

253.04 

32  70 

277.43 

27  5. 65 

262. 66 

2^9.82 

252.35 

3275 

277.15" 

■■779773' 

2" 577*+*“ 

252775" 

3230 

277.67 

27  5.  e.2 

262.76 

269.98 

252. 91 

32  85 

277 .85 

275.99 

262.86 

2T3 .14 

25  2.96 

32  5 0 

279.20 

277.13 

263.55 

271.29 

253.31 

3295 

780 . 16 

278. C4 

264.03 

. 09 

253.55 

3300 

280 . 79 

273.53 

264.33 

272. 60 

253.69 

T3R5 

281 . 23 

' 77-8  ;97‘ 

"264;  57  ' 

~27-r-7\T- 

' 753T7S- 

3 310 

7 8 1 . 22 

278.93 

264.52 

272 ,Q6 

25  3. 77 

3 315 

281 . 07 

278.81 

264. 44 

272.83 

253.  73 

3320 

281. C9 

278. 88 

264 .45 

272 . 85 

253. 74 

3725 

280 .96 

278.73 

764 . 4 3 

277. 75 

253. T2 

3330 

291.14 

278.69 

264. 5C 

272 .90 

253.7 7 

73  35 

■■281  .15" 

273TF9- 

""76V.  49  - 

2^7 .90 

' 25T.T7 ' 

3340 

231 .13 

273. 87 

264.47 

272 . 88 

253.75 

3345 

281.09 

27  8. 8 3 

264.44 r- 

277.-65 

25T.T2 

1350 

280.73 

278.51 

264.23 

272.53 

253.60 

3 3 55 

280 .52 

27  8 . 37 

264.10 

272  .-35 

253752 
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2 

7 

4 

5 

5 3 C- 

7 7 ->.71 

277.61 

963 , 6r 

971 .67 

253.70 

5 3 -j  3 

??'■>  . 75 

277,65 

267.72 

271  . 71 

253.73 

7 3 70 

7 6 C . 78 

278. 56 

264.26 

7?? . 58 

257.62 

5 *75 

7 81,63 

2 7 9 • 3 3 

264. 7f 

77  7 . ? 9 

253.36 

7 3 3 0 

7"  2 . 27 

279.66 

265 . 3 9 

277 . 05 

254.03 

> 3 (5 

.46 

28  0 . 0 5 

265.14 

274 . 01 

254.11 

* * 

9S? . 21 

279.8’ 

265. 07 

273 . 80 

254.06 

5 315 

7 81.93 

27 «. 57 

26 *.90 

?73 . 56 

253.91 

3 40  0 

9ai . 32 

279.13 

264. 6j 

277 . n& 

253.86 

74  05 

7 '■  1 . 33 

779. 06 

76 * . 6 1 

?73 .07 

053. 85 

34io 

9r  * . 87 

278.66 

264.39 

?7? . 69 

257.74 

3415 

?o: .08 

c 7 7 , 9 7 

264.00 

979  . u 3 

253.54 

34  ?n 

97(=,27 

977. ? 6 

267.54 

271 . 36 

253.7* 

3425 

?y<  ,?< 

776.39 

263.07 

2'"'  . 51 

253.0“ 

34  3D 

?76 . 93 

77  5.  16 

262.36 

269 . 74 

252.70 

74  35 

776.41 

274.61 

262.08 

9f  8 . 89 

252.55 

3 *•  -♦  0 

7 7 6 . 2 5 

274 . 5* 

261.99 

268 . 74 

252.50 

34*6 

2 7 6 . n7 

97  4.  7 5 

761.8  9 

9 6 3 . r'  9 

252.44 

3450 

775.71 

274.  C3 

261.66 

?r'?  . 27 

252.33 
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Program  SPCTRA 


Purpose: 


Program  Flow: 


SPCTRA  computes  the  spectral  radiances  in  the  5 to 
25//m  interval  corresponding  to  the  band  radiances 
computed  by  BCKGND.  The  input  data  to  SPCTRA  are 
the  molecular  band  line  constants,  molecular  vibrational 
level  and  band  constants  used  by,  and  vibrational  population 
and  band  radiance  data  generated  by  BCKGND.  In  itr 
present  form,  SPCTRA  computes  spectral  radiances 
between  2.  5 and  25  JJ.m  for  limb  views  and  endo-atmos- 
pheric  zenith  angles  of  0 and  90  degrees  at  5 km  intervals 
between  60  and  200  km. 

When  program  execution  begins,  values  are  assigned  to 
the  integer  variables  NOUT,  NLN  and  NALT  and  the  sub- 
routine ROTATE,  which  reads  in  data  for  the  molecular 
lines,  is  called.  The  outer  DO  range  ending  with  statement 
500  is  entered.  The  DO  parameters  are  fixed  at  present  to 
produce  two  sets  of  spectra,  using  night  and  day  input  from 
program  BCKGND.  The  program  may  be  made  more 
flexible  by  reading  in  DO  parameters  prior  to  execution  of 
the  DO  statement.  The  variable  BLAM  is  bp  t to  5,0  and 
the  outer  DO  range  ending  with  statement  1 J5  is  entered. 

The  effect  of  the  first  two  statements  is  to  set  up  a table 
of  wavenumbers  at  2 cm"1  intervals  between  2.5  and  25/U.m 
in  real  array  AVLAM.  These  are  the  finest  wavenumbers 
at  which  spectral  radiances  may  be  printed  out.  The  inner 
DO  range  ending  with  statement  105  sets  the  values  of  out- 
put spectral  radiances  equal  to  zero. 


The  DO  range  ending  with  statement  110  sets  temperature 
values  in  array  TEMP  to  those  corresponding  to  the  U.S. 
Standard  Atmosphere,  1962.  Temperature  values  may 
also  be  read  in  prior  to  execution  of  the  first  DO  statement. 
Following  final  execution  of  this  DO  range,  temperatures 
are  printed  out. 

Next,  the  number  of  molecular  species  is  read  in 
and  stored  in  integer  variable  NSPTOT.  The  outer  DO 
range  ending  with  statement  400  is  entered  and  executed 
NSPTOT  times.  Within  this  DO  range,  molecular 
vibrational  level  and  band  data  and  "ibrational  population 
data  are  read  in  and  mean  temperatures  and  total  number 
of  molecules  per  unit  area  in  the  vertical  and  horizontal 
directions  are  computed  for  each  vibrational  state.  Total 
band  radiances  corresponding  to  limb  viewing  and  vertical 
and  horizontal  endo-atmospheric  viewing  are  read  in. 
Relative  line  radiances  are  computed  and  normalized  to 
the  total  band  radiances.  The  spectra  corresponding  to 
the  normalized  line  radiances  are  computed  and  stored  in 
arrays  SPECA,  SPECB  and  SPECC.  Details  of  these 
computations  are  given  below. 

The  code  number  of  the  species  ISP  is  read  in. 

The  number  of  vibrational  levels  NX.EVEL,  the  number 
of  bands  NBAND  and  the  molecular  weight  RM  are  read 
in.  In  the  DO  range  ending  with  statement  210,  the  level 
description  codes  LSG  and  LCC,  the  energy  and  the  mean 
rotational  constant  for  each  vibrational  level  are  read  in. 

In  the  DO  range  ending  with  statement  220  the  band 
description  codes  LBC  and  LBU  and  the  band  strength 
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for  each  band  are  read  in.  The  number  of  altitudes  NALT 
for  which  vibrational  populations  are  computed  for  the 
species  and  the  number  of  altitudes  NR  LEV  for  which  band 
radiances  are  computed  by  BCKGND  are  read  in.  Vibra- 
tional level  populations  are  read  in  and  stored  in  array  CON. 

The  mean  temperatures  and  total  number  of  mole- 
cules per  square  centimeter  in  the  vertical  and  horizontal 
directions  at  5 km  intervals  between  60  and  200  km  are 
determined  in  the  DO  range  ending  with  statement  129. 

The  product  of  temperature  and  vibrational  level  popu- 
lation is  computed  in  the  DO  range  ending  with  statement 
112.  The  integrals  of  total  number  of  molecules  and  the 
product  of  temperature  and  number  density  in  the  vertical 
direction  are  computed  in  the  DO  range  ending  with  state- 
ment 120.  Simpson's  rule  integration  is  used.  The  total 
number  of  molecules  in  a given  vibrational  level  above  a 
given  altitude  is  stored  in  array  COLCNT.  The  corre- 
sponding mean  temperatures,  obtained  by  dividing  the 
integral  of  the  product  of  temperature  and  number  density 
by  the  integral  of  number  density,  are  stored  in  array 
AVTEMP.  After  exit  from  DO  range,  the  column  counts 
and  temperatures  are  printed  out. 

The  DO  range  ending  with  statement  128  performs 
a similar  computation  for  column  counts  and  mean 
temperatures  in  the  horizontal  direction.  The  integrations 
required  are  performed  in  function  QINT.  Upon  exit  from 
this  DO  range,  column  counts  stored  in  array  HORCNT 
and  mean  temperatures  stored  in  array  BVTEMP  are 
printed  out. 
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The  inner  DO  range  ending  with  statement  400  is 
now  entered  and  executed  for  each  band.  This  DO  range 
controls  the  computation  of  the  spectra. 

The  altitude  HT,  the  limb  band  radiance  RADC, 
and  endoatmospheric  radiances  RADA  (upward  direction) 
and  RADB  (horizontal  direction)  are  read  in  the  DO  range 
ending  with  statement  230.  Band  information  is  computed 
from  band  and  level  constants.  A test  is  made  on  the 
species  code  ISP  to  determine  whether  the  molecule  is 
linear  or  nonlinear. 

Following  the  test  of  ISP,  if  the  molecule  is  linear 
and  is  not  in  McClatchey' s tabulation,  line  positions  are 
computed  in  subroutine  ROTPOS.  The  DO  range  ending 
with  statement  250  is  entered  and  executed  for  each  of 
the  20  altitudes  for  .vhich  spectra  are  computed.  Line 
radiances  are  computed  in  subroutine  ROTSTR  for  vertical 
endo-atmospheric  viewing.  Successive  calls  to  subroutine 
FILTER  compute  the  spectral  contributions  from  lines  in  the 
P,  Q and  R branches.  The  sequence  of  calls  to  ROTSTR 
and  FILTER  is  rt  peated  for  endo-atmospheric  viewing  in 
the  horizontal  direction  and  for  limb  viewing.  Upon  exit 
from  this  DO  range,  control  is  transferred  to  statement 
280. 

If  the  test  on  ISP  determines  that  the  molecular 
line  data  is  available  control  is  transferred  to  statement 
260  and  the  DO  range  ending  with  statement  270  is  executed 
for  each  of  the  28  altitudes  for  which  spectra  are  computed. 
Spectra  are  generated  from  the  line  radiances  by  calls  to 
subroutine  WATSTR  and  FILTER,  for  oach  of  the  three 
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viewing  conditions. 

Following  the  final  execution  of  the  DO  ranges 
ending  with  statement  400,  the  spectra  are  printed  out  in 
the  order  endo-atmospheric  viewing  in  the  upward 
direction,  endo-atmospheric  viewing  in  the  horizontal 
direction  and  limb  viewing. 

Upon  final  execution  of  the  DO  range  ending  with 
statement  500,  the  program  halts  with  a STOP  statement. 


Mnemonic  Variables: 

AVLAM  Real  array  in  which  output  wavenumbers  are  stored. 

AVTEMP  Real  array  in  which  mean  temperatures  in  the  vertical 

direction  are  stored. 

BV  Real  array  in  which  rotational  constants  for  each 

vibrational  level  are  stored. 

BVTEMP  R«al  array  in  which  mean  temperatures  in  the  hori- 
zontal direction  are  stored. 

COLCNT  Real  array  in  which  total  number  of  molecules  per 

square  centimeter  above  60,  65,  . . . 200  km  in  each 
vibrational  level  are  stored. 

CON  Real  array  in  which  number  densities  for  each 

vibrational  level  are  stored. 

HHOLAM  Real  array  in  which  wavenumbers  of  band  lines  are 
stored. 

HHOLSE  Real  array  in  which  lower  state  energies  of  band  lines 

are  stored. 

HHOSTR  Real  array  in  which  strengths  of  lines  are  stored. 

HORCNT  Real  array  in  which  total  number  of  molecules  per 

square  centimeter  in  each  vibrational  level  in  the  hori- 
zontal direction  at  altitudes  60,  65,...  200  km  are  stored. 
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HR  AD 
P 

PRAD 

PS 

Q 

QS 

R 

RADA 

RADB 

RADC 

RCM 

ROTEA 

ROTLAM 

ROTLSE 

ROTST 

RS 


Real  array  in  which  line  radiances  are  stored. 

Real  array  in  which  wavelengths  of  P branch  of  linear 
molecule  are  returned  by  subroutine  ROTPOS. 

Real  array  in  which  products  of  number  density  and 
temperature  are  stored. 

Real  array  in  which  P branch  line  radiances  are  re- 
turned by  subroutine  ROTSTR. 

Real  array  in  which  Q branch  line  wavelengths  are  re- 
turned by  subroutine  ROTPOS. 

Real  array  in  which  Q branch  line  radiances  are  re- 
turned by  subroutine  ROTSTR. 

Real  array  in  which  R branch  line  wavelengths  are  re- 
turned by  subroutine  ROTPOS. 

Real  array  into  which  vertical  direction  and  endo- 
atmospheric  band  radiances  are  read. 

Real  array  into  which  horizontal  direction  endo-atmospheric 
band  radiances  are  read. 

Real  array  in  which  limb  view  band  radiances  are  read. 

Real  array  into  which  vibrational  level  energies  are  read. 
Real  array  used  by  subroutines  ROTATE  and  WATSTR, 
stores  factors  needed  to  compute  line  radiances  in 
water  vapor  pure  rotational  spectrum. 

Real  array  in  which  water  vapor  pure  rotational  line 
wavelengths  are  stored. 

Real  array  in  which  lower  state  energies  of  water  vapor 
pure  rotational  lines  are  stored. 

Real  array  in  which  absorption  strengths  of  water  vapor 
pure  rotational  lines  are  stored. 

Real  array  in  which  R branch  line  radiances  are  re- 
turned by  subroutine  ROTSTR. 
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SPECA 

SPECB 

SPECC 

STR 

TEMP 

BLAM 

CENTER 

HT 

HI,  H2,  H3, 
H4,  H5 

RBASE 

RB2 

RM 

RTOP 

SNB 

STH 

SVH 

TC 


Real  array  in  which  spectral  radiances  for  vertical 
viewing  endo-atmospheric  geometry  are  generated. 

Real  array  in  which  spectral  radiances  for  endo-atmos- 
pheric  horizontal  viewing  are  generated. 

Real  array  in  which  spectral  radiances  for  limb  viewing 
are  generated. 

Real  array  used  to  store  band  strengths. 

Real  array  used  to  store  atmospheric  temperatures. 

Real  variable  used  as  temporary  storage  during 
generation  of  values  in  array  AVLAM. 

Real  variable  used  to  store  wavenumber  of  band  center. 
Real  variable  into  which  altitudes  at  which  band  radiances 
are  computed  by  BCKGND  are  read. 

Real  variables  used  to  store  distance  increments  used 
in  computing  total  column  counts  and  mean  temperatures 
in  the  horizontal  direction. 

Real  variable  used  to  store  distance  of  base  altitude  from 
center  of  earth  in  computing  horizontal  column  counts. 
Real  variable  used  to  store  value  of  RBASE  2. 

Real  variable  used  to  store  value  of  molecular  weight. 

Real  variable  used  to  store  distance  from  center  of  earth 
of  given  altitude  level. 

Real  variable  used  to  store  current  value  of  band  strength. 
Real  variable  used  as  temporary  storage  is  computing 
horizontal  column  counts. 

Real  variable  used  as  temporary  storage  in  computing 
horizontal  column  counts. 

Real  variable  used  to  accumulate  integral  of  number 
density  in  computing  vertical  column  counts. 


124 


TCH 

TP 

TPH 


LBC 

LBU 

LCC 

LSC 

I 

ISP 

J 

K 


KB 

KK 

KL 

NA 

NALT 


Real  variable  used  to  accumulate  integral  of  number 
density  in  computing  horizontal  column  counts. 

Real  variable  used  to  accumulate  integral  of  product  of 
temperature  in  the  vertical  direction. 

Real  variable  used  to  accumulate  integral  of  product  of 
temperature  and  number  density  in  computing  mean 
temperature  in  the  horizontal  direction. 

Integer  array  storing  level  numbers  of  initial  and  final 
states  of  a band. 

Integer  array  storing  band  type  and  symmetry  information. 
Integer  array  storing  vibrational  level  information. 

Integer  array  storing  vibrational  level  information. 

Integer  variable  used  as  index  in  DO  loops. 

Integer  variable  into  which  species  identification  code 
is  read  (2  for  CO^,  3 for  H^O,  4 for  NO,  5 for  N^O  and 
6 or  7 for  ozone). 

Integer  variable  used  as  index  in  DO  loops. 

Integer  variable  (set  equal  to  28  in  current  program) 
used  as  final  parameter  in  DO  statement  controlling 
computation  of  column  counts  and  mean  temperatures. 

The  value  of  K is  the  number  of  altitudes  for  which  spectra 
are  computed. 

Integer  variable  used  to  store  band  type  code. 

Integor  variable  used  as  final  parameter  in  DO  range 
ending  with  statement  125,  computed  from  K and  I. 

Integer  variable  used  to  store  band  symmetry  code. 

Integer  variable  used  as  a .titude  level  index  in 
computing  horizontal  column  counts. 

Integer  variable  storing  total  number  of  altitude  levels. 
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NB 

NBAND 

NBC 

NCASES 

NF 

NI 

NL 

N.LEVEL 

NLN 

NRLEV 

NSP 

NSPTOT 

NZ 


Integer  variable  used  as  index  for  number  of  current 
band. 

Integer  variable  storing  total  number  of  bands  considered 
for  a given  species. 

Integer  variable  used  to  store  current  value  of  LBC. 
Integer  variable  used  as  index  for  spectrum  set. 

Integer  variable  used  as  index  for  lower  vibrational  state 
of  band. 

Integer  variable  used  as  index  for  upper  vibrational  state 
of  band. 

Integer  variable  used  as  index  for  vibrational  level 
number. 

Integer  variable  in  which  is  stored  total  number  of 
vibrational  levels  for  a molecular  species. 

Integer  variable  controlling  number  of  lines  treated  in 
each  branch  of  a band  of  a linear  molecule,  set  to  100  in 
current  program. 

Integer  variable  into  which  is  read  number  of  altitudes 
at  which  band  radiances  for  a given  species  are  computed 
by  BCKGND. 

Integ  r variable  used  as  index  for  DO  range  executed 
NSPTOT  times. 

Integer  variable  into  which  is  read  total  number  of 
species  to  be  included  in  a spectrum. 

Integer  variable  used  as  altitude  index  in  DO  loops. 
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CA.L  FILTER  ( SPEC , QS , Q , NL N ) 


Subroutine  FILTER 


Purpose  FILTER  computes  the  contribution  to  spectral  radiance 

from  each  line  in  a band  or  sub- band,  in  the  interval 
2.  5//m  to  25//m.  Provision  is  made  for  varying 
resolution  in  either  wavelength  or  wavenumber  space. 
Sufficient  storage  is  provided  to  compute  radiances  at 
2 wavenumber  intervals  with  triangular  slit  resolution 
of  5 wavenumbers. 

Method;  FILTER  searches  an  array  of  emission  line  strengths 

and  assigns  a weighted  strength  to  an  output  wavenumber 
or  wavelength. 

Program  Flow;  The  outer  DO  range  ending  with  statement  150  controls 
the  distribution  of  spectral  radiance  from  the  Ith  of  N 
given  lines,  and  computes  the  indexes  JS  and  JF  of  the 
first  and  last  output  points  in  array  SPEC  to  which  a given 
line  contributes.  The  inner  DO  range  ending  with  state- 
ment 120  computes  the  contribution  of  a given  line  to 
spectral  radiance  at  output  wavelengths  or  wavenumbers. 

Inputs;  RAD,  WVLTH,  N,  KOMPUT 

Outpuc;  SPEC 

Mnemonic  Variables; 

SPEC  Real  array  into  which  spectra  radiances  are  added 

at  regular  intervals  between  2.5  and  25//m  or  400 
and  4000  wavenumbers. 
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RAD 


WULTH 

N 

KOMPUT 

RADT 

RLAM 

RT 

I 

J 

JiT 

JS 


Real  array  containing  radiances  of  individual  lines. 
Real  array  containing  wavelengths  of  individual 
lines. 

Integer  variable  containing  number  of  lines  in  band 
or  sub-band  for  which  spectral  radiance  is  to  be 
computed. 

Integer  variable  which  determines  number  of  output 
spectral  points  and  whether  wavelength  or  wavenumber 
output  is  made. 

Real  variable  used  as  temporary  storage  for  radia.  e 
value  of  a single  line. 

Real  variable  used  temporary  storage  for  wavelength 
value  at  which  spectral  radiance  is  computed. 

Real  variable  used  as  temporary  storage  for  wave- 
length of  a single  line. 

Integer  variable  used  as  index  for  outer  DO  loop  which 
treats  each  line  in  turn. 

Integer  variable  used  as  index  for  inner  DO  loop  which 
assigns  spectral  contribution  of  each  line  to  output 
wavelength. 

Integer  variable  used  as  final  parameter  in  inner  DO 
loop,  index  of  greatest  output  wavelength  for  a given 
line. 

Integer  variable  used  as  initial  parameter  in  inner 
DO  loop,  index  of  shortest  output  wavelength  for  a 
given  line. 


133 


OOOOOOOCJOOOOOOOOC5000 
HMKJ’lAiDNetT'OWWMJ'IAiOK  05  O'  n 
W H W r4  H H H rH  rt  m 


or 

o 

O' 

X 

in 

-J 

1 

— Z 

*— ■ 

Z UJ 

X 

T. 

• > 

c 

ID 

X c 

z 

2 

=>  3 

■> 

a: 

Z 

3C 

M 

w 

■>  o 

r 

fy 

l/) 

S 3 

a 

V 

CD 

- 1- 

U 

-» 

■H 

< 

O 

in 

* 

CD 

o 

C* 

«r  •« 

-*■ 

UJ 

• 

• 

CO 

o 

O 1 

ck  z 

*rr 

o 

CVj 

W 

in 

• 

•>o 

o 

• 

H 

o c« 

O M 

<r 

rl 

4 

4* 

II 

o • 

UJ  M 

a 

H 

O 

in 

a o 

•» 

II 

in 

IT 

Ll 

►— 

00  _J 

• 

• 

II  *D 

4 

w o 

H CL 

M 

o 

O 

o 

4 

V 

n t 2 

i uj  o ~ 

0.0  » M 
l/)0.  rl- 
00  X 
Zhii  3 
C V z 

M *-  •> 
1/1  7 •» 

Z O o 
li  T IT  II 
T T rl 
MO  M 
C O O O' 
O 


in  in  ^ 
n O'  rl 

J-  fo 


I I t—  * — »— 


■H  CO  "5 

O U.  * — 
« ll  T r>  o 
rH  ► rl  UJ 

lt  a 
• • “5  I 00 


_l  09  09  II  II 

• • "r  UJ  UU 

•0  r,  25  T 

U-  V II  - 7 I 57 
— ) "3  o w i—  I— ; Oi 
v ^ j-  (.'  I-  1-0 

rl  3 UJ  7 Z I-  Cl 

u_  u_  z a o c uj  z 

hmc  arcnooc^uj 

Cl 


0.  O' 

W W C/' 
— i 

ii  ii  ~ 


U.  V 
~>  -5  c. 


C/J  u_  u. 
“5  -5  M 


II  iillll  I III  II U,  PMPi^tWV  m « 


iSPP 


PMMip 


F 


I. 

' 

i 

I 


i 

: 

I 


: 

r 

: 


Subroutine  ROTATE 

Purpose:  Subroutine  ROTATE  is  used  to  read  in  spectral  line 

data  from  the  compilation  of  McClatchey,  et  al.  The 
input  data  consists  of  line  positions  and  strengths 
for  vibrational  - rotational  bands  of  water,  carbon 
dioxide,  ozone  and  nitrous  oxide,  and  for  the  pure 
rotational  spectrum  of  water.  This  data  is  processed 
for  use  by  other  subroutines. 

Program  Flow:  The  DO  range  ending  with  statement  110  reads  in 

data  on  the  strongest  lines  in  a vibration  - rotation 
band.  This  data  ''.onsists  of  the  wave  number  of  the 
line,  the  line  strength  at  296  K and  the  wavenumber 
of  the  lower  rotational  state.  ROTATE  converts  the 
wavenumber  to  wavelength  in  //m  stored  in  array 
HHOLAM,  converts  the  line  strength  to  units  cm”1 /atm  - 
cm  at  273.  15k  stored  in  array  HHOSTR  and  stores  the 
wavenumber  of  the  lower  state  in  array  HHOLSE. 

The  rotational  partition  function  for  water  vapor 
is  computed  next  if  water  vapor  data  is  being  read. 

The  DO  range  ending  at  statement  130  is  used  only 
for  water  vapor  pure  rotational  lines  and  computes 
the  wavelength  (stored  in  array  ROTLAM),  Einstein 
A coefficient  (stored  in  array  EA),  relative  line 
strength  in  absorption  at  273. 15k  (stored  in  array 
ROTST),  relative  line  strength  in  emission  at 
273.  15k  (stored  in  array  ROTEA)  and  the  wave- 
number  of  the  lower  rotational  state  (stored  in  array 
ROTLSE). 


135 


Mnemonic  Variables : 


EA 

EALC 

EALR 

ENERGY 

ENL 

ENY 

HHOLAM 

HHOLSE 

HHOSTR 

ROTEA 

ROTLAM 

ROTLSE 

ROTST 


R«al  variable  used  to  store  Einstein  A coefficient 
for  water  vapor  pure  rotational  spectrum  line. 

Real  variable  used  to  store  intermediate  results  in 
computing  pure  rotational  line  strengths  in 
emission. 

Real  variable  used  to  store  intermediate  results 
in  computing  pure  rotational  line  strengths  in 
absorption. 

Real  variable  into  which  is  read  wavenumbers  of 
pure  rotational  lines. 

Real  variable  into  which  is  read  wavenumbers  of 
lower  states  of  pure  rotational  lines. 

Real  variable  into  which  is  read  wavenumber  of 
upper  state  of  pure  rotational  line. 

Real  array  in  which  wavelengths  of  vibration  - 
rotation  band  lines  are  stored. 

Real  array  in  which  wavenumbers  of  lower  states 

of  a vibration  - rotation  band  are  stored. 

Real  array  in  which  line  strengths  of  vibration  - 

-2  1 

rotation  band  are  stored  (cm  - /XT m~  STP). 
Real  array  in  which  pure  rotational  line  strengths  in 
emission  (273.  15k)  are  stored. 

Real  array  in  which  pure  rotational  line  wave- 
lengths are  stored. 

Real  array  in  which  lower  state  wavenumber  of 
pure  rotational  lines  are  stored. 

Real  array  in  which  pure  rotational  line  strengths 
in  absorption  are  stored. 
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WVLTH 

EAFAC 

QR 

QRFAC 

RA 

RB 

RC 

TEFAC 

TQFAC 

I 

NLINES 


Real  array  in  which  wavelengths  of  vibration  - 
rotation  band  lines  are  stored  (redundant). 

Real  variable  used  to  store  pre-factor  in  computing 
Einstein  A coefficients  for  pure  rotational  lines. 

Real  variable  in  which  is  stored  approximate 
rotational  partition  function  for  water  vapor 
molecule . 

Real  variable  into  which  is  read  prefactor  for 
computing  rotational  partition  function. 

Real  variable  into  which  is  stored  first  water  vapor 
rotational  constant. 

Real  variable  into  which  is  stored  second  water  vapor 
rotational  constant. 

Real  variable  into  which  is  stored  third  water  vapor 
rotational  constant. 

Real  variable  in  which  is  stored  exponential  factor 
used  in  converting  line  intensities  at  296K  to  line 
strengths  at  273.  15K. 

Real  variable  in  which  is  stored  factor  relating 
rotational  function  at  296K  to  that  at  273.  15K. 

Integer  variable  used  as  index  in  DO  loop. 

Integer  variable  into  which  is  read  number  of  lines 
included  in  pure  rotational  spectrum  and  in  any 
single  vibration  - rotation  band. 
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Subroutine  WATSTR 


Purpose:  WATSTR  computes  the  radiances  of  individual  lines 

in  vibration  - rotation  bands  and  in  the  pure  rotational 
spectrum  of  water  vapor.  The  radiances  are  computed 
with  the  assumption  that  the  atmosphere  in  the  line  of 
sight  may  be  treated  as  having  uniform  temperature. 
The  radiances  computed  for  the  vibration  - rotation 
bands  are  scaled  to  match  the  total  band  radiances 
produced  by  the  BACKGROUND  program.  Absolute 
radiance  values  are  computed  for  the  water  vapor  pure 
rotational  spectrum. 

Program  Flow:  Upon  entry  to  WATSTR,  a conversion  factor  TQFAC 

is  computed  for  later  use  in  computing  the  effect  of 
charging  rotational  partition  function  with  temperature. 
Control  is  then  transferred  to  the  appropriate  part  of 
the  subroutine  as  determined  by  the  value  of  N. 

For  computation  of  the  radiances  of  lines  in  the  water 
vapor  pure  rotational  spectrum  (N  = 1),  control  is 
transferred  to  statement  100.  A radiance  factor 
RADFAC  and  an  optical  thickness  factor  TAUFAC 
are  computed.  Then  the  optical  thickness  and 
radiance  values  corrected  for  optical  thickness 
and  induced  emission  are  computed  for  each  line  in 
the  DO  range  ending  with  statement  110.  Control  is 
returned  to  the  main  program.  Line  constants 
generated  by  subroutine  ROTATE  are  available  in  the 
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arrays  of  common  area  ROTCOM.  Radiance  values 
are  returned  in  array  HRAD. 


Mnemo  nic 
RAD 

RM 

T 

DEN 

N 


Control  is  transferred  to  statement  200  for  compu- 
tation of  radiances  of  lines  in  the  other  bands.  An 
optical  thickness  factor  TAUFAC  and  factor  TEFAC 
used  in  computing  populations  of  rotational  states  at 
temperature  T are  computed.  The  variable  SUM, 
used  to  accumulate  the  sum  of  relative  radiances,  is 
set  to  zero.  In  the  DO  range  ending  with  statement 
210,  the  optical  thickness  of  each  line  is  computed 
and  the  relative  radiance  value,  corrected  for  optical 
thickness,  is  computed  and  stored  in  array  HRAD. 

The  sum  relative  radiances  is  accumulated  in  SUM. 
Upon  exit  from  this  DO  range,  the  total  radiance  RAD 
is  divided  by  SUM  and  stored  in  RADFAC.  In  the  DO 
range  ending  with  statement  220,  absolute  radiance 
values  are  obtained  and  stored  in  array  HRAD. 

Control  is  then  returned  to  the  main  program. 

Variables: 

Real  variable  containing  total  band  strength  for 
vibration  - rotation  band. 

Real  variable  containing  molecular  weight  of  molecule. 
Real  variable  containing  value  of  temperature. 

Real  variable  containing  total  number  of  ground  state 
molecules  per  square  centimeter  in  line  of  sight. 
Integer  variable  determining  whether  entry  to  WATSTR 
results  in  computation  for  water  vapor  pure  rotational 
spectrum  (N  = 1)  or  other  band  (N  = 2). 
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HHOLAM 

HHOLSE 

HHOSTR 

HRAD 

ROTEA 

ROTLAM 

ROTLSE 

ROTST 

HSTR 

RADFAC 

SUM 

TAU 

TAUFAC 

TEFAC 

TQFAC 


Real  array,  values  as  in  ROTATE. 

Real  array,  values  as  in  ROTATE. 

Real  array,  values  as  in  ROTATE. 

Real  array  used  to  return  radiance  values  of 
individual  lines. 

Real  array,  values  as  in  ROTATE. 

Real  array,  values  as  in  ROTATE. 

Real  array,  values  as  in  ROTATE. 

Real  array,  values  as  in  ROTATE. 

Real  variable  used  as  temporary  storage  in  computing 
optical  thickness. 

Real  variable  in  which  is  stored  prefactor  used  in 
computing  radiance  values. 

Real  variable  in  which  partial  sums  of  relative 
radiances  are  accumulated. 

Real  variable  in  which  optical  thickness  at  line  center 
is  generated. 

Real  variable  in  which  prefactor  for  optical  thickness 
is  stored. 

R®al  variable  in  which  factor  used  in  computing 
rotational  populations  is  stored. 

Real  variable  in  which  factor  converting  rotational 
partition  function  is  stored. 
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Subroutine  ROTSTR 


Purpose:  Subroutine  ROTSTR  computes  the  relative  radiances 

of  lines  in  the  vibration  - rotation  bands  of  those 
linear  molecules,  which  are  not  included  in  the  line 
parameter  compilation  of  McClatchey,  et  al.  The 
radiances  are  corrected  for  effects  of  optical  thick- 
ness, and  absolute  values  then  determined  by 
normalizing  to  total  band  radiances  computed  by  the 
BACKGROUND  program.  As  for  nonlinear  molecules, 
a single  temperature  is  assumed  for  the  molecules  in 
the  line  of  sight. 


Program  Flow: 


Upon  entry  to  ROTSTR,  the  reduced  rotational 
temperature  SIG  is  computed,  relative  strengths  of 
the  P,  Q and  R lines  for  J = 0 are  assigned  and 
values  assigned  to  the  increment  of  J values  DJ 
and  the  variable  J value  AJ.  Depending  on  the  value 
of  ND,  relative  strengths  of  rotational  lines  are 
computed  in  either  the  DO  range  ending  with  state- 
ment 110  (parallel  transition)  or  the  DO  range  ending 


with  statement  120  (perpendicular  transition).  Control 
is  transferred  to  statement  300  and  the  prefactor 
TAUFAC  used  in  computing  optical  thickness  is 
evaluated.  The  variable  SUM  which  accumulates  the 
sum  of  relative  radiances  is  set  to  zero.  In  the  DO 
range  ending  with  statement  310,  optical  thicknesses 
for  lines  in  the  P,  Q and  R branches  and  relative 
radiances  for  these  lines,  corrected  for  the  effects 
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of  optical  thickness,  are  computed.  In  the  DO  range 
ending  with  statement  320,  absolute  line  radiances 
normalized  to  the  total  band  radiance  RAD  are  computed. 
Control  is  then  returned  to  the  main  program. 


Mnemonic  Variables: 


RAD 

RM 

BVPP 

AJ 

DJ 


DEN 


Real  variable  containing  total  band  strength. 

Real  variable  containing  molecular  weight  of  molecule. 
Real  variable  containing  mean  rotational  constant  of 
lower  vibrational  state. 

Real  variable  containing  varying  value  of  rotational 
quantum  number  J of  lower  vibrational  state. 

Real  variable  containing  value  of  increment  of 
rotational  quantum  number. 

Real  variable  containing  mean  temperature  of  emitting 
layer. 

Real  variable  containing  number  of  molecules  in  lower 


STR 

PS 


state  per  square  centimeter  in  line  of  sight. 

Real  variable  containing  band  strength  at  300K. 

Real  array  in  which  radiances  of  lines  in  P branch 
are  returned. 


QS 

RS 

NP 


Real  array  in  which  radiances  of  lines  in  Q branch 
are  returned. 

Real  array  in  which  radiances  of  lines  in  R branch 
are  returned. 

Integer  variable  determining  whether  band  is  parallel 
(NP  = 1)  or  perpendicular  )NP  = 2). 
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ND 

P 

Q 

R 

FAC 

PTAU 

QTAU 

RADFAC 

RTAU 

SIG 

SUM 

TAUFAC 


Integer  variable  determining  whether  all  lines  are 
present  (ND  = 1)  or  alternate  lines  are  missing 
(ND  = 2). 

Real  array  containing  wavelengths  of  lines  in  P branch. 
Real  array  containing  wavelengths  of  lines  in  Q branch. 
Real  array  containing  wavelengths  of  lines  in  R branch. 
Real  variable  used  to  store  exponential  factor  in 
computation  of  line  strengths. 

Real  variable  storing  optical  thickness  in  a P branch 
line. 

Real  variable  storing  optical  thickness  in  a Q branch 
line. 

Real  variable  storing  ratio  of  total  band  radiance  RAD 
and  sum  of  relative  line  radiances. 

Real  variable  storing  optical  thickness  of  R branch 
line. 

Real  variable  storing  reduced  rotational  temperature. 
Real  variable  which  accumulates  sum  of  relative  line 
radiances. 

Real  variable  storing  prefactor  used  in  computing 
optical  thickness  at  line  centers. 
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Subroutine  ROTPOS 


Purpose: 


Method : 


ROTPOS  computes  the  wavelength  of  lines  in  the  P, 
Q and  R branches  of  the  infrared  vibrational  - 
rotational  bands  of  linear  molecules  neglecting  spin 
and  A orAdoubling.  It  is  used  for  bands  not  in- 
cluded in  the  compilation  of  McClatchey,  et  al. 

The  line  positions  are  computed  using  equations 
(IV,  19),  (IV,  20)  and  (IV,  22)  of  Herzberg  (1945). 

It  should  be  noted  that  use  of  these  simple  formulae 
results  is  an  error  in  line  positions  which  may 
approach  0.  l//m  for  large  J values  in  transitions 
between  doublet  states. 


Program  Flow:  Upon  entry  to  ROTPOS,  the  auxiliary  values  BDIF, 

BSUM,  RCENT  and  TRI  used  in  the  compution  are 
computed  from  the  band  center  position  and  the 
rotational  constants.  The  initial  J value  is  set  to 
zero  and  provisions  made  for  omitting  alternate 
lines  by  setting  DJ  equal  to  ND.  This  latter  starting 
condition  means  that  the  alternate  lines  omitted  are 
always  those  with  odd  J values.  The  wavelengths  of 
100  lines  ir  each  branch  in  units  j-L  m are  then 
computed  in  the  DO  range  ending  with  statement  110. 
The  wavelengths  are  printed  out  and  control  returned 
to  the  main  program. 
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Mnemonic  Variable s : 


AJ 

BDIF 

BSUM 

DJ 

RCENT 

TRI 

CENTER 

BVP 

BVPP 

P 

Q 

R 

ND 


Real  variable  in  which  current  rotational  quantum 
number  J is  stored. 

Real  variable  storing  difference  of  rotational  constants 
BVP  and  BVPP. 

Real  variable  storing  sum  of  rotational  constants. 

Real  variable  storing  increment  of  rotational  quantum 
number. 

Real  variable  containing  constant  used  in  computing 
R branch  line  positions. 

Real  variable  containing  constant  ubed  in  computing 
R branch  line  positions. 

Real  variable  giving  the  wave  number  of  the  band 
center. 

Real  variable  giving  the  mean  rotational  constant 
of  the  upper  vibrational  state. 

Real  variable  giving  the  mean  rotational  constant 
of  the  lower  vibrational  state. 

Real  array  in  which  wavelengths  of  the  P branch 
are  returned. 

Real  array  in  which  wavelengths  of  the  Q branch  are 
returned. 

Real  array  in  which  wavelengths  of  the  R branch  are 
returned. 

Integer  variable  determining  whether  all  lines  in  a 
band  are  present  (ND  = 1)  or  whether  alternate  lines 
are  missing  (ND  = Z). 
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PLOTTING  PROGRAM 


The  plotting  program,  TRY  CRT,  which  was  developed  under 
the  previous  contract,  has  been  revised  and  adapted  to  use  the  pen 
and  paper  plotter  in  addition  to  the  microform  plotter.  The  version 
presented  here  is  designed  to  read  the  output  file  produced  by  the 
spectral  radiance  program,  SPCTRA.  It  consists  of  the  main  program, 
TRY CRT,  and  six  subroutines,  LINSCA,  LOGSCA,  LINAXS,  LOGAXS, 
GRID  and  NULINE.  In  addition,  it  uses  the  AFCRL/CDC  6600  system 
subroutines  which  are  needed  to  produce  plots. 

Array  IDENT  contains  information  required  by  the  computer 
installation  to  identify  the  plotted  graphs.  The  CDC  6600  subroutines 
are  initialized  by  calls  to  subroutines  PLTID3  and  PLOT.  The 
number  of  graphs  to  be  plotted,  NCASES,  number  of  characters,  NXL, 
in  x-axis  caption,  x-axis  caption,  XBCD,  number  of  characters,  NYL, 
in  y-axis  caption,  and  y-axis  caption,  Y BCD,  are  read  in  from  cards 
and  printed  on  the  line  printer.  The  DO  loop  ending  with  statement 
400  is  then  executed  for  each  graph  to  be  produced.  Minimum  wave- 
length, WVMIN,  and  maximum  wavelength,  WVMAX,  are  read. 

Plotting  constants  INF,  NX  and  NY  are  computed.  Wavelengths  are 
from  the  SPCTRA  output  tape  (channel  2)  and  stored  in  array 
AVLAM.  An  appropriate  subset  is  then  stored  in  array  X.  These  are 
then  scaled  by  subroutine  LINSCA.  Radiance  data  are  then  read  from 
channel  2 and  stored  in  array  SPCTRA.  An  appropriate  subset  is 
selected  and  stored  in  array  Y.  Values  in  Y are  scaled  for  plotting 
by  a call  to  subroutine  LOGSCA.  Axes  are  drawn  on  all  four  sides 
of  the  graph  by  calls  to  subroutines  LINAXS  and  LOGAXS.  Calls  to 
subroutine  GRID  place  vertical  and  horizontal  grid  lines  on  the 
graph.  Radiance  values  are  further  scaled  in  the  DO  loop  ending 
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with  statement  150.  Radiance  values  are  placed  on  the  graph  by  calls 
to  subroutine  NULINE  in  the  loop  ending  with  statement  200.  After 
each  set  of  curves  on  a graph  is  produced,  a call  to  subroutine  PLOT 
prepares  for  a new  graph.  The  DO  loop  ending  with  statement  250  is 
repeated  three  times,  the  first  and  second  cycles  plotting  vertical 
and  horizontal  endoatmospheric  radiances  while  the  third  plots  limb 
viewing  radiances.  When  only  limb  viewing  radiances  are  desired, 
statement  250  should  be  placed  immediately  after  the  statement 
"READ  (2)  SPCTRA"  on  line  TORT  540. 

Subroutines  LINSCA  and  LOGSCA  were  written  to  provide 
automatic  scaling  of  data  points  in  a manner  which  would  utilize 
maximum  space  on  the  graph.  In  addition  to  scaling  data  points, 
these  subroutines  generate  values  KONX,  DX,  MINY  and  MAXY  which 
provide  data  which  is  used  in  the  axis  drawing  subroutines  LINAXS 
and  LOGAXS  and  the  grid  drawing  subroutine  GRID.  The  subroutine 
NULINE  was  written  to  replace  the  standard  AFCRL/CDC  6600  sub- 
routine LINE  which  draws  straight  lines  between  plots.  NULINE 
permits  points  to  be  plotted  a limited  distance  outside  the  graph 
boundaries  if  the  data  for  a graph  covers  more  than  9 orders  of 
magnitude. 

The  following  listing  includes  both  the  original  main  program 
using  microfilm  plotter  and  a main  program  TRYPEN  which  uses  the 
pen  plotter. 
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Program  BCKGND  (BGND) 


Purpose: 


Program 

Flow: 


The  main  program  BCKGND  (short  for  BACKGROUND) 
serves  primarily  to  control  the  operation  of  subroutines 
which  compute  vibrational  populations  and  radiative  rates. 
BCKGND  reads  and  checks  the  input  data,  and  transfers 
control  to  Subroutine  BANRAD  which  computes  vibrational 
populations,  local  optical  thicknesses  and  volume  radi- 
ances for  all  included  bands  of  a single  molecular  species. 
When  control  is  returned,  BCKGND  computes  and  outputs 
integrated  band  radiances  for  limb  viewing  (exoatmospheric) 
and  seven  endoatmcspheric  viewing  angles  at  selected 
tangent  heights  or  endoatmo spheric  altitudes. 

The  initial  step  in  operation  of  BCKGND  is  the  reading  of 
all  input  data.  The  number  IGAS  of  molecular  species 
including  total  number  density  is  read  in.  Then  in  the  DO 
range  ending  with  statement  170,  the  following  operations 
are  performed.  The  title  card,  containing  alphanumeric 
identification,  an  index  check  number  and  maximum  alti- 
tude considered  for  the  species  is  read  in  and  its  check 
number  checked  against  the  index  of  the  DC'  range.  If  the 
check  succeeds,  altitude  and  number  densities  are  read 
into  arrays  ALT  and  STORE  until  an  altitude  correspond- 
ing to  the  maximum  altitude  is  found.  When  this  altitude 
is  found,  the  logarithm  of  input  number  densities  is  formed 
in  the  DO  range  ending  with  statement  130  and  control  is 
passed  to  Subroutine  INTERP  which  interpolates  to  form 
values  at  1 km  intervals  up  to  160  km  and  at  5 rm  intervals 
at  higher  altitudes.  The  interpolated  values  are  converted 
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to  number  densities  in  the  DO  range  ending  with  state- 
ment 160.  In  the  DO  range  ending  with  statement  175, 
total  atmosphere  number  densities  are  stored  in  array 
TOT  DEN.  Temperature  data  are  then  read  in  similarly 
in  the  DO  range  ending  with  statement  245  and  interpo- 
lated values  found  by  Subroutine  INTERP.  The  complete 
set  of  number  densities  and  temperatures  is  then  pointed 
out.  Next,  vibrational  level  and  band  constants  are  read 
in  for  each  species  in  the  DO  range  ending  with  statement 
331.  The  first  card  for  each  species  contains  an  alpha- 
numeric identification,  number  of  vibrational  levels  NLEV, 
number  of  bands  NBAN,  molecular  mass  AMAS,  and  two 
vibrational  energy  transfer  coefficients  DEXT  and  DEXB. 
The  alphabetic  code  is  checked  and  if  it  is  a valid  code 
stored  in  array  DATATY,  these  constants  are  transferred 
to  the  appropriate  arrays.  In  the  DO  range  ending  with 
statement  313,  constants  for  each  vibrational  level  are 
read  in,  and  card  sequence  numbers  checked.  In  the  DO 
range  ending  with  statement  316,  constants  for  each  band 
are  read  in  and  card  sequence  numbers  checked.  For 
^^2  anc^  ^2^'  tabular  values  of  temperatures  and  associ- 
ated vibrational  energy  exchange  coefficients  are  read 
before  the  vibrational  level  and  molecular  band  data. 

The  computation  of  molecular  band  radiances  then  pro- 
ceeds in  the  following  manner.  The  outer  DO  range 
index  NLTE  and  ending  with  statement  900  controls  the 
factors  considered  in  a complete  set  of  calculations  for 
any  of  five  possible  sets.  The  corresponding  DO  state- 
ment may  be  varied  depending  on  what  factors  are  desired. 
See  the  description  of  Subroutine  BANRAD  for  the  effect 


of  different  values  of  the  index  NLTE.  Each  cycle 
through  the  n ited  DO  range  with  Index  I and  ending 
with  statement  50  results  in  the  computation  of  band 
radiance  values  >r  a single  molecule.  Action  within 
this  DO  range  is  s fellows. 

Subroutine  BANRAD  is  called  and  returns  volume 
emission  radiance  values  in  array  RAD  and  local 
optical  thickness  values  in  array  TAU.  Then  for  each 
band,  line  of  sight  band  radiances  are  computed  in  the 
DO  range  ending  with  statement  840.  Within  this  DO 
range,  the  value  of  NLTE  is  examined  and.  appropriate 
headings  are  printed.  Then  trapezoidal  rule  integration 
is  performed  to  determine  the  radiance  arising  along 
each  of  the  eight  lines  of  sight,  including  the  modifica- 
tion necessary  for  the  effects  of  optical  thickness  along 
the  line  of  sight.  The  DO  range  ending  with  statement 
830  controls  the  integration  for  each  tangent  height  and 
endoatmc spheric  viewing  level.  In  the  current  version, 
radiance  values  are  computed  at  5 km  intervals  of  tan- 
gent height  up  to  150  km  and  at  25  km  intervals  above 
this  level,  ^he  integration  of  radiance  values  is  per- 
formed for  jach  endoatmospheric  viewing  level  in  the  DO 
range  ending  with  statement  810.  The  integration  begins 
at  the  viewing  altitude.  At  each  point  on  the  line  of  sight, 
the  contribution  of  the  local  volume  emission  -ate  to  that 
received  at  the  viewing  point  is  computed  by  Subroutine 
SMI,  which  takes  into  account  the  effects  of  optical  thick- 
ness between  the  emitting  point  and  the  viewing  point. 

The  equivalent  integration  for  the  limb  viewing  case  is 
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performed  in  the  DO  ranges  ending  with  statements 
815  and  818.  Upon  completion  of  the  integration  for 
each  tangent  height  and  endoatmospheric  viewing  level, 
band  radiances  for  that  level  are  punched  and  printed. 
This  punched  output,  together  with  that  pioduced  by 
Subroutine  BANRAD,  is  used  as  input  data  to  the  special 
radiance  program  SPECTRA. 

Tables  V-2  through  V-4  list  the  molecular  bands  in 
the  program.  These  are  follo\,  J d by  a sample  set  of 
input  data  and  the  program  listing. 
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{ 

Table  V-2 


Water  Vapor 

Band  Data 

Vibrational 

Transition 

Band  Center 

Band  Strength 

Solar  Flux 

Upper 

State 

Lower 
Sta  te 

Wavelength 

{Um) 

at  296  K 

1 2 - 1 
(mol-1  cm  cm  ) 

Band  Center 
(photons  cm“^//m-l) 

010 

000 

6.27 

1.  06(-17)* 

5.60(15) 

020 

000 

3.17 

6. 5 8 ( - 20 ) 

3. 50(16) 

100 

000 

2.  73 

3. 62(-19) 

5.  10(16) 

001 

000 

2.  66 

CO 

r-H 

O' 

O' 

r- 

5.40(16) 

020 

010 

6.42 

8.  61 (-21 ) 

5. 70(15) 

100 

010 

4.85 

8.  90(-23) 

1. 20(16) 

001 

101 

4.63 

4.  90(-22) 

1. 50(16) 

030 

010 

3.  26 

7.  99(-23) 

3. 30(16) 

110 

010 

2.75 

1.  50(-22) 

5.  00(16) 

Oil 

010 

2.  68 

3. 30(-21 ) 

5.40(16) 

030 

020 

6.  60 

7. 00(-24) 

4. 90(15) 

030 

000 

2.  14 

2.  00(-22) 

9. 30(16) 

no 

000 

1.91 

1. 83(-20) 

1. 20(17) 

oil 

000 

1.88 

9.  1 6 ( - 1 9 ) 

1.  25(17) 

* Number  in  parentheses  is  power  of  ten. 
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Table  V-3 


Carbon  Dioxide  Band  Data 

Vibrational 


Transition 

Upper  Lower 
State  State 

Band  Center 

Wavelength 
(Mm) 

Band  Strength 
at  296  K 

/ ■!  “ 1 2 -1 
(mol  cm  cm  ) 

Solar  Flux 
Band  Center 
-2 

(photons  cm  iJm 

01101 

00001 

14.  98 

8.  26  (-18) 



4.  60  (14) 

10002 

01101 

16.  18 

1.  44  (-19) 

3.  70  (14) 

02201 

01101 

14.  98 

6.  49  (-19) 

4.  60  (14) 

11102 

10002 

15.45 

2.  22  (-20) 

4.  20  (14) 

1102 

02201 

16.  74 

5.  21  (-21) 

3.  35  (14) 

03301 

02201 

14.  97 

3.  82  (-20) 

4.  60  (14) 

10001 

01101 

13.  87 

1.  85  (-19) 

5.  80  (14) 

00011 

10001 

10.41 

4.  91  (-22) 

1.  31  (15) 

00011 

00001 

4.  62 

9.  60  (-17) 

1.  60  (16) 
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Table  V-4 


Ozone  Band  Data 


Vibrational 


Transition 

Upper  Lower 

State  ' State 

Band  Center 

Wavelength 
(U.  m) 

Band  Strength,  296  K 
(mol- 1 cm^  cm"') 

Solar  Flux 
(photons  crri  / 

010 

000 

14.  27 

4.  1 3(- 1 9)* 

5.  30(14) 

001 

000 

9.  60 

1 . 29  ( - 1 7) 

1. 70(15) 

100 

000 

9.  06 

3.  47 (-19) 

1.90(15) 

101 

000 

4.  74 

1.  33(- 1 8) 

1.20(16) 

111 

000 

3.  59 

2.  32(-20) 

2. 50(16) 

003 

000 

3.  29 

1.  10(-19) 

3.20(16) 

101 

100 

9.92 

9.  49(-21 ) 

1.50(15) 

011 

010 

9.79 

4.  23(-l  9) 

1. 60(15) 

002 

001 

9.74 

1. 62(-l  9) 

1.60(15) 

003 

002 

10.  29 

1. 66(-21) 

1.40(14) 

011 

000 

5.  81 

3. 1 C(-20) 

6.  90(15) 

* Number  in  parentheses  is  power  of  ten. 
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120  CONTINUE 


Subroutine  BANRAD 


Purpose: 


Method: 


BANRAD  computes  the  population  of  vibrational  levels 
of  infrared  active  species,  including  the  effects  of  col- 
lisional  excitation;  absorption  of  radiation  from  the 
lower  atmosphere,  sun  and  molecules  radiating  at 
high  altitude.  The  values  returned  to  the  main  pro- 
gram are  volume  radiance  and  unit  optical  thickness 
at  each  altitude  for  each  band  included  in  the  computation. 

Subroutine  BANRAD  computes  the  population  of  vibra- 
tional levels  of  infrared  active  species,  including  the 
effects  of  collisional  excitation  and  absorption  radiation 
from  the  lower  atmosphere,  sun  and  molecules  radiat- 
ing at  high  altitude.  The  values  returned  to  the  main 
program  are  volume  radiance  and  unit  optical  thickness 
at  each  altitude  for  each  band  included  in  the  computat.^n. 

The  collisional  processes  which  are  assumed  to  effect 
the  vibrational  state  of  a molecule  are  (1)  tranlational- 
vibrational  (T-V)  interactions  in  which  translation 
energy  is  transformed  into  vibrational  energy  or  vice- 
versa  during  a collision,  (2)  intermolecular  vibrational 
exchange  (V-V)  in  which  a molecule  of  one  species  ex- 
changes a quantum  of  vibrational  energy  with  a molecule 
of  another  species,  and  (3)  intramolecular  vibrational 
exchange  in  polyatomic  molecules  in  which  a collision 
transfers  the  energy  in  one  mode  of  vibration  into 
another  mode  within  the  same  molecule. 
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When  a molecule  in  the  ground  state  collides  with 
another  molecule  there  is  a finite  probability  P^, 
that  a portion  of  the  relative  translational  energy  will 
be  transformed  to  vibrational  energy  if  the  translational 
energy  is  equal  to  or  exceeds  the  vibrational  energy. 
Similarly,  if  the  molecule  is  in  an  excited  vibrational 
state,  there  is  a finite  probability  P ^ that  the  vibra- 
tional energy  will  be  converted  to  relative  translational 
energy  during  the  collision.  For  a Boltzmann  distribu- 
tion of  translational  energies  at  a given  temperature 

■n  T3  gu  . he  , 

PTV  PVT  g GXP  ( " X kT  ) 

1 o 

where  g^  and  gj  are  the  statistical  weights  of  the  upper 

and  lov/er  vibrational  levels,  h is  Planck's  constant 
-27 

(6.6256  x 10  erg-sec),  c is  the  velocity  of  light 
(2.9979  x 10^  cm/sec),  k is  the  Boltzmann  constant 
(1.  38054  x 10  ^ erg/K),  T is  absolute  temperature  (K) 
and  XQ  is  the  wavelength  (cm)  corresponding  to  the 
energy  difference  of  the  two  vibrational  levels.  Similar 
expressions  hold  for  inter  molecular  and  intramolecular 
Y-V  transfer. 

The  rate  at  which  a vibrational  level  of  a single  mole- 
cule is  excited  (K^,^.)  or  de -excited  (K^.^,)  is  given  by 
the  product  of  the  collision  frequency  Z and  the  trans- 
ition probability.  The  subroutine  uses  temperature 
dependent  rate  coefficients  k^y,  etc. , such  that 

if  [M]  is  the  concentration  of  the  collision  partner, 

^yp  — ^yp[^^]>  etc. 


201 


The  radiative  processes  which  are  assumed  to  affect 
the  vibrational  state  of  a molecule  are  (1)  absorption 
of  solar  radiation,  (2)  absorption  of  radiation  from  the 
lower  atmosphere,  (3)  absorption  of  radiation  emitted 
by  molecules  in  the  upper  atmosphere  and  (4)  spontane- 
ous emission  at  a rate  A^^(sec  *),  the  Einstein  A co- 
efficient. Because  absorbing  molecules  are  present 
between  each  source  of  radiation  and  the  molecule  for 
which  the  vibrational  state  is  to  be  computed,  radiative 
transfer  functions  must  be  used  to  take  into  account  the 
effects  of  absorption.  The  model  assumes  that  individ- 
ual lines  have  a Doppler  (thermal  velocity  broadened) 
spectral  contour. 

3 

At  a given  altitude,  there  will  be  Nj  molecules /cm 

in  the  ground  vibrational  state  of  a molecule  and 
3 

molecules/cm  in  an  upper  state.  Molecular  nitro- 
gen and  oxygen  are  assumed  to  be  the  only  collision 
partners  effective  in  transfer  of  vibrational  quanta,  and 
their  concentrations  are  used  to  determine  total  vibra- 
tional excitation  rates  KTy,  KVT,  etc.  The  treatment 
below  refers  to  a two  state  molecule. 

In  unit  volume,  molecules  in  the  ground  state  are 
excited  to  the  upper  Gtate  by: 

(1)  Collisional  T-V  excitation  at  a rate 

(2)  Collisional  V-V  excitation  at  a rate  N.K.  ... 

1 MV 

(the  subscript  MV  is  used  for  excitation,  VM 
for  de -excitation,  by  transfer  of  vibrational 
quanta). 


(3) 


Absorption  of  solar  flux  at  a rate  N I G'M(k  ) 

1 o I s 

where  G'  is  molecular  cross  section  for  ab- 


sorption of  radiative  flux 


(4)  Absorption  of  lower  atmospheric  radiance  at 

a rate  N.Nv  G'N^  (k) 

A A)  t 1 6 t 


(5)  Absorption  of  radiance  emitted  at  other  levels 

of  the  atmosphere  at  a rate  N f Nv  G'Nb  (k  ) dz 

1/  A , D 21  D 


Molecules  in  the  upper  state  are  de-excited  by: 


(1)  Colli sional  V-T  de- excitation  at  a rate  N K 

u VT 


(2)  Colli  sional  V-V  de- excitation  at  a rate  N 


(3)  Spontaneous  radiative  de -excitation  at  a rate 


N A , 
u ul 


Excitation  rates  are  assumed  to  balance  de -excitation 
rates  and  the  equation 


dN 

"dT  = Nl(kTV+KMV  + 1c»C,M^ks>  + NX,tG'NU(kt) 


+/NX,DG'N21<kD>dZ  * N1<KVT+KVM+Aul>  * ° 


is  solved  by  iteration  to  determine  N and  N at  all 
7 u 1 


altitudes. 
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Program 

Flow: 


The  first  step  in  the  computation  is  determination  of 
the  thermal  equilibrium  populations  of  vibrational 
levels.  This  is  performed  in  the  DO  range  ending  with 
statement  130.  Each  passage  through  this  DO  range 
computes  populations  for  one  altitude.  Within  this  DO 
range  are  two  others.  The  DO  range  ending  with  state- 
ment 110  computes  the  vibrational  pop\ilation  of  each 
level  relative  to  the  ground  vibrational  level,  using  the 
energies  (wave-number)  stored  in  WAVE  ant  the  statis- 
tical weights  stored  temporarily  in  real  array  VIBPOP. 
The  sum  of  relative  populations  is  stored  in  SUM.  The 
final  ground  state  population  is  then  computed  by  divid- 
ing total  concentration  by  SUM,  and  stored  in  real  array 
GNDSTE.  The  DO  range  ending  with  statement  120  then 
computes  the  vibrational  population  of  each  level  by  multi- 
plying total  concentration  by  relative  concentration 
divided  by  SUM. 

The  DO  range  ending  with  statement  135  selects  the 
vibrational  level  which  is  coupled  with  the  first  excited 
nitrogen  vibrational  level. 

The  DO  range  ending  with  statement  160  computes  the 
rate  at  which  a s ngle  molecule  in  a given  level  is  ex- 
cited or  de -excited  by  collisions  at  each  altitude. 

The  outer  DO  range  ending  with  statement  180  computes 
wavelengths  (BNDLAM),  Einstein  A coefficients  (EA), 
rates  for  absorption  of  radiation  from  the  sun  (GDOWN) 
and  from  the  earth  and  lower  atmosphere  (GUP)  and  the 
G factor  used  in  computing  radiative  transfer  within  the 
upper  atmosphere  (GFAC).  The  inner  DO  range  ending 
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with  Statement  180  computes  the  unit  optical  thicknesses 
for  each  band  at  each  altitude  using  Function  TAUMAX, 
and  stores  the  result  in  array  BNDTAU.  In  addition, 
it  computes  the  volume  radiances  and  stores  them  in 
array  ENDRAD. 

At  this  point,  if  NLTE  is  1,  computation  for  thermal 
equilibrium  is  complete  and  control  is  transferred  to 
Statement  900.  If  NLTE  is  not  1,  the  computation  for 
non- equilibrium  is  performed  iteratively  in  the  DO 
range  ending  with  Statement  400.  The  DO  range  ending 
with  Statement  190  computes  various  integrated  optical 
..hicknesses  for  each  band.  The  DO  range  ending  with 
Statement  182  computes  the  total  thickness  between  each 
two  altitudes  in  the  vertical  direction  and  stores  the 
values  in  array  TTH. 

The  DO  range  ending  with  Statement  184  (1)  computes 
the  optical  thickness  between  the  base  altitude  and  each 
other  altitude,  storing  the  result  in  array  TTUP  and  (2) 
computes  the  total  optical  thickness  between  the  top 
altitude  and  each  other  altitude,  storing  the  results  in 
array  TTDOWN.  If  the  computation  is  not  for  daytime 
conditions,  the  DO  range  ending  with  Statement  186  sets 
the  values  in  TTDOWN  equal  to  0.  0,  for  efficiency  in  the 
later  computations.  If  the  computation  is  for  the  effects 
of  collisional  excitation  only,  the  DO  range  ending  with 
Statement  188  sets  the  values  in  TTUP  equal  to  0.0 

The  computation  of  vibrational  populations  now  begins. 
Control  is  transferred  to 

Statement  210  if  I is  1,  CH^ 

Statement  220  if  I is  2,  CO^ 
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Statement  230  if  I is  3,  FLO 

Statement  240  if  I is  4,  NO 

Statement  250  iflis5,  NO 

2 

Statement  230  if  I is  6,  pre-dawn  Ozone 

Statement  230  if  I is  7,  noon  Ozone 

Following  each  of  these  statements  is  a DO  range  in  which 

vibrational  populations  are  computed.  A discussion  of  two 

statements  in  the  computation  for  CO?  will  illustrate  the 

computational  details.  First,  look  at  the  statement  in  the 

DO  range  immediately  following  Statement  221  that  begins 

VIBPOP(NZ,  8)  - which  computes  the  population  of  the  first 

level  of  the  \)  vibrational  mode.  This  vibrational  level  is 

coupled  to  the  nitrogen  first  vibrational  level  by  collision, 

to  the  upward  flux  from  the  earth  and  lower  atmosphere, 

to  solar  flux,  and  to  lower  vibrational  levels  by  radiation. 

Excitation  by  absorption  of  radiation  fromwithin  the  upper 

atmosphere  is  neglected  for  this  level.  A ground  state 

molecule  is  excited  by  collision  with  nitrogen  at  the  rate 

VVNSIG(NZ)  (collision  exchange),  by  absorption  of  solar 

flux  at  a rate  GDOWN(9)*SMl(KB(9),  TDOWN2(NZ)) 

(equals  I^G'M^k^))  and  by  absorption  of  flux  from  the 

earth  and  lower  atmosphere  at  a rate  GUP(9)*SLM  1(KB(9), 

TTUP2(NZ))  (equals  Ni  G'N  (k  )).  An  excited  molecule 

'N,  t It 

is  de-excited  by  radiation  at  the  rate  EA(9)  (other  Einstein 
A coefficients  are  negligible),  by  vibrational  exchange 
with  nitrogen  at  a rate  VVRSIG(NZ),  and  by  transfer  of 
vibrational  energy  to  other  modes  within  the  C02  molecule 
at  a rate  VVDSIG(NZ).  The  vibrational  population  is  then 
the  ground  state  population  times  the  ratio  of  excitation 
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and  de  - excitation  rates.  Second,  look  at  the  statement 
beginning  POP2  = . . . , which  computes  the  vibrational 
population  of  level  2,  the  first  excited  level  of  the  0 

2 

vibrational  mode.  This  level  is  excited  from  the  ground 
state  by  collisional  excitation  at  a rate  TVS  = TVSIG, 
by  absorption  of  radiation  from  the  lower  atmosphere  and 
earth  at  a rate  GUP(1)*SLM1(KB(1),  TTUP(NZ,  1))  (equals 

NA,tG  Nl  ^kt^'  by  absorPtion  °f  radiation  from  other  points 

in  the  upper  atmosphere  in  the  15.  0 fundamental  band 

at  a rate  21F*GFAC(  1 )*(XFER(VIBPOP(  1 , 2),  TTH(1,1),  NZ, 

1,  1,  KB(1))  + XFER(VIBPOP(NZ,  2),  TTH(NZ,  1),  NALT, 

NZ,  2,  KB(  1 ) ) ) -'-EA(  1 ) ^equal s ^ ^G 'N^(k^)dz^  and  by 

absorption  of  solar  radiation  at  the  rate  GDOWN(l)*SMI 

(KB(  1 ),  TTDOWN(NZ,  1))  (equals  I G'M,(k  )).  It  is 

o Is 

excited  from  levels  3,  4 and  7 by  collisional  excitation 
VTS  and  by  radiation  EA(2),  EA(3)  and  EA(7).  De- 
excitation occurs  by  radiation  EA(1),  and  by  collisional 
de -excitation  to  the  ground  state  at  a rate  VTS  = VTSIG(NZ), 
and  by  collisional  excitation  to  higher  states  at  a rate  TVS. 
The  approximate  vibrational  population  is  then  the  ratio 
of  the  sums  of  population  lutes  times  higher  state  popu- 
lations divided  by  the  sum  of  de-excitation  rates  for 
single  molecule. 

After  each  set  of  vibrational  populations  is  completed,  new 
values  for  volume  band  radiances  are  computed  in  the 
nested  DO  ranges  ending  with  Statement  360. 

After  three  iterations  to  compute  vibrational  populations 
control  is  returned  to  the  main  program. 
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Inputs : 


NLEVEL,  NBAND,  NALT,  CDEXT,  CDEXB,  RM,  LSC, 
LCC,  I,  LBC,  LBU,  BY,  STR,  WAVE,  CONC,  SFLUX, 
TFLUX,  KCHEM,  NR,  NP,  NW,  NLTE,  TEMP,  BNDLAM, 
DENTOT,  NBAND1. 

Outputs:  VIBPOP,  BNDTAU,  BND.RAD 


Mnemonic  Variables: 


AA 

AAA 


BNDL 

BNDLAM 

BNDRAD 

BNDTAU 

BV 

CDEXB 

CDEXT 


CHEML 

CONC 

C02 


DENTOT 


intermediate  f ,ctor  in  computing  vibrational 
population  c 

quantity  used  in  NO  chemiluminescence  calculation 

in  original  Subroutine  BANRAD  application;  it  has 

no  effect  at  present 

wavelength  of  current  interest  (^m) 

wavelength  array  (/on) 

volume  radiance 

unit  optical  thickness 

rotational  constant  (cm  *) 

rate  coefficient  at  300  K for  vibrational  de- 

excitation  of  lowest  NO  vibrational  level  (cm  sec  ) 

rate  coefficient  at  300  K for  vibrational  exchange 

of  NO  with  N or  O for  the  most  significant 

^ 3 - 1 

vibrational  level  (cm  sec  ) 
currently  inactive 

_3 

array  of  NO  concentrations  (cm  ) 
real  array  used  to  store  reduced  temperature  value 
and  rate  coefficients  to  generate  vibrational  energy 
exchange  rates  for  CO^  currently  inactive 

array  of  total  number  densities  at  altitude  mesh 

, -3. 

(cm  ) 
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EA 

FAC 

FACTOR 

GDOWN 

GFAC 

GNDSTE 

GUP 

I 

ITER 

IX 

JF 

JI 

JL 

JX 

KB 

KBN 

KCHEM 

KL 

KLN 

LBC 

LBU 

LCC 


array  of  Einstein  A coefficients 
used  as  temporary  storage 

temporary  storage  in  computing  vibrational  exchange 

real  array  in  which  are  stored  rates  of  absorption 

of  solar  radiation,  I G'  for  each  band 

o 

real  array  in  which  are  stored  values  of  G1  for 
each  band 

real  array  in  which  are  stored  ground  state  vibrational 
populations,  equivalenced  to  VIBPOP(l,  1) 
real  array  in  which  are  stored  rates  of  absorption 
from  the  lower  atmosphere,  2TTN  v G'  for  each  band 

L 

index  of  species  of  interest  (equals  4 for  NO) 
integer  used  as  index  in  iterative  computation  of 
vibrational  level  population 
index  for  altitude  level 

index  for  final  vibrational  level  in  a band  transition 

index  for  initial  vibrational  1 evel  in  a band  transition 

integer  used  to  compute  JI,  JF 

index  for  vibrational  level 

integer  array  used  to  store  band  types 

index  of  band  type 

chemiluminescense  index;  currently  has  nc  effect 
integer  arr^y  storing  level  types 
current  level  type 

integer  array  containing  information  on  which 
levels  produce  a band 

integer  array  containing  band  type  information 
integer  array  containing  vibrational  level  coupling 
go  N2  or  02 
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LSC 

NALT 

NB 

NBAND 

NBAND 1 

NBC 

NBX 

NF 

NI 

NL 

NLEVEL 

NLTE 


NMIN 

NP 

NR 

NW 

NY,  NZ 

POP2,  POP3, 
POP4,  POP5, 
POP6,  POP7, 
POP8 

POWER 


RM 


integer  array  containing  level  type  information 
number  of  altitudes  for  which  results  are  to  be 
returned 

index  for  band  number 
number  of  bands  for  this  species 
equal  to  NBAND 
temporary  storage  index 
temporary  r.uorage  index 
temporary  storage  index 
temporary  storage  index 
level  number  index 

number  of  vibrational  levels  for  this  species 
integer  variable  used  to  indicate  conditions  of 
computation,  1 for  local  thermodynamic  equilibrium, 
2 for  collisional  excitation  only,  3 for  night 
conditions,  4 for  noon  conditions,  5 to  include 
high  vibrationax  temperature  for  nitrogen 
integer  variable  storing  (NALT  - 1) 
output  device  index 
input  device  index 
output  device  index 

index  for  altitude  level 
real  variables  used  as  temperature 

storage  for  vibrational  populations 

real  array  in  which  are  stored  conversion 

factors  between  emission  rates  in  photons 

-3  -1  j -3  -1 

cm  sec  and  watts  cm  ster  for 

each  band 

molecular  mass  for  this  molecule  (g/mole) 
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SFLUX 

SNB 

STR 

SUM 

TDOWN2 

TEMP 

TFLUX 

TMP 

TNNV 

TTDOWN 

TTH 

TTUP 

TTUP2 

TVS 

TVSIG 

VIBPOP 

VTS 

VTSIG 

VVDSIG 

VVNSIG 

VVRSIG 


array  of  solar  flux  values  for  each  band 

temporary  storage  of  band  strength  (cm.  \atm-cm)  *) 

array  of  band  strength  for  each  band  (cm  \atm-cm)  *) 

used  to  store  partial  sums 

real  array  equivalenced  to  TTDOWN(l,  9) 

array  containing  temperature  profile  (K) 

array  containing  fluxes  for  lower  altitudes  for 

each  band 

temporary  storage  for  temperature  (K) 

nitrogen  vibrational  temperature  (K) 

array  storing  values  of  total  optical  thickness 

from  top  of  atmosphere  for  each  band 

array  storing  values  of  total  optical  thickness 

between  adjacent  altitudes  fo^  each  band 

array  storing  values  of  total  optical  thickness 

from  bottom  of  upper  atmosphere  for  each  band 

array  equivalenced  to  TTITP(9,  1) 

used  to  store  local  value  of  TVSIG 

array  storing  values  of  vibrational  populations 

of  each  state  for  each  altitude 

array  storing  values  of  vibrational  populations 

of  each  state  for  each  altitude 

local  value  of  VTSIG 

array  of  vibrational  de- excitation  rates  at 
each  altitude 

array  of  estimated  values  of  transfer  of 
vibrational  energy  (downward) 
array  of  rales  of  transfer  of  vibrational  energy 
from  to  this  molecule 

array  of  rates  of  transfer  of  vibrational  energy 
to  N2 
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VVUSIG 


VVWAVE 

WATVA 


WAVE 

XYM 

XYMM 

XYN 


array  of  estimated  values  of  rate  of  transfer 

of  vibrational  exchange  (upward) 

difference  between  vibrational  energy  and  a 

molecular  vibrational  energy 

array  storing  reduced  temperature  values  a' d 

rate  coefficients  to  generate  vibrational  energy 

exchange  rates  for  CO^ 

array  containing  energy  level  of  each  band  (cm*1) 
temporary  storage  variables 
temporary  storage  variable 
temporary  storage  variable 


Non- Library  Externals: 

EVAL,  SIMP,  SLM1,  SMI,  TAUMAX,  XFER 


Called  By: 


BCKGND 
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Function  XFER 


Function  XFER  is  used  to  compute  the  excitation 
of  a vibrational  level  caused  by  molecules 
radiating  in  a given  band. 

Function  XFER  computes  the  contribution  to  the 
source  function  from  the  radiation  received  in  a 
given  vibration- rotation  band  by  a molecule  at  a 
given  altitude.  A single  entry  to  Function  XFER 
computes  thib  number  for  radiation  received 
from  molecules  located  either  below  or  above 
this  altitude.  Equation  (2-60)  is  summed  over 
each  altitude  interval  involved. 

Program  Flow:  Function  XFER  assumes  that  the  values  in  RAD(l) 

and  TTH(l)  are  the  values  corresponding  to  the 
lowest  altitude  over  which  the  integration  is  to 
be  performed.  Thus  when  called  from  Subroutine 
BANRAD  to  determine  the  effect  of  molecules 
above  a given  altitude  with  index  NZ,  the  calling 
argument  list  has  the  form  (VIBPOP(NZ,  2), 
TTH(NZ,  1),  NALT,  NZ,  2,KB(1))  for  a band 
transition  between  vibrational  levels  1 and  2 
and  band  number  1. 

Upon  entry  to  Function  XFER,  the  number  of 
altitude  intervals  is  found  from  the  difference  of 
NTOP  and  NBOT,  and  stored  in  NINT.  If  NTOP  = 
NBOT,  the  value  0.  0 is  returned  by  Function  XFER. 


Purpose: 


Method: 
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Otherwise,  K is  tested  and  control  transferred  to 
the  appropriate  part  of  the  function. 

When  K is  one,  XFER  computes  the  radiative 
contribution  to  excitation  from  layers  below  the 
altitude  level  for  which  the  computation  is  made. 
When  K is  two,  the  contribution  from  layers  above 
is  computed. 

The  computation  from  a single  level  is  found  in  each 
passage  through  the  loop  ending  with  statement  120 
(if  K is  one)  or  statement  220  (if  K is  two).  Before 
entry  to  either  loop,  RA  is  set  to  one,  TB  to  zero, 
and  EA3  to  zero.  RFAC,  the  ratio  of  ground  state 
to  upper  state  vibrational  populations  is  computed 
and  SUM  is  set  to  one-half.  Within  the  appropriate 
loop,  TB  is  incremented  by  the  optical  thickness 
obtained  from  array  TTH.  1/2  N (t)  for  the 

X w 

band  is  obtained  from  function  SLM1  and  the  value 

stored  in  EB2.  The  band  function  l/2L3(f)  is 

computed  from  the  value  of  L (X)  obtained  by  a 

o ' 

call  to  function  SL3  and  1/2N12(C)  through  use 
of  the  recurrence  relation 

2 L3(T)  = 4 tLo{C)  +tN12CC)  1 

Equation  (2-60)  is  then  applied  and  the  result 
added  to  the  value  in  SUM. 
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Mnemonic  Variables: 


EA3 


EB2 


EB3 


K 

KB 

NBOT 

NTOP 

NF 

NINT 

NPL 

NY 

NZ 

RA 

RAD 

RB 

RFAC 

SUM 


Temporary  storage  for  value  of  the  band  function 
Lg  for  the  nearer  of  two  adjacent  integration 
points 

Temporary  storage  for  value  of  the  band  function 
1/2  for  the  farther  away  of  two  adjacent 
integration  points 

Temporary  storage  for  value  of  the  band  function 
f°r  the  farther  away  of  two  adjacent  integration 
points 

integer  parameter  used  to  choose  appror-'ate 
section  of  this  function 

integer  parameter  used  to  choose  appropriate 
sections  of  functions  SL.M1  and  SL3 
indices  in  altitude  array  of  top  and  bottom 
altitudes  in  layer  of  interest 

index  used  to  select  values  from  array  VIBPOP 
(NTOP  --  NBOT),  number  of  intervals 
(NINT  + 1) 

index  for  altitude  level 
index  for  altitude  level 

relative  vibrational  population  of  nearer  altitude 
point 

real  array  containing  vibrational  population  of  an 
upper  state 

relative  vibrational  population  of  farther  altitude 
point 

ratio  of  ground  state  population  to  upper  vibrational 
level  population 

immediate  result  of  integration 
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TB 

TM 

TTH 

VIBPOP 

Non-Library 

Externals: 

Called  by: 


optical  distance  from  level  for  which  calculation 
is  being  made  to  farther  of  two  levels 
increment  of  optical  thickness  between  two  adjacent 
levels 

real  array  storing  incremental  optical  thicknesses 
real  array  storing  vibrational  populations  of  gromd 
and  excited  levels  of  a molecule 


SLM1,  SL3 
BANRAD 
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FUMCTTON  XFER  (R AO, TTH, NTOP,NBOT , K, KB ) 
COMMON  /POPCOM/  VI BPOP ( 201 , 25 > 
DIMENSION  RAO (201) , TTH ( 20 1 ) ,H OLD { 20 1 ) 
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Function  SL3 


Purpose: 


Method: 


Program 


Inputs: 

Output: 

Mnemonic 

K 


SL3 


This  function  yields  the  total  absorption  of  radiation 
from  a continuum  source  passing  through  a line 
absorbing  medium.  The  result  of  this  computation 
is  used  by  function  XFER  to  compute  the  band 
function  L^(C). 

SL3  computes  the  function  L (r)  for  parallel  and 
perpendicular  transitions  in  linear  molecules  The 
algorithm  uses  a polynomial  approximation  to  values 
of  L (r)  computed  for  a reduced  rotational  tempera- 
ture hcB^/kT  = 0.  002.  The  values  computed  by 
SL3  are  within  five  percent  of  the  correct  values 
for  any  reduced  rotational  temperature  encountered 
in  running  the  program. 

Qow:  Upon  entry,  the  absolute  value  of  the  argument  TAU 

is  defined  as  T.  Some  logic  then  is  used  to  choose 
the  appropriate  polynomial  curve  fit  to  the  function 
which  is  computed  and  return  is  executed. 

K,  TAU 

SLM1 

Variables: 

input  parameter  used  to  choose  appropriate  branch 
of  the  function,  K equals  1 for  parallel  band,  K 
equals  2 for  perpendicular  band, 
computed  value  of  L^(c). 
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* ^,Vj- 


-1— 


T 

TAU 

X 

Y 

Non-Library 

Externals: 

Called  by: 


absolute  value  of  argument  of  TAU 
input  variable,  optical  depth  kfc. 
argument  used  to  compute  curve  fit 
natural  logarithm  of  T. 


None 


XFER 
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FUMCTTON  SL  3 (K,TAU) 
= 5B9  (TAU) 
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Function  TAUMAX 


Purpose: 


Method: 


Program 


The  purpose  of  this  function  is  to  compute  the  optical 
thickness  of  the  strongest  line  in  an  infrared  vibrational  - 
rotational  band. 


Function  TAUMAX  computes  the  optical  thickness  kQ  at  the 
center  of  the  strongest  line  of  a vibrational- rotational  band 
of  a linear  molecule,  assuming  Doppler  line  shape. 


Flow:  The  optical  thickness  at  the  center  of  a line  with  Doppler 

shape  is 


k = 5.  3612  x 10' 
o 


18 


I M 

V T 


n,  x s;ine 

lc  o lu 


where  M is  the  molecular  mass  (gm/mole),  T is  the 

absolute  temperature  (k),  N,  is  the  columnar  number 

-2 

density  (molecules /cm  ),  is  the  line  center  wavelength 
(^m)  and  sj^ne  is  the  line  strength  in  absorption  (cm  ^ 
(atm-cm)  *)  at  300K.  Function  TAUMAX  computes  the 
reduced  rotational  temperature  hcZWkT  and  transmits 
this  as  an  argument  to  Function  SJMAX  which  computes 
the  fraction  of  total  band  strength  in  the  strongest  line 


SJ  1 Slu 
max 


band 


TAUMAX  then  computes  the  optical  thickness 


at  the  center  of  the  strongest  line,  using  S. 


S;;ne  in  the  above  expression. 


in  place  of 


max 


Inputs : 

KL,  KB,  RM,  BV,  T,  FLAM,  S,  DEN 

Output: 

TAUMAX 

Mnemonic  Variables: 

BV 

rotational  constant  (cm  *) 

DEN 

columnar  number  density  of  molecule 

FLAM 

wavelength  (jpn) 

KB 

integer  parameter 

KL 

integer  parameter 

RM 

mass  of  molecule  (gm/mole) 

S 

band  strength  (cm'^atm-cm)'1) 

SIG 

the  quantity  hcB  /kT 

V 

T 

temperature  (K) 

TAUMAX  value  of  optical  thickness  returned  by 

Non- Library 

Externals: 

SJMAX 

Called  By: 

BANRAD 

-2. 
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TAUMAX  = (5.  3613  x 10 


*FLAM*DEN*S*SJMAX(KL,  KB,  SIG) 


18  lRM 

N T 


O RETURN 


: 
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FU'lr  TT  IN  TAUMAX  ( KL  , K0 , RM , BV,  T , F L AM  , S , DE'J ) 
SIS  = 1.43879*8V/T 


V 7 

«T  CL 

i r 
*-  o 

« UJ  7 
t—  CL  U 
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Function  SJMAX 


Purpose: 


The  purpose  of  this  function  is  to  compute  relative  line 
strengths. 


Method: 


Function  SJMAX  computes  the  relative  strength  of  the 
strongest  line  in  a parallel  or  perpendicular  band  of  a 
linear  molecule,  assuming  that  all  parallel  bands  are 
^51  - ^ transitions  and  all  perpendicular  bands  are 

1 |Tf  1 5^"* 

7f  - 2.  transitions.  The  effect  of  neglecting  the  actual 

transition  involved  is  negligible  for  the  uses  of  Function 
SJMAX. 


Program  Flow:  Function  SJMAX  computes  the  J values  of  the  strongest 

line  in  the  band,  computes  the  strength  from  the  appropriate 
Honl- London  formula,  and  divides  this  by  the  rotational 
partition  function  returned  by  Function  QJPART. 

Inputs:  KL,  KB,  SIG 

Output:  SJMAX 


Mnemonic  Variables: 


A intermediate  result  in  calculation 

KB  parameter  to  choose  type  of  transition 

(KB  = 1^  - 51  transition)  (KB  = 2 77  -5l  transition) 
KL  parameter  passed  to  QJPART 
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SIG 

the  quantity  (hcB^/kT) 

SJM 

strength  of  strongest  line  in  band 

SJMAX 

relative  strength  of  strongest  line  in  band 

TJM 

effective  rotational  quantum  number 

Non-Librarv 

Externals: 

QJPART 

Called  By: 

TAUMAX 
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FLOW  CHART: 


SJMAX 
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Function  QJPART 


Purpose: 


Method: 


Program 

Flow: 


The  purpose  of  this  function  is  to  compute  rotational 
partition  functions  for  linear  molecules.  Refer  to 
Subroutine  BANRAD  for  definition  of  symbols. 

Function  QJPART  computes  the  rotational  partition  function 
of  a single  vibrational  level  of  a rigid  linear  molecule, 
assuming  that  the  energy  Ej  of  a rotational  state  with 
quantum  number  J is 

E = hcB  J(J  + 1) 

J v 

Provision  is  mu.de  for  computing  the  rotational  partition 
functions  for  vibrational  levels  in  which  rotational  states 
with  only  even  or  odd  J values  are  present. 

T rotational  levels  are  assumed  populated. 

If  K = 2,  only  even  rotational  levels  are  assumed 
populated.  If  K = 3,  only  odd  rotational  level  are 
assumed  populated. 


The  computation  uses  polynomial  approximates  to  the 
partition  function 

Qj  =^j(2J  + l)  exp  ( - Ej/kT) 

with  SIG  as  the  parameter.  The  range  of  applicability 
is  limited  to  values  of  SIG  less  than  0.  2,  and  the  largest 
value  of  SIG  encountered  for  the  molecules  presently 
included  is  0.  015. 
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K,  SIG 


Inputs : 
Outputs: 


OJPART 


Menemonic  Variables: 


K 

QJPART 

SIG 

X 


parameter  used  to  choose  appropriate  fit 
value  of  rotational  partition  function 


the  quantity  (hcB^/kT) 


the  value  of  SIG,  defined  for  convenience  in 
the  curve  fit 


Non-Library 
Externals : 


None 


Called  By; 


SJMAX 


| 

'I 

j| 
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Function  SIMP 


Purpose: 


The  purpose  of  this  function  is  to  compute  an 
integral  using  Simpson's  rule. 


Method: 


Function  SIMP  computes  the  integral  of  a function 
between  any  two  of  three  equally  spaced  points. 
Simpson's  rule  integration  is  used. 


Program 

Flow: 


Upon  entry,  the  value  of  K is  tested  to  choose  the 
appropriate  branch  of  the  function  with  the  following 


Inputs : 
Output: 


results : 

K < 2 

SIMP  = (h  / 12)  (8A  (2)  + 5 A (1)  - A (3) 

00 

II 

SIMP  = (h/  12)  (8A  (2)  + 5A  (3)  - A (1) 

K > 2 

SIMP  = (h/3)  (4A  (2)  + A (1)  + A (3)) 

A,  H,  K 

SIMP 


Mnemonic  Variables: 


A 

H 

K 


array  of  values  to  be  integrated 

integration  step  size 

parameter  to  choose  method  desired 


Non-Library 

Externals: 


None 


Called  By : 


BANRAD,  XFER 
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Function  SMI 


Purpose : This  function  yields  a factor  needed  to  compute  the  true 

rate  of  absorption  at  an  optical  depth,  k^,  of  radiation 
(from  a spectrally  constant  source)  in  a resonant  medium. 
Refer  to  Subroutine  BANRAD  for  further  explanation  and 
definition  of  symbols  used  below. 

Method:  SMI  computes  the  function  M (k  ) for  parallel  and 

1 o 

perpendicular  transitions  in  linear  molecules.  The 
algorithm  employs  a polynomial  approximation  to 
values  of  M^  (k^)  computed  at  a reduced  rotational 
temperature  hcB  /kT  = 0.  002. 


Program 

Flow-  Upon  entry  the  absolute  value  of  the  argument  TAU  is 

defined  as  T.  Some  logic  then  is  used  to  choose  the 
appropriate  polynomial  curve  fit  to  the  function  which 
is  computed  and  return  is  executed. 


Inputs:  K,  TAU 


Output:  SMI 


Mnemonic  Variables: 


K 


SMI 

T 


input  parameter  used  to  choose  appropriate  branch 

of  the  function;  K equals  1 for  parallel  band, 

K equals  2 for  perpendicular  band 

rate  of  absorption  at  optical  thickness  k , relative 

o 

to  zero  optical  thickness 

absolute  value  of  argument  TAU 

input  variable  - optical  depth  k 
250  ° 


TAU 


Mnemonic  Variables  (Cont'd) 


X argument  used  to  compute  curve  fit 

Y natural  logarithm  of  T 

Non-  Library 

Externals:  None 

Called  By;  BANRAD 
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FLOW  CHART: 


SMI 
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Function  SLM1 


Purpose: 


Method: 


Program 

Flow: 


Inputs; 

Output: 


This  function  yields  the  factor  used  to  obtain  the  actual 
absorption  of  radiation  from  an  extended  black  body 
surface  passing  through  a resonant  medium.  Refer 
to  Subroutine  BANRAD  for  further  explanation  and 
definition  of  symbols  used  below. 

SLM 1 computes  the  function  (kfc)  for  parallel 

and  perpendicular  transitions  in  linear  molecules.  The 
algorithm  uses  a polynomial  approximation  to  values  of 

2 7 r NI  (V  comPuted  for  a reduced  rotational  temperature 
hcB^/kr  = 0.  002.  The  values  computed  by  SLM  1 are 
within  five  percent  of  the  correct  values  for  any  reduced 
rotational  [temperature  encountered  in  running  the 
program. 

l 

1 

Upon  entry  the  absolute  value  of  the  argument  TAU  is 
defined  as  T.  Some  logic  then  is  used  to  choose  the 
appropriate  polynomial  curve  fit  to  the  function  which 
is  computed  and  return  is  executed. 

K,  TAU 


SLM  1 


Mnemonic  Variables: 


K 


SLM1 

T 

TAU 

X 

X 


ir.put  parameter  used  to  chocse  appropriate 
branch  of  the  function,  K equals  1 for  parallel 


band,  K equals  2 for  perpendicular  band 
computed  value  of  ^ 


27r  Nib<kt> 


absolute  value  A argument  TAU 
input  variable,  optical  depth  k 
argument  used  to  compute  curve  fit 
natural  logarithm  of  T 


Non-Library 

Externals: 


None 


Called  By: 


BANRAD 


FUNCTION  SLMi  (K » T AU ) 
T = AflS(TAU) 
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APPENDIX 


CHEMISTRY  PROGRAM  SNAPS2 

The  chemistry  program  developed  by  Degges  (1972)  ,s  listed  here 
for  reference  purposes.  Its  main  use  has  been  to  estimate  diurnal 
variations  in  chemical  abundances,  given  initial  molecular  species 
concentrations. 

The  reaction  rates  originally  used  are  listed  in  Table  A-l.  The 
new  values  of  Garvin  and  Hampson  (1974)  are  listed  in  Table  A-2  when 
the  values  vary  by  fifty  per  cent  or  more. 

The  chemistry  program  consists  of  a main  program,  SNAPS2, 
and  6 subroutines;  ROMTST,  TKEL,  SOLLY,  CONCEN,  COLUMN, 
and  SUNPHI.  The  main  program  performs  the  chemical  computations, 
including  the  effects  of  vertical  transport.  ROMTST  embodies  the  ex- 
trapolation method  of  Bulirsch  and  Stoer  (1966).  TKEL  is  used  to  com- 
pute temperatures  at  each  altitude.  SOLLY  computes  the  solar  zenith 
angle  and  fraction  of  the  solar  disk  visible  at  each  altitude.  COLUMN 
integrates  the  number  densities  of  molecular  nitrogen,  oxygen  and 
ozone  above  a given  altitude  in  the  direction  of  the  sun.  CONCEN  is 
called  by  COLUMN  to  obtain  number  densities  required  for  the  integration. 
SUNPHI  computes  the  attenuation  of  solar  flux  by  absorption  between  a 
given  altitude  and  the  assumed  top  of  the  atmosphere  in  the  solar  direc- 
tion and  returns  the  photodissociation  rates  required  for  the  chemistry. 

The  operation  of  the  main  program  is  described  briefly  below. 
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Table  A- 2 


Revised  Values  for  Reaction  Rates 
From  Garvin  and  Hampson  (1974) 


R5a 

r-H 

rH 

1 

o 

rH 

K 

rg 

• 

R6 

8.  0 x 10"  11  exp  (-500/T) 

R8 

- 3 3 

6.  7 x 10  exp  (+290/T) 

R9 

1.  6 x 10“12  exp  (-1000/T) 

RIO 

1.  0 x 10"  13  exp  (-  1250/T) 

R12 

should  vary  between  2x10 

R13 

4.  2 x 10-11  exp  (-350/T) 

R14 

3.  0 x 10“14  exp  (-4480/T) 

R15 

3.  0 x 10-11  exp  (-500/T) 

R16 

1.  7 x 10-11  exp  (-910/T) 

R17 

2.  8 x 10"12  exp  (-2125/T) 

R18 

„ „ - 12 

2.  8 x 10  exp  (-1900/T) 

R29 

2. 9 x 10" 16 

R30 

3.5x10  ^ 

R33 

- 13 

1.  4 x 10 

R36 

3.  0 x 10'12  exp  {-  1770/T) 

R47 

- 1 7 

6.5x10 

R48 

- 12 

9.  1 x 10 

R57 

- 13 

5. 7 x 10  1J 
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Photod  3sociation  Processes  Considered  for  Chemistry  Program 
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Program  Flow: 


Execution  of  Program  SNAPS  begins  with  the 
initialization  of  constants  used  by  the  program.  The 
values  of  eddy  diffusion  coefficients  used  at  altitudes 
above  40  km  are  stored  in  array  EDCOEF.  The  peak 
value  set  by  the  DATA  statement  occurs  at  104  km.  The 
eddy  diffusion  coefficient  profile  is  lowered  by  9 km  in  the 
DO  loop  ending  with  statement  90  and  the  low  altitude  eddy 
diffusion  coefficients  in  array  EDCF  are  transferred  to 
arrad  EDCOEF  in  the  DO  loop  ending  with  statement  95. 

This  permits  use  of  separate  sets  of  high  and  low  altitude 
diffusion  coefficients  by  changing  the  appropriate  data 
statements,  and  easy  variation  of  the  altitude  of  maximum 
turbulent  transport.  The  data  described  above  is  then 
read  in.  The  information  on  each  data  card  is  written  to 
printer  output  channel  6 as  it  is  read.  Columns  13  through 
80  of  the  second  through  fourteenth  data  cards  are  printed 
in  addition  to  the  number  placed  in  columns  2 through  12  so 
that  comments  may  be  added.  The  first  14  data  cards  are 
printed  on  a page  numbered  -1.  The  input  concentrations 
require  about  10  pages  to  print,  the  first  of  which  is  numbered 
page  0.  Further  execution  is  determined  by  the  input  data. 

A brief  flow  chart  is  given  on  the  next  page. 

The  outermost  control  DO  loop  ending  with  statement 
450  is  executed  IHRMAX  times.  Within  this  loop,  the 
computer  clock  is  called  and  computer  running  time  since 
the  beginning  of  the  program  is  stored  in  variable  SECA. 
Variable  locations  used  in  the  DO  loop  ending  with  state- 
ment 110  are  set  to  zero  and  this  loop  is  entered  and 
executed  for  each  altitude  between  NBASE  and  NTOP. 
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The  mean  molec  Tar  mass  RMWT  at  each  altitude  is 
computed  by  the  statements  associated  with  DO  loop  102. 
Temperatures  are  computed  by  calls  to  function  TKEL. 
The  graviational  acceleration  at  each  altitude  is  computed 

pj 

and  loial  mixing  scale  heights  av  are  computed. 
Coefficients  for  the  eddy  diffusion  contributions  to  g 
and  f]  terms  in  Equation  (3-4)  are  determined  from 
a finite  difference  approximation,  and  stored  in  arrays 
BTED,  QJEDL  and  QJEDU.  Factors  needed  for  compu- 
tation of  molecular  diffusion  contributions  are  computed 
and  stored  in  arrays  AFAC,  BFAC,  CFAC,  DFAC  and 
RDTEMP.  The  molecular  concentrations  at  each  altitude 
above  the  base  altitude  are  scaled  to  be  consistent  with 
the  finite  difference  scheme  in  the  DO  loop  ending  with 
statement  104.  An  effective  molecular  scale  temperature 
is  used  for  the  local  scale  height.  Before  exit  from  DO 
loop  110,  and  O^  concentrations  and  sample  computed 
values  are  written  to  the  printer  output  channel  for 
altitudes  between  but  not  including  those  associated 
with  NBASE  and  NTOP. 

The  DO  loop  ending  with  statement  400  is  entered 
and  executed  MMAX  times.  Each  execution  advances 
the  integration  time  by  DLTIME  seconds.  The  time  at 
the  end  of  the  step  is  computed  in  floating  and  fixed  point 
form.  If  the  value  of  NPRNT  permits,  a new  page  is 
started  and  the  time  and  column  headings  are  printed. 
Molecular  diffusion  coefficients  are  computed  for  the 
two  lowest  altitudes  i.  the  statements  associated  with 


The  DO  loop  ending  with  statement  350  is  entered 
and  executed  for  each  altitude.  Molecular  diffusion 
coefficients  are  computed  for  the  next  higher  altitude, 
and  the  molecular  diffusion  contributions  to  the  J3  and 
g terms  are  computed  and  combined  with  the  eddy 
diffusion  contributions  in  the  DO  loop  ending  with  state- 
ment 140.  The  solar  zenith  angle  CHI  and  the  fraction 
ETA  of  the  sun's  disc  visible  are  computed  by  Sub- 
routine SOLLY.  If  the  sun  is  visible,  column  counts 
of  N^,  O^,  O and  O^  in  the  solar  direction  are  computed 
in  Subroutine  COLUMN  and  photodissociation  rates  are 
computed  in  Subroutine  SUNPHI  fc the  time  at  the  mid- 
point of  the  time  interval  of  integration.  Reaction  rate 
coefficients  are  computed  for  the  temperature  at  the 
given  altitude. 

The  integration  for  a global  time  increment  DLTIME 
at  a single  altitude  is  performed  in  the  DO  lo.»p  ending 
with  statement  300.  The  operations  within  this  loop  are 
repeated  either  until  convergence  to  relative  accuracy 
EPS  is  achieved  or  until  MLMAX  repetitions  have  been 
performed  with  no  convergence.  Each  time  through  this 
loop,  a new  number  of  local  time  intervals,  MMF,  is 
computed  to  conform  to  the  extrapolation  procedure  of 
Subroutine  ROMTST.  Values  of  the  inital  concentrations 
at  the  integration  altitude  in  array  CONSP  are  placed  in 
array  DCX.  The  DO  loop  ending  with  statement  290  then 
integrates  over  the  MMF  local  time  steps.  The  values  in 
array  DCX  are  copied  into  array  TCX  and  the  time  in- 
crement DT  IME  is  set  to  one-half  the  length  of  the  local 
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increment.  An  index  ITER  is  set  to  one  and  appropriate 
third  body  concentrations  are  computed.  The  algorithm 
of  Equation  2-10  is  then  applied  for  each  species  in  the 
DO  loop  ending  with  statement  255,  the  results  being 
placed  in  array  RCX.  ITER  is  incremented  and  the 
results  in  array  RCX  are  placed  in  array  TCX.  The 
time  increment  is  doubled  and  the  DO  loop  ending  with 
statement  255  is  again  executed,  using  the  intermediate 
results  from  the  first  execution  to  compute  the  <=<  ’ J3 
and  g terms.  The  results  of  this  second  cycle, 
computed  in  array  RCX,  are  transferred  to  array  DCX 
and  the  statements  in  the  DO  loop  ending  with  statement 
290  repeated  until  the  global  time  interval  DLTIME  has 
been  covered.  Subroutine  ROMTST  is  called  to  test  the 
values  in  array  DCX  with  previous  results  produced  in 
the  DO  range  ending  with  statement  290.  The  improved 
concentration  values  are  returned  in  array  TCX.  If  con- 
vergence is  achieved,  control  is  transferred  to  statement 
310.  Otherwise,  the  DO  loop  ending  with  statement  300 
is  repeated.  If  this  process  does  not  result  in  convergence 
in  MLMAX  tries,  the  end  point  values  in  array  RCX  are 
placed  in  array  TCX  in  the  DO  loop  ending  with  statement 
315.  The  endpoint  results  for  the  next  lower  altitude, 
which  have  been  stored  in  array  SCX,  are  written  into 
array  CONSP  and  the  contents  of  array  TCX  are  placed 
in  array  SCX.  Production  rates  for  OH  by  two  different 
reactions  are  computed  and  stored  in  QCOHA  and  QCOHB. 


I 


"Rectangular"  sums  of  total  column  counts  are  incremented 


in  DZOHA  and  DZOHB.  These  values  are  computed  to 
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provide  a rough  check  on  the  hydrogen  compounc.  chemistry 
and  may  be  compared  with  airglow  observations.  If  the 
value  of  NPRNT  permits,  the  altitude,  QCOHA,  QCOHB, 
DZOHA,  DZOHB  and  the  concentrations  of  16  molecules 
are  written  to  the  printer  output  channel.  The  statements 
in  the  DO  loop  ending  with  statement  350  are  then  repeated 
until  concentrations  at  all  altitudes  between  (but  not  in- 
cluding) NBASE  and  NTOP  have  been  integrated  over  the 
global  time  step  DLTIME. 

When  all  altitudes  have  been  integrated,  a summary  of 
number  of  iterations  is  sent  to  the  printer  output  channel 
and  the  time  value  TIME  is  incremented  by  DLTIME. 

After  the  integration  has  been  advanced  MMAX  times 
with  time  increment  DLTIME,  all  concentrations  are  sent 
to  the  punch  output  channel  by  the  statements  in  DO  loop 
420  and  to  the  printer  channel  by  the  statements  in  DO  loop 
440.  The  internal  computer  clock  is  again  accessed  and 
the  starting,  end,  and  elapsed  time  required  in  DO  loop 
450  are  sent  to  the  printer  output  channel.  Program 
execution  ceases  after  this  outer  DO  loop  has  been 
executed  IHRMAX  times. 

The  scaling  of  concentrations  performed  in  the  DO 
loop  ending  with  statement  104  permits  use  of  output  data 
from  runs  at  another  latitude  or  with  a different  temperature 
profile  to  be  used  as  input  data.  While  there  may  be  some 
difficulty  in  reaching  convergence  at  altitudes  below  90  km 
in  the  first  time  step,  this  is  not  a serious  problem.  When 
using  a global  time  step  DLTIME  of  600  seconds  and  a 
value  of  0.01  for  EPS,  computation  of  cnemical  changes 
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over  a day  for  150  altitudes  requires  on  the  order  of 
20  minutes  running  time  on  the  CDC  6600.  This  is 
longer  than  would  be  required  by  the  differencing  methods 
usually  employed  (e.g.,  Shimazaka  and  Lairo,  1970)  but 
the  error  control  compensates  for  this.  It  is  not  necessary 
to  alter  time  increments  at  dawn  or  sunset  since  the  con- 
vergence algorithm  of  Subroutine  ROMTST  automatically 
forces  the  use  of  smaller  time  subdivisions  of  the  global 
time  interval  DLTIME. 


The  necessary  input  data  to  the  chemistry  program 

is  a set  of  14  data  cards  containing  identification  and 

variable  parameters  and  a set  of  603  cards  containing 

3 

concentrations  of  N2,  O 2>  0(  P),  C>3,  OH,  H,  H O,  HC>2, 

H2°2’  H2’  CH4’  C°2'  CO’  °(1d)’  N<4s)’  no>  N02’  N2°’ 
N(  D)  and  Ar  at  one  km  intervals  between  ground  level  and 

200  km.  The  data  is  read  in  the  following  order: 

1.  IDENT  60  characters  of  identifying  comments 

which  are  used  in  output  page  headings. 

The  first  character  should  be  a blank. 

Co-latitude,  degrees,  of  the  computation 

location.  A value  of  zero  corresponds  to 

the  nort  h pole. 

East  longitude,  degrees,  of  computation 
location. 

Solar  declination  angle,  degrees,  used  to 
determine  time  of  year. 

5.  GREF  Surface  graviational  acceleration, 

/ 2 

cm/sec  , at  computation  location. 


2.  DTHETA 


3.  DPH1 


4.  SOLDEC 
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6.  REFF 

7.  EPS 

8.  MLMAX 

9.  NBASE 

10.  NTOF 

1 1 . IHRMAX 

12.  MMAX 

13.  DLTIME 


Effective  value  of  earth  radius,  km,  for 
use  in  computing  graviational  accelerations 
at  points  above  gror  ::d  level. 

Relative  error  permitted  for  each  time 
step  in  the  integration. 

Maximum  number  of  times  integration 
is  to  be  attempted  for  each  time  step  to 
obtain  relative  error  EPS. 

Index  of  base  altitude  for  computations, 
altitude  in  kiyi  plus  one. 

Index  of  highest  altitude  in  computations, 
altitude  in  km  plus  one,  restricted  to  201 
or  less  in  current  program. 

Number  of  times  at  which  complete  output 
is  required. 

Number  of  integration  steps  between  complete 
output  points. 

Intagration  interval,  seconds;  the  total 
simulated  time  for  a given  run,  in  seconds, 


is  IHRMAX  * MMAX  * DLTIME. 


14.  TIME  Initial  time,  seconds,  measured  from  noon; 

AM  values  are  negative,  PM  values  are 
positive. 

15.  CONSP  The  603  cards  containing  concentrations  and 

identifying  indices;  concentrations  are  read 
into  array  CONSP. 


Outputs:  A complete  set  of  data,  identical  to  the  input  data,  is 

written  on  tape  4 at  the  end  of  each  major  time  step. 
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Mnemonic  Variables: 


Real  Arrays;  AFAC 


AMWT 


BFAC 


BTED 


BTEMD 


CARD 

CFAC 


CMD 


CMDL 

CMDU 

COL 


Stores  values  of  numerical  approximation 

to  — -Al 
t dz  * 

Stores  values  of  molecular  weights, 

gm/mole;  values  are  set  by  data  state- 
ment. 


Stores  values  of  the  coefficient  of  molecular 
concentration  in  the  numerical  approximation 


to 


\ 2 
O n 


d 2 


r + 


_l  -At  n _ _i_  I 


Z " bz~  T* 

Stor«ss  values  of  coefficients  used  in 
computing  contribution  to  J]  terms  from 
eddy  diffusion. 


Stores  values  of  contribution  to  J3  term 
from  both  eddy  and  molecular  diffusion 
for  each  molecule  at  a given  altitude  during 
integration. 

Used  as  temporary  storage  for  comments  in 

columns  13  through  80  of  first  14  data  cards. 

Stores  values  of  numerical  approximation 

used  to  compute  for  diffusive 

H 

scale  heights. 


Stores  values  of  molecular  diffusion  co- 
efficients for  each  molecular  species  at  a 
given  altitude  during  integration. 

Stores  values  of  CMD  at  next  lower  altitude. 
Stores  values  of  CMD  at  next  higher  altitude. 
Stores  values  of  integrated  column  counts 
computed  by  Subroutine  COLUMN. 
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CONSP 

DCX 

DFAC 

EDCF 

EDCOEF 

QJEDL 

QJEDU 

QJEMD 

RATE 

RCX 


Stores  values  of  molecular  number 
densities;  initial  values  are  read  in  from 
cards;  values  are  altered  at  each  inte- 
gration step. 

Stores  initial  number  densities  for  one 
altitude  during  local  time  step  in  inte- 
gration. 

Stores  values  of  g/RT  used  in  computing 
atmospheric  scale  heights  and  contribution 
of  molecular  diffusion. 

Stores  values  of  eddy  diffusion  coefficients 
for  altitudes  below  40  km,  set  by  DATA 
statement. 

Stores  values  of  eddy  diffusion  coefficients, 
initially  contains  high  altitude  values  set 
by  DATA  statement. 

Stores  coefficients  of  eddy  diffusion  con- 
tribution to  g terms  from  next  lower 
altitude. 

Stores  coefficients  of  eddy  diffusion  contri- 
butions to  g terms  from  next  higher 
altitude. 

Stores  values  of  contribution  to  g term 
from  eddy  and  molecular  diffusion  for  each 
molecular  species  at  a given  altitude 
during  integration. 

Stores  values  of  photodissociation  rates 
computed  by  Subroutine  SUNPHI. 

Stores  intermediate  values  of  number  den- 
sities for  a single  aH’tude  computed  during 
local  time  step  in  integration. 
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RDTEMP 


sex 


TCX 


TEMPK 


Stores  values  of  temperature  coefficients 
used  in  computing  molecular  diffusion 
coefficients . 

Stores  values  of  number  densities  computed 
at  end  of  global  time  step  for  next  lower 
altitude;  values  eventually  replace 
corresponding  numbers  in  CONSP. 

Stores  intermediate  values  of  number 
densities  used  in  computing  e><  , J3  and 
g values  during  integration  of  local  time 
step  for  a single  altitude. 

Stores  values  of  temperatures  at  each 
altitude,  values  are  computed  by  function 


TETOT 


TFMT 


TINO 


TIO 


TI02 


TST 


1 


Stores  values  of  total  electron  densities 
at  each  altitude,  values  are  set  by  DATA 
statement. 

Stores  format  statement  for  computer 
running  time  printout,  set  by  DATA  state- 
ment. 


Stores  values  of  NO  concentrations  at  each 


altitude,  set  by  DATA  statement. 


Stores  values  of  concentrations  at  each 


altitude,  set  by  DATA  statement. 


Stores  values  of  concentrations  at  each 


altitude,  set  by  DATA  statement. 

Used  as  workspace  by  Subroutine  ROMTST. 


The  names  of  integer  arrays  and  their  uses  are: 

IDENT  Stores  60  character  comment  read  from 
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- - 


HiMUMM 


- ‘ - 


ini 


IDNT 


LSP 


MLPRNT 


Real  Variables: 


AFACNA 


BFACNA 


BTEDNA 


BTSQ 

BTT 


CFACNA 


CHI 


first  data  card,  which  is  used  in  page 
headings. 

Stores  image  of  last  20  characters  of 
first  data  card. 

Stores  Hollerith  constants  used  in  printing 
column  headings. 

Stores  data  on  number  of  iterations  re- 
quired for  convergence  of  integration  at 
a given  altitude,  and  (largest)  index  of 
species  which  did  not  converge  in  MLMAX 
iterations. 


Temporary  storage  of  value  of  AFAC  for 
altitude  at  which  integration  is  being 
performed . 

Temporary  storage  of  value  of  BFAC  for 
altitude  at  which  integration  is  being 
performed. 

Temporary  storage  of  value  of  BTED  for 
altitude  at  which  integration  is  being  per- 
formed . 

Squire  of  J2  term  used  in  integration. 
Value  (partial  or  complete)  of  JJ  term 
used  in  integration  for  a single  molecular 
specie  s . 

Temporary  storage  of  value  of  CFAC  for 
altitude  at  which  integration  is  being  per- 
formed. 

Solar  zenith  angle,  radiance,  computed  by 
Subroutine  SOLLY. 
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- - 


- — — ■ 


CHISAV 

CMDNSP 


CON 


CONFAC 


CONSPN 


Cl 


C2 


C3 


DCT 


DDDZ 


DELZ 

DFACNA 


DHDZ 

DKDZ 


Temporary  storage  for  value  of  CHJ. 
Temporary  storage  lor  value  of  mole- 
cular diffusion  coefficient  for  single 
molecular  species. 

Temporary  storage  for  factor  used  in 
scaling  molecular  concentrations  at 
altitude  NTOP  to  fit  scale  height  defined 
by  molecular  scale  temperature. 
Temporary  storage  for  factor  used  to 
scale  molecular  concentrations  to  fit 
local  scale  heights  defined  by  molecular 
scale  temperature. 

Temporary  storage  for  concentration  of  a 
single  molecular  species  at  a given  altitude. 
Temporary  storage  for  molecular  nitrogen 
concentration. 

Temporary  storage  for  molecular  oxygen 
concentration. 

Temporary  storage  for  atomic  oxygen  con- 
centration. 

Temporary  storage  for  single  value  from 
array  DCX. 

Temporary  storage  for  value  of  ~v  D for 

o z 

single  molecular  species. 

Altitude  increment,  cm. 

Temporary  storage  for  value  of  DFAC  at 
altitude  for  which  integration  is  being 
performed. 


Temporary  storage  for  value  of 

<3 


d H 


Temporary  storage  for  value  of  _a_K  . 

<5  z 
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DLT 

DLTIME 

DPHI 


DTDZ 
D THETA 

DTIME 

DZOHA 

DZOHB 

ED 

EPS 

ETA 

FAC 

FACB 

FACDMB 

FACM 

FACRT 

FAC2AX 

FAC4AQ 

G 


Value  of  0 used  in  integration  for  a single 
molecular  species. 

Time  increment  for  integration  step,  read 
in  as  data. 

East  longitude  of  point  of  integration,  read 
in  as  data;  not  needed  in  present  program 
which  assumes  times  are  true  solar  times. 
Temporary  storage  for  value  of 
Colatitude  of  point  at  which  integration  is 
performed,  read  in  as  data. 

Storage  for  value  of  local  time  step  in  in- 
tegration. 

Values  of  integrated  production  rates  of  OH. 

Temporary  storage  for  value  of  eddy 
diffusion  coefficient  at  a given  altitude. 

Value  of  relative  error  allowed  in  integration. 
Fraction  of  solar  disc  visible  at  given 
altitude. 

Temporary  storage  for  intermediate  results, 
use  varies  in  different  parts  of  program. 
Temporary  storage  for  either  J3 1 or  6t  in 
Equations  (2-10). 

Temporary  storage  for  <6  - /?  in  Equatio  (3-5). 

Temporary  storage  for  value  of  l-exp'-C^t). 

2 

Temporary  storage  for  4cXg /JJ  . 

Temporary  storage  for  2o<  x. 

Temporary  storage  for  4o<g. 

Value  of  graviational  acceleration  at  altitude 
NALT. 


GL 

GREF 

GU 

QCOHA 

QCOHB 

QJEDLN 

QJEDUN 

QJT 


RDELZ 

RE 

REFF 

RFAC 

RI1,  ... 

rmwt 

RMV/TL 

RMWTU 

ROBS 

RTIME 

R1 , ... 

R2DELZ 


Value  of  graviational  acceleration  at 
altitude  NALT  -1. 

Surface  value  of  graviational  acceleration, 
read  in  as  data. 

Value  of  graviational  acceleration  at 
altitude  NALT  +1. 

Values  of  OH  production  rates  at  altitude 
of  integration. 

Temporary  storage  for  value  of  QJEDL  at 
altitude  of  integration. 

Temporary  storage  for  value  of  QJEDU  at 
altitude  of  integration. 

Temporary  storage  for  partial  or  complete 
value  of  g term  used  in  integration  of  a 
single  molecular  species. 

Reciprocal  of  twice  the  altitude  increment 
DELZ. 

Effective  radius  of  earth,  km,  used  in 
geometrical  computations. 

Effective  radius  of  earth,  km,  used  in 
gravitational  calculations,  read  in  as  data. 
Temporary  storage  of  intermediate  results. 
Storage  for  values  of  ionic  reaction  rates. 
Mean  molecular  mass  at  altitude  NALT. 
Mean  molecular  mass  at  altitude  NALT  -1. 
Mean  molecular  mass  at  altitude  NALT  +1, 
Altitude  of  given  point,  km. 

Variable  used  in  printing  time  values. 

Values  of  forward  neutral  reaction  rates. 
Square  of  reciprocal  of  altitude  increment, 
DELZ. 


SEC  A 
SECB 
SECC 
SFAC 
SH 

SHL 

SHU 


SI1,  ... 
SOLDEC 
SUM 
SUML 
SUMU 
SI , ... 

TCH4,  ... 


TEMP 

TEMPI 

TEMPL 

TEMPU 

TEMP  12 

TEMP32 


Storage  for  values  of  computation  of 
computer  running  time. 

Temporary  storage  of  intermediate  results. 
Value  of  mixing  scale  height  at  altitude 
NALT. 

Value  of  mixing  scale  height  at  altitude 
NALT  -1. 

Value  of  mixing  scale  height  at  altitude 
NALT  +1. 

Values  of  reverse  ionic  reaction  rates. 

Solar  declination  angle,  read  in  as  data. 

Temporary  storage  used  in  forming  mean 
molecular  masses. 

Values  of  reverse  neutral  reaction  rate 
coefficients. 

Values  of  individual  molecular  concentrations 
in  array  TCX,  named  individually  for  con- 
venience in  forming  c<,  j]  and  g coefficients. 
TCH4  corresponds  to  CH^,  TCO  corresponds 
to  CO,  etc. 

Value  of  temperature  at  altitude  NALT. 
Reciprocal  of  temperature. 

Value  of  temperature  at  altitude  NALT  - 1. 
Value  of  temperature  at  altitude  NALT  +1. 
Value  of  minus  one-half  power  of  temperature. 
Value  of  minus  three-halves  power  of 
temperature. 
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TEMP52 


TFAC 

TIME 


XM4,  . . . 


XXX 12 


Value  of  minus  five-halves  power  of 
temperature. 

Temporary  storage  of  intermediate  results. 
Time,  expressed  in  seconds  from  noon, 
read  in  as  data;  altered  in  program  to  be 
the  time  at  midpoint  of  a global  integration 
step  so  that  approximate  mean  values  of 
photodissociation  rates  will  be  computed  by 
Subroutine  SUNPHI. 

Third  body  concentrations  for  appropriate 
three-body  reactions. 

Storage  for  factor  used  in  computing  g term 
for  CO  and  JJ  term  for  CO. 


Integer  Variables: 


IHRMAX 


ITER 


Index  in  DO  loop  95,  auxiliary  index  in  DO 
loop  110,  and  index  in  aimplied  input/out- 
put DO  loops. 

Final  parameter  in  innermost  implied  DO 
loop  used  in  printing  output  concentrations 
within  DO  loop  440. 

Index  in  DO  loop  450. 

Final  parameter  in  DO  loop  450,  input  as 
data. 

Index  in  print  control  DO  loop  440. 

Initial  parameter  in  implied  DO  loop,  see  IF. 
Control  index  in  integration,  used  between 
statements  200  and  260  to  select  proper 
values  for  array  TCX. 
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IX 


KCH4, 


KTEST 


Index  in  DO  loop  435,  incremented  for 
each  page  of  output. 

Index  in  implied  DO  loops  used  in  reading 
data  into  array  CONSP  within  DO  loop  100. 
Index  variables  used  to  select  appropriate 
values  of  quantities  related  to  given  species; 
most  uses  replaced  by  use  of  TCH4,  etc. 
Parameter  returned  by  Subroutine  ROMTST, 
zero  if  integration  has  converged,  non- 
zero otherwise. 

LALT,  LCHI,  LCOLN2,  LC0L02,  L0OLO3,  LDZOHA, 
LDZOHB,  LJ02,  LJ03,  LML,  LQCOHA,  LQCOHB,  LTIME 
Store  Hollerith  constants  used  in  printing 
column  headings  and  are  set  by  DATA  state- 
ments. 

Index  of  DO  loop  400,  incremented  for  each 
global  time  step. 

Index  of  DO  loop  300,  incremented  for 
each  iteration  of  integration  over  global 
time  step  at  a given  altitude. 

Final  parameter  of  DO  loop  300,  read  in  as 
data. 

Temporary  storage  for  value  of  ML,  used 
after  exit  from  DO  loop  300  to  print  number 
of  iterations  required  for  convergence. 

Final  parameter  of  DO  loop  400,  read  in  as 
data. 

Final  parameter  of  DO  loop  290,  computed 
from  ML  and  is  number  of  local  subdivisions 
of  global  time  step  DLTIME  for  the  ML 
iteration  of  integration. 


M 


ML 


MLMAX 


MLSAVE 


MMAX 


MMF 
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MMS 

NA 

NALT 

NAMIN 

NB 

NBASE 

NF 

NMIN 

NPAGE 

NPLU 

NPRNT 


NS 

NSP 


NSPM 


NT 


Index  for  DO  loop  290,  see  MMF. 
Auxiliary  index,  usually  set  to  NALT  -1. 
Index  for  DO  loops  90,  100,  110,  350, 

420  and  430;  denotes  a given  altitude. 
Value  of  NALT  -1  in  DO  loop  350. 

Value  of  NALT  -2  in  DO  loop  110. 

Index  of  lowest  altitude  for  integration, 
read  in  as  data. 

Final  parameter  of  DO  loop  430. 

Stores  input  value  of  NBASE,  which  is 
changed  by  program. 

Stores  output  page  number. 

Value  of  NALT  +1  in  DO  loop  350. 

Variable  available  to  control  output  of 
concentrations  at  integration  steps  inter- 
mediate to  fully  documented  output  at  end 
of  DO  loop  450,  set  to  1 at  present,  so 
values  are  output  for  each  interval 
DLTIME. 

Initial  parameter  of  DO  loop  430. 

Index  in  DO  loops  102,  104,  130,  140,  180, 
190,  255,  258,  280,  315,  and  320;  a value 
of  NSP  selects  a value  of  a variable  corre- 
sponding to  a single  molecular  species. 
Number  of  molecular  species  included  in 
program,  20  at  present,  used  as  final 
parameter  of  DO  loops  with  NSP  as  index. 
Final  parameter  of  DO  loop  420,  in  which 
output  cards  are  punched. 
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NTIME 


NTOP 

NW 

NXX 


N1 

N2 

N3 


Used  in  printing  time  values. 

Index  of  highest  altitude  used  in  inte- 
gration, input  as  data. 

Channel  number  for  printing  computer 
running  tim®  information,  set  to  6 here. 
Temporary  storage  for  altitude  index 
read  from  cards  containin,  concentration 
data. 

Storage  space  used  for  sequence  numbers 
of  cards  read  in  or  punched  out  at  a given 
altitude. 


Flow  Chart  for  Chemistry  Program 


^SNAPS^) 


SET  CONSTANTS 
READ  AND  WRITE  DATA 


COMPUTE  EDDY  DIFFUSION  CONTRIBUTIONS 
COMPUTE  FACTORS  FOR  MOLECULAR  DIFFUSION 
SCALE  CONCENTRATIONS 


COMPUTE  BASE  ALTITUDE  MOLE< 
M PLACE  BASE  CONCENTRATIONS 

r 

3ULAR  DIFFUSION  COEFFICIENTS 
IN  TEMPORARY  STORAGE,  SCX 

nsr 

i 

COMPUTE  DIFFUSION  COEFFICIENTS  FOR  GIVEN  ALTITUDE 
COMPUTE  MOLECULAR  DIFFUSION  CONTRIBUTIONS 
COMPUTE  SOLAR  ZENITH  ANGLE,  SUBROUTINE  SOLLY 
COMPUTE  COLUMN  COUNTS  IN  SOLAR  DIRECTION,  SUBROUTINE  COLUMN 
COMPUTE  PHOTODISSOCIATION  RATES,  SUBROUTINE  SUNPHI 
COMPUTE  RATE  COEFFICIENTS 


© 


SET  LOCAL  TIME  STEP 

PLACE  INITIAL  CONCENTRATIONS  IN  TEMPORARY  STORAGE,  DCX 


PLACE  CONCENTRATIONS  AT  START  OF  LOCAL  TIME 
STEP  IN  ARRAY  TCX 

SE1  TIME  INCREMENT  TO  1/2  LOCAL  TIME  STEP 
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Flew  Chart  for  Chemistry  Program  (continued) 


w 

^ COMPUTE  THIRD  BODY  CONCENTRATIONS 


COMPUTE  q TERMS  FOR  SINGLE  SPECIES 

APPLY  INTEGRATION  ALGORITHM,  RESULT  TO  ARRAY  RCX 


NO 


ALL  SPECIES 
INTEGRATED^ 


YES 


LOCAL  TIME  STEP" 
RESULT  AT 
ENDPOINT 
" — OBTAINED 


NO 


-YJSS 


SET  TIME  INCREMENT  TO  LOCAL  TIME  STEP 
PLACE  RESULTS  (ARRAY  RCX)  IN  ARRAY  TCX 


PLACE  RESULTS  (ARRAY  RCX)  IN  ARRAY  DCX  k~ 


GLOBAL  INTERVAL 
COMPLETED? 


YES 


apply  extrapolation,  SUBROUTINE  ROMTST 
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Flow  Chart  for  Chemistry  Program  (continued) 


Flow  Chart  for  Chemistry  Program  (continued) 
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